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ABSTRACT 


The Fermi Large Area Telescope {Fermi- LAT, hereafter LAT), the primary 
instrument on the Fermi Gamma-ray Space Telescope (Fermi) mission, is an 
imaging, wide field-of-view, high-energy 7 -ray telescope, covering the energy 
range from 20 MeV to more than 300 GeV. During the first years of the mission 
the LAT team has gained considerable insight into the in-flight performance of 
the instrument. Accordingly, we have updated the analysis used to reduce LAT 
data for public release as well as the Instrument Response Functions (IRFs), the 
description of the instrument performance provided for data analysis. In this 
paper we describe the effects that motivated these updates. Furthermore, we 
discuss how we originally derived IRFs from Monte Carlo simulations and later 
corrected those IRFs for discrepancies observed between flight and simulated 
data. We also give details of the validations performed using flight data and 
quantify the residual uncertainties in the IRFs. Finally, we describe techniques 
the LAT team has developed to propagate those uncertainties into estimates of 
the systematic errors on common measurements such as fluxes and spectra of 
astrophysical sources. 

Subject headings: instrumentation: detectors - instrumentation: miscellaneous - 
methods: data analysis - methods: observational - telescopes 
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1. INTRODUCTION 


The Fermi Gamma-ray Space Telescope (Fermi) was launched on 2008 June 11. Com- 


missioning of the Fermi Large Area Telescope (LAT) began on 2008 June 24 (Abdo et ah 
2009f). On 2008 August 4, the LAT began nominal science operations. Approximately 


one year later the LAT data were publicly released via the Fermi Science Support Center 
(FSSCfl 


The LAT is a pair-conversion telescope; individual 7 rays convert to e + e~ pairs, which 
are recorded by the instrument. By reconstructing the e + e~ pair we can deduce the energy 
and direction of the incident 7 ray. Accordingly, LAT data analysis is entirely event-based: 
we record and analyze each incident particle separately. 


In the first three years of LAT observations (from 2008 August 4 to 2011 August 4), 
the LAT read out over 1.8 x 10 11 individual events, of which ~ 3.4 x 10 10 were transmitted 
to the ground and subsequently analyzed in the LAT data processing pipeline at the LAT 
Instrument Science Operations Center (ISOC). Of those, ~ 1.44 x 10 8 passed detailed 7-ray 
selection criteria and entered the LAT public data set. 


The LAT team and the FSSC work together to develop, maintain and publicly distribute 
a suite of instrument-specific science analysis tools (hereafter ScienceTool^ | that can be 
used to perform standard astronomical analyses. A critical component of these tools is the 
parametrized representations of instrument performance: the instrument response functions 
(IRFs). I11 practice, the LAT team assumes that the IRFs can be factorized into three parts 
(the validity of this assumption is studied in § 7.1.4): 


1. Effective Area, A e s(E , v , s), the product of the cross-sectional geometrical collection area, 
7 ray conversion probability, and the efficiency of a given event selection (denoted by s ) for 
a 7 ray with energy E and direction v in the LAT frame; 

2. Point-spread, Function (PSF), P(v'\E,v,s ), the probability density to reconstruct an 
incident direction v' for a 7 ray with (E, v) in the event selection s; 

3. Energy Dispersion , D(E'] E,v, s), the probability density to measure an event energy E' 
for a 7 ray with (E, v ) in the event selection s. 

The IRFs described above are designed to be used in a maximum likelihood analyst^] as 


^http:/ /fermi. gsfc.nasa.gov/ssc 

2 ht t p : / /fermi . gsfc . nasa. gov / ssc /data/ analysis / software 

^http:/ /fermi. gsfc. nasa.gov/ssc/data/analysis/scitools/likelihood_tutorial.html 
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described in Mattox et al. (1996). Given a distribution of 7 rays S(E,p), where p refers 


to the celestial directions of the 7 rays, we can use the IRFs to predict the distribution of 
observed 7 rays M(E',p',s ): 


( 1 ) 


M(E',p',s)=J I J S(E,p)A eS (E,v(t;p),s)x 

P(v'(t,p '); E, v(t] p), s)D(E'; E, v(t;p ), s)dEdttdt. 


The integrals are over the time range of interest for the analysis, the solid angle in the LAT 
reference frame and the energy range of the LAT. 

Note that the IRFs can change markedly across the LAT field-of-view (FoV) . Therefore, 
we define the exposure at a given energy for any given energy and direction in the sky S(E,p) 
as the integral over the time range of interest of the effective area for that particular direction; 

£{E,p,s) = J A eS (E,v(t,p),s)dt. (2) 

Another important quantity is the distribution of observing time in the LAT reference frame 
of any given direction in the sky (henceforth referred to as the observing profile, and written 
Gbs), and which is closely related to the exposure: 

£(E,p,s) = J A eS (E,v,s)t ohs (v-p)dQ. (3) 


The absolute timing performance of the LAT has been described in detail in Abdo et~ah 


(2009f) and Smith et al. (2008) and will not be discussed in this paper. 


To allow users to perform most standard analyses with minimum effort, the LAT team 
also provides, via the FSSC, a spatial and spectral model of the Galactic diffuse 7-ray 
emission and a spectral template for isotropic 7-ray emission)^] In this prescription, con- 
tamination of the 7-ray sample from residual charged cosmic rays (CR) is included in the 
isotropic spectral template. Although not part of the IRFs, this background contamination 
is an important aspect of the instrument performance. 


From the instrument design to the high level source analysis, the LAT team has relied 
heavily on Monte Carlo (MG) simulations of 7-ray interactions with the LAT to characterize 
performance and develop IRFs. The high quality data produced since launch have largely 
validated this choice. However, unsurprisingly, the real flight data exhibited unanticipated 


4 


http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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features that required modifications to the IRFs. After years of observations the LAT data set 
itself is by far the best source of calibration data available to characterize these modifications. 

LAT event analysis has substantially improved since launch. We have applied the accu- 
mulated updates in occasional re-releases of data, corresponding to reprocessing the entirety 
of LAT data to make use of the latest available analysis. In addition to being a resource- 
consuming task, re-releases require that users download the newly available data and replace 
all science analysis products. In addition, during the mission we have also made minor im- 
provements in the IRFs based on a better understanding of the properties of the event 
analysis and an improved description of the LAT performance. These incremental IRFs 
can be computed and released without modifying existing data, and many of the analysis 
products remain valid. 

We have released two major iterations of the data analysis since launch: 


• Pass 6 indicates the event analysis scheme designed prior to launch. As such, it was 
based exclusively on our informed estimates of the cosmic-ray environment at the orbit of 
Fermi and a MC-based evaluation of the LAT performance. After the commissioning phase, 
as data started accumulating, we observed phenomena that were not reproduced in the 
MC simulations (see § 2.5 and § 5.2). Without modifying the event analysis in any way, 


we opted to reduce systematic errors by adding these effects to the MC simulations, and 
we re-evaluated the LAT performance (in particular we calculated new IRFs, see § 5.2). 


While this did not allow us to recover any of the lost LAT performance, it ensured that 
real and simulated data were subject to the same effects and the MC-estimated performance 
was therefore adequate for science analysis. We have described the initial Pass 6 release 


(P6_V1) in Atwood et al. (2009), and the corrected IRFs (P6_V3) in Rando & the Fermi LAT 


Collaboration (2009). We will discuss some improvements that were incorporated into the 


later P6_V11 IRFs in § A4 and § H 


• Pass 7 indicates an improved version of the event analysis, for which we updated parts of 
the data reduction process to account for known on-orbit effects by making use of the large 
number of real events the LAT collected in 2 years of operation. The event reconstruction and 
the overall analysis design were not modified, but the event classification was re-optimized 
on simulated data-sets including all known on-orbit effects. Large samples of real events were 
used to assess the efficiency of each step and the systematics involved. Particular attention 
was paid to increasing effective area below ~ 300 MeV where the impact of on-orbit effects 
was large, while maintaining tolerable rates of CR contamination at those energies. Event 
class definitions were optimized based on comparisons of MC events and selected samples of 
real LAT data. See § [3] for a description of Pass 7. 


14 


All data released prior to 2011 August 1 were based on Pass 6. On 2011 August 1 
we released Pass 7 data for the entire mission to date, and since then all data have been 
processed only with Pass 7. 

This paper has two primary purposes. The first is to describe Pass 7 (§ |3j), quanti- 
fying the differences with respect to Pass 6 when necessary. The second is to detail our 
understanding of the LAT, and toward that end we describe how we have used flight data to 
validate the generally excellent fidelity of our simulations of particle interactions in the LAT, 
as well as the resulting IRFs and residual charged particle contamination. In particular, we 
describe the methods and control data samples we have used to study the residual charged 
particle contamination (§|4]), effective area (§ @>, PSF (§§, and energy dispersion (§ 01 of 
the LAT. Furthermore, we quantify the uncertainties in each case, and discuss how these 
uncertainties affect high-level scientific analyses (§[8]). 

For convenience, we have included lists of the acronyms and abbreviations (Appendix [B| 
and notation conventions (Appendix [C]) used in this paper. 


2. LAT INSTRUMENT, ORBITAL ENVIRONMENT, DATA PROCESSING 

AND SIMULATIONS 


In this paper, we focus primarily on those aspects of the LAT instrument, data, and 
analysis algorithms that are most relevant for the understanding and validation of the LAT 
performance. Additional discussion of these subjects was provided in a dedicated paper 


(Atwood et al. 

2009 

). The calibrations of the LAT subsystems are described in a second 

paper ( 

Abdo et al. 

2009f 

)• 


2.1. LAT Instrument 


The LAT consists of three detector subsystems. A tracker/converter (TKR), comprising 
18 layers of paired x-y Silicon Strip Detector (SSD) planes with interleaved tungsten foils, 


which promote pair conversion and measure the directions of incident particles (Atwood 


et al. 2007). A calorimeter (CAL), composed of 8.6 radiation lengths of CsI(Tl) scintillation 


crystals stacked in 8 layers, provides energy measurements as well as some imaging capa- 


bility (Grove & Johnson 2010). An Anticoincidence Detector (ACD), featuring an array 


of plastic scintillator tiles and wavelength-shifting fibers, surrounds the TKR and rejects 
CR backgrounds (Moiseev et al. 2007). In addition to these three subsystems a triggering 
and data acquisition system selects and records the most likely y-ray candidate events for 
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transmission to the ground. Both the CAL and TKR consist of 16 modules (often referred 
to as towers ) arranged in a 4 x 4 grid. Each tower has a footprint of ~ 37 cm x 37 cm and 
is ~ 85 cm high (from the top of the TKR to the bottom of the CAL). A schematic of the 
LAT is shown in Figure [lj and defines the coordinate system used throughout this paper. 
Note that the z-axis corresponds to the LAT boresight, and the incidence ( 6 ) and azimuth 
(0) angles are defined with respect to the z- and x- axes respectively. 



Fig. 1. — Schematic of the LAT, including the layout of the 16 CAL modules and 12 of the 
16 TKR modules (for graphical clarity the ACD is not shown). This figure also defines the 
(6*,0) coordinate system used throughout the paper. 


2.1.1. Silicon Tracker 


The TKR is the section of the LAT where 7 rays ideally convert to e + e _ pairs and their 
trajectories are measured. A full description of the TKR can be found in Atwood et al. 
(2007) and Atwood et al. (2009). Starting from the top (farthest from the CAL), the first 


12 paired layers are arranged to immediately follow converter foils, which are composed of 
~ 3% of a radiation length of tungsten. Minimizing the separation of the converter foils from 
the following SSD planes, and hence the lever-arm between the conversion point and the first 
position measurements, is critical to minimize the effects of multiple scattering. This section 
of the TKR is referred to as the thin or front section. The next 4 layers are similar except 
that the tungsten converters are ~ 6 times thicker; these layers are referred to as the thick 
or back section. The last two layers have no converter; this is dictated by the TKR trigger, 
which requires hits in 3 x-y paired adjacent layers 
to 7 rays that convert in the last two layers. 


see 


3.1.1) and is therefore insensitive 


Thus the TKR effectively divides into two distinct instruments with notable differences 
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CAL 


Fig. 2. — Schematic of a LAT tower (including a TKR and a CAL module). The layout of 
the tungsten conversion planes in the TKR is illustrated. 

in performance, especially with respect to the PSF and background contamination. This 
choice was suggested by the need to balance two basic (and somewhat conflicting) require- 
ments: simultaneously obtaining good angular resolution and a large conversion probability. 
The tungsten foils were designed such that there are approximately the same number of 
7 rays (integrated over the instrument FoV) converted in the thin and thick sections. In 
addition to these considerations, experience on orbit has also revealed that the aggregate of 
the thick layers (~ 0.8 radiation lengths) limits the amount of back-scattered particles from 
the CAL returning into the TKR and ACD in high-energy events (i.e. , the CAL backsplash ) 
and reduces tails of showers in the TKR from events entering the back of the CAL. These 
two effects help to decrease the background contamination in front-converting events. 

After three years of on-orbit experience with the TKR we can now assess the validity 
of our design decisions. The choice of the solid-state TKR technology has resulted in neg- 
ligible down time and extremely stable operation, minimizing the necessity for calibrations. 
Furthermore, the very high signal-to-noise ratio of the TKR analog readout electronics has 
resulted in a single hit efficiency, averaged over the active silicon surface, greater than 99.8%, 
with a typical noise occupancy smaller than 10” 5 for a single readout channel. (We note for 
completeness that the fraction of non-active area presented by the TKR is ~ 11% at normal 
incidence). As discussed below, this has yielded extremely high efficiency for finding tracks 
and has been key to providing the information necessary to reject backgrounds. 

The efficiency and noise occupancy of the TKR over the first three years of operation are 
shown in Figure [3] The variations in the average single strip noise occupancy are dominated 
by one or a few noisy strips, which have been disabled at different times during the mission. 
The baseline of 4 x 10~ 6 is dominated by accidental coincidences between event readouts and 


17 


charged particle tracks (see below, and § 2.1.4) and corresponds to an upper limit of 


noise hits per event in the full LAT on average. Since these noise hits are distributed across 
16 towers and 36 layers per tower, their effect on the event reconstruction is insignificant. 




Fig. 3. — (a) Average TKR hit efficiency and (b) single-strip noise occupancy through the 
first three years of the mission. Each data point is the average value for the full LAT over a 
week of data taking. The shaded background regions mark the first two years of operation, 
corresponding to the data selection used to calibrate the instrument performance. 


The TKR readout is digital , i.e. , the readout is binary, with a single threshold discrim- 
inator for each channel, and no pulse height information is collected at the strip level. The 
individual electronic chains connected to each SSD strip consist of a charge-sensitive pream- 
plifier followed by a simple CR-RC shaper with a peaking time of ~ 1.5 ps. Due to the 
implementation details, the baseline restoration tends to be current-limited, and the signal 
at the output of the shaper is far from the exponential decay characteristic of a linear net- 
work, with the discriminated output being high for ~ 10 p,s for Minimum Ionizing Particles 
(MIPs) at the nominal ~ 1/4 MIP threshold setting. As a consequence, the latched TKR 
strip signals are typically present for that amount of time after the passage of a MIP. If the 
LAT is triggered within this time window, these latent signals will be read out and become 
part of the TKR event data. The rate of occurrence of these ghost signals (which may result 
in additional tracks when the events are reconstructed) depends directly on the charged par- 
ticle background rate. Mitigation against this contamination in the data is discussed below 
[3] and § |5.2 


m 


2.1.4 


Each detector subsystem contributes one or more trigger primitive signals that the LAT 
trigger system uses to adjudicate whether to read out the LAT detectors (see § 3.1). The 


TKR trigger is a coincidence of 3 x-y paired adjacent layers within a single tower (hence a 
total of 6 consecutive SSD detector planes). Due primarily to power constraints, the TKR 
electronic system does not feature a dedicated fast signal for trigger purposes. Rather, the 
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logical OR of the discriminated strip signals from each detector plane is used to initiate a 
non-retriggerable one-shot pulse of fixed length (typically 32 clock ticks or 1.6 (i s), that is 
the basic building block of the three-in- a-row trigger primitive (see § 3.1.1). In addition, 


the length (or time over threshold ) of this layer-OR signal is recorded and included in the 
data stream. Since the time over threshold depends on the magnitude of the ionization in 
the SSD, which in turn depends on on the characteristics of the ionizing particle, it provides 
useful information for the background rejection stage. 


The efficiency of the three-in-a-row trigger is > 99%. This is due in part to the redun- 
dancy of this trigger for the vast majority of events (i.e. , by passing through many layers of 
Si, most events have multiple opportunities to form a three- in-a row). The trigger efficiency, 
measured in flight using MIP tracks crossing multiple towers, is shown in Figure [4] to be very 
stable. 



Fig. 4. — Average efficiency of the TKR three-in-a-row trigger through the first three years 
of the mission. Each data point is the average value for the full LAT over a week of data 
taking. The shaded background region marks the first two years of operation, corresponding 
to the data selection used to calibrate the instrument performance. 


Perhaps the most important figure of merit for the TKR is the resulting PSF for re- 
constructed y-ray directions. At low energy the PSF is dominated by multiple scattering, 
primarily within the tungsten conversion foils (tungsten accounts for ~ 67% of the material 
in the thin section and ~ 92% of the material in the thick section). At high energy the com- 
bination of the strip pitch of 228 /mi, the spacing of the planes and the overall height of the 
TKR result in a limiting precision for the average conversion of ~ 0.1° at normal incidence. 
MC simulations predict that the transition to this measurement precision dominated regime 
should occur between ~ 3 GeV and ~ 20 GeV. We will discuss the PSF in significantly more 
detail in § [6j 
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2.1.2. Electromagnetic Calorimeter 


Details of the CAL can be found in Atwood et al. (2009) and Grove & Johnson (2010). 
Here, we highlight some key aspects for understanding the LAT performance. The CAL 
is a 3D imaging calorimeter. This is achieved by arranging the Csl crystals in each tower 
module in 8 layers, each with 12 crystal logs (with dimensions 326 mm x 26.7 mm x 19.9 mm) 
spanning the width of the layer. The logs of alternating layers are rotated by 90° about the 
LAT boresight, and aligned with the x and y axes of the LAT coordinate system. 


Each log is read out with four photodiodes, two at each end: a large photodiode covering 
low energies (< 1 GeV per crystal), and a small photodiode covering high energies (< 70 GeV 
per crystal). Each photodiode is connected to a charge-sensitive preamplifier whose output 
drives (i) a slow (~ 3.5 /zs peaking time) shaping amplifier for spectroscopy and (ii) a fast 
shaping amplifier (~ 0.5 /zs peaking time) for trigger discrimination. In addition, the output 
of each slow shaper is connected to two separate track-and-hold stages with different gains 
(xl and x8). 


The outputs of the four track-and-hold stages are multiplexed to a single analog-to- 
digital converter. The four gain ranges (two photodiodes x two track-and-hold gains) span an 
enormous dynamic range, from < 2 MeV to 70 GeV deposited per crystal, which is necessary 
to cover the full energy range of the LAT. A zero-suppression discriminator on each crystal 
sparsihes the CAL data by eliminating the signals from all crystals with energies < 2 MeV. 
To minimize CAL data volume, each log end reports only a single range, the lowest-gain 
unsaturated range (the one-range, best-range readout) for most events. For likely heavy 
ions, each log end reports all four ranges (the four-range readout) for calibrating the energy 
scale across the different ranges (see § 3.1.1 for details of how the readout mode is selected). 


The CAL provides two inputs to the global LAT trigger. At each log end the output 
of each fast shaper (for both the large and the small diode) is compared to an adjustable 
threshold by a discriminator to form two separate trigger request signals. In the standard 
science configuration, the discriminator thresholds are set at 100 MeV and 1 GeV energy 
deposition. The 1 GeV threshold is > 90% efficient for incident 7 rays above 20 GeV. 


For each log with deposited energy, two position coordinates are derived simply from 
the geometrical location of the log within the CAL array, while the longitudinal position is 
derived from the ratio of signals at opposite ends of the log: the crystal surfaces were treated 
to provide monotonically decreasing scintillation light collection with increasing distance 
from a photodiode. Thus, the CAL provides a 3D image of the energy deposition for each 
event. 


Since the CAL is only 8.6 radiation lengths thick at normal incidence, for energies 
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greater than a few GeV shower leakage becomes the dominant factor limiting the energy 
resolution, in particular because event-to-event variations in the early shower development 
cause fluctuations in the leakage out the back of the CAL. Indeed, by ~ 100 GeV about 
half of the total energy in showers at normal incidence escapes out the back of the LAT on 
average. The intrinsic 3D imaging capability of the CAL is key to mitigating the degradation 
of the energy resolution at high energy through an event-by-event 3D fit to the shower 


profile. This was demonstrated both in beam tests and in simulations (Baldini et al. 2007 


Ackermann et al. 2010a), achieving better than 10% energy resolution well past 100 GeV. 
The imaging capability also plays a critical role in the rejection of hadronic showers (see 
3.3.7). Furthermore, for events depositing more than ~ 1 GeV in the CAL, imaging in 


the CAL can be exploited to aid in the TKR reconstruction and in determining the event 
direction 


see 


3.2). 


We have monitored the performance of the CAL continuously over its three years of 
operation. We observe a slow (~ 1% per year) decrease in the scintillation light yield of the 
crystals due to radiation damage in low Earth orbit, as we anticipated prior to launch (see 
also § 7.3.2). Although we do not yet correct for this decreased light yield, we have derived 


calibration constants on a channel-by-channel basis for the mission to date. In 2012 January 
we started reprocessing the full data set with these updated calibrations, and in the future 
expect to maintain a quarter-yearly cadence of updates to ensure that the calibrated values 
do not drift by more than 0.5%. 


2.1.3. Anticoincidence Detector 


The third LAT subsystem is the ACD, critically important for the identification of 
LAT-entering charged cosmic rays. Details of its design can be found in Moiseev et al. 


(2007) and Atwood et al. (2009). From the experience of the LAT predecessor, the Energetic 
Gamma Ray Experiment Telescope (EGRET), came the realization that a high degree of 
segmentation was required in order to minimize self-veto due to hard X-ray back-scattering 
(often referred to as backsplash) from showers in the CAL (Esposito et al. 1999). 


The ACD consists of 25 scintillating plastic tiles covering the top of the instrument and 
16 tiles covering each of the four sides (89 in all). The dimensions of the tiles range between 
561 and 2650 cm 2 in geometrical surface and between 10 and 12 mm in thickness. By design, 
the segmentation of the ACD does not match that of the LAT tower modules, to avoid lining 
up gaps between tiles with gaps in the TKR and CAL. The design requirements for the ACD 
specified the capability to reject entering charged particles with an efficiency > 99.97%. To 
meet the efficiency requirement, careful design of light collection using wavelength-shifting 
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fibers, and meticulous care in the fabrication to maintain the maximum light yield were 
needed. The result was an average light yield of ~ 23 photo-electrons (p.e.) for a normally- 
incident MIP for each of the two redundant readouts for each tile. 


The required segmentation inevitably led to less than complete hermeticity, with con- 
struction and launch survival considerations setting lower limits on the sizes of the gaps 
between tiles. Tiles overlap in one dimension, leaving gaps between tile rows in the other. 
The gaps are ~ 2.5 mm and coverage of these gaps is provided by bundles of scintillating 
fibers (called ribbons ), read out at each end by a photomultiplier tube (PMT). The light yield 
for these ribbons, however, is considerably less than for the tiles: it is typically ~ 8 p.e. /MIP 
and varies considerably along the length of a bundle. Therefore, along the gaps the efficiency 
for detecting the passage of charged particles is lower. However, the gaps comprise a small 
fraction of the total area (< 1%) and accommodating them did not require compromising 
the design requirements. In addition to the gaps between ACD tile rows, the corners on the 
sides of the ACD have gaps that are not covered by ribbons and must be accounted for in 
the reconstruction and event classification (see § 3.4). 


As with the TKR and CAL, the ACD provides information used in the hardware trigger 
as well as in the full reconstruction of the event. The output of each PMT is connected to 
(i) a fast shaping amplifier (with ~ 400 ns shaping time) for trigger purposes and (ii) two 
separate slow electronics chains (with ~ 4 /rs shaping time and different gains) to measure 
the signal amplitude. The use of the fast signal in the context of the LAT global trigger 
will be discussed in more detail in § [3j Although the main purpose of the ACD is to signal 
the passage of charged particles, this subsystem also provides pulse height information. The 
two independent slow signals and the accompanying circuitry for automatic range selection 
accommodate a large dynamic range, from ~ 0.1 MIP to hundreds of MIPs. 


As for both the TKR and the CAL, the performance of the ACD has been continuously 
monitored over the past three years. The stability of the MIP signal is shown in Figure [5] and 
in summary shows very little degradation. Note that slight changes in the selection criteria 
and spectra of the MIP calibration event sample cause small (< 0.5%) variations in mean 
deposited energy per event. 


2.1. 4- Trigger and Data Acquisition 

The LAT hardware trigger collects information from the LAT subsystems and, if certain 
conditions are fulfilled, initiates event readout. Because each readout cycle produces a min- 
imum of 26.5 /is of dead time (even more if the readout buffers fill dynamically), the trigger 
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Fig. 5. — Relative position of the MIP peak, calibrated for the tile response and corrected 
for the incidence angle, through the first three years of the mission. Each data point is the 
average value for the 89 ACD tiles over a week of data taking, with the value of the Erst 
point being arbitrarily set to 1. The shaded background region marks the first two years of 
operation, corresponding to the data selection used to calibrate the instrument performance. 
Data from the first several months of the mission are missing because this quantity had not 
yet been included as part of the automated data monitoring and trending. 


was designed to be efficient for 7 rays while keeping the total trigger rate, which is dominated 
by charged CRs, low enough to limit the dead-time fraction to less than about 10% (which 
corresponds to a readout rate of about 3.8 kHz). The triggering criteria are programmable 
to allow additional, prescaled event streams for continuous instrument monitoring and cal- 
ibration during normal operation. We will defer discussion of the actual configuration used 
in standard science operations and of the corresponding performance to the more general 
discussion of event processing in 


3.1 


To limit the data volume to the available telemetry bandwidth, collected data are passed 
to the on-board filter. The on-board filter (see § 3.1.2) consists of a few event selection 
algorithms running in parallel, each independently able to accept a given event for inclusion 
in the data stream to be downlinked. Buffers on the input to the on-board filter can store 
on the order of 100 events awaiting processing. Provided that the on-board filter processes 
at least the average incoming data rate no additional dead time will be accrued because 
the on-board filter is busy. The processors used for the on-board filter must also build and 
compress the events for output, and the time required to make a filter decision varies widely 
between events. Therefore, the event rate that the on-board filter can handle depends on the 
number of events passing the filter. In broad terms, the processors will saturate for output 
rates between 1 kHz and 2.5 kHz, depending on the configuration of the on-board filter. In 
practice, the average output rate is about 350 Hz, and the amount of dead time introduced 
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by the on-board filter is negligible. 

Soon after launch, it became apparent that the LAT was recording events that included 
an unanticipated background: remnants of electronic signals from particles that traversed 
the LAT a few /is before the particle that triggered the event. We refer to these remnants 
as ghosts. An example of such an event is shown in Figure [6j 
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Fig. 6. — Example of a ghost event in the LAT ( y-z orthogonal projection). In addition 
to an 8.5 GeV back-converting y-ray candidate (on the right) there is additional activity in 
all the three LAT subsystems, with the remnants of a charged-particle track crossing the 
ACD, TKR and CAL. The small crosses represent the clusters (i.e. , groups of adjacent hit 
strips) in the TKR, while the variable-size squares indicate the reconstructed location of the 
energy deposition for every hit crystal in the CAL (the side of the square being proportional 
to the magnitude of the energy release). The dashed line indicates the y-ray direction. For 
graphical clarity, only the ACD volumes with a signal above the zero suppression level are 
displayed. 

It is important to re-emphasize a point made in the previous sub-sections: many of 
the signals that are passed to the hardware trigger and the on-board filter are generated by 
dedicated circuits whose shaping times are significantly shorter than for the circuits that read 
out the data from the same sensors for transmission to the ground. Although the signals with 
longer shaping times used for the ground event processing measure the energy deposited in 
the individual channels far more precisely, they also suffer the adverse consequence of being 
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more susceptible to ghost particles. Table [T| gives the characteristic times for both the fast 
signals (i.e., those used in the trigger) and the slow signals (i.e. , those used in the event-level 
analysis) for each of the subsystems. 


Subsystem 

Fast signal 

Slow signal 


(trigger) 

(event data) 

ACD 

0.4 {is 

4 {is 

CAL 

0.5 {is 

3.5 {is 

TKR 

1.5 {is 

10 {is 


Table 1: Characteristic readout time windows for LAT subsystems. The TKR subsystem 
does not provide a dedicated fast signal for the trigger: the peaking time for the shaped 
TKR signal is ~ 1.5 {is (which is the relevant number for the trigger) but the decay time is 
much longer and the average time over threshold for the discriminated output is of the order 
of 10 {is for a MIP at normal incidence. 


The ghost events have been the largest detrimental effect observed in flight data. They 
affect all of the subsystems, significantly complicate the event reconstruction and analysis, 
and can cause serious errors in event classification: 

• they can leave spurious tracks in the TKR that do not point in the same direction as the 
incident 7 ray and might not be associated with the correct ACD and CAL information due 
to the different time constants of the subsystems; 

• they can leave sizable signals in the CAL that can confuse the reconstruction of the 
electromagnetic shower, degrading the measurement of the shape and direction of the shower 
itself, and that can cause the energy reconstruction to mis-estimate the energy of the incident 

7 r ay; 

• they can leave significant signals in the ACD that can be much larger than ordinary 
backsplash from the CAL and that can cause an otherwise reconstructable 7 ray to be 
classified as a CR. 

We will discuss the effects on the LAT effective area in more detail in § [5j In brief: 
above ~ 1 GeV, where the average fractional ghost signal in the CAL is small, the loss of 
effective area is of the order of 10% or less. This loss is largely due to the fact that valid 7 
rays can be rejected in event reconstruction and classification if ghost energy depositions in 
the CAL cause a disagreement between the apparent shower directions in CAL and TKR. 
At lower 7 -ray energies ghost signals can represent the dominant contribution to the energy 
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deposition in the instrument, and the corresponding loss of the effective area can exceed 20%. 
We also emphasize that, while these figures are averaged over orbital conditions, the fraction 
of events that suffer from ghost signals, as well as the quantity of ghost signals present in 
the affected events, varies by a factor of 2-3 at different points in the orbit due to variation 
in the geomagnetic cutoff and resulting CR rates. Recovering the effective area loss due to 
ghost signals is one of the original and main motivations of ongoing improvements to the 
reconstruction (Rochester et al. 2010). 


Finally, during extremely bright Solar Flares (SFs) the pile-up of several > 10 keV X- 
rays within a time interval comparable with the ACD signal shaping time can also cause 7 
rays to be classified as CRs (see Appendix [A]) . 


2.2. Orbital Environment and Event Rates 

Fermi is in a 565- km altitude orbit with an inclination of 25.6°. The orbit has a period 
of ~ 96 minutes, and its pole precesses about the celestial pole with a period of ~ 53.4 days. 
At this inclination Fermi spends about 15% of the time inside the South Atlantic Anomaly 
(SAA). Science data taking is suspended while Fermi is within the SAA because of the high 


flux of trapped particles (for details, see Abdo et al. 2009f) 


When Fermi is in standard sky-survey mode the spacecraft rocks north and south about 
the orbital plane on alternate orbits. Specifically, the LAT boresight is offset from the 
zenith toward either the north or south orbital poles by a characteristic rocking angle. On 
2009 September 3 this rocking angle was increased from 35° to 50° in order to lower the 
temperature of the spacecraft batteries and thus extend their lifetime. As a result of this 
change, the amount of the Earth limb that is subtended by the FoV of the LAT during 
survey-mode observations increased substantially. The most noticeable consequence is a 
much larger contribution to the LAT data volume from atmospheric 7 rays. This will be 
discussed more in § [4], and further details about atmospheric 7 rays can be found in Abdo 


et al. (2009a). 


The flux of charged particles passing through the LAT is usually several thousand times 
larger than the 7-ray flux. Accordingly, several stages of event selection are needed to purify 
the 7-ray content (see § 3.1, § 2.4 and § [3]) and some care must be taken to account for 


contamination of the 7-ray sample by charged particles that are incorrectly classified as 7 
rays (see § [4]) . Here we distinguish 4 stages of the event classification process: 


1. hardware trigger request : the LAT detects some traces of particle interaction and starts 
the triggering process (§ 3.1); 
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2. hardware trigger accept : the event in question generates an acceptable trigger pattern 
and is read out and passed to the on-board filter (§ 0; 

3. on-hoard filter: the event passes the on-board 7-ray filter criteria 


3.1); 


4. standard 7- ray selection: the event passes more stringent selection criteria (§ 3.4), such 
as P7S0URCE, which is the selection currently being recommended for analysis of individual 
point sources, or P6_DIFFUSE, the selection recommended for analysis of point sources in the 
Pass 6 iteration of the event selections. 


Figure [7] shows an example of both the orbital variations of the charged particle flux, 
and how the initially overwhelming CR contamination is made tractable while maintaining 
high efficiency for 7 rays by several stages of data reduction and analysis. Figure [8] shows how 
the charged particle rate varies with spacecraft position. In particular, data taken during 
the parts of the orbit with the highest background rates are more difficult to analyze for 
two reasons: (i) there are simply more non-7-ray events that must be rejected, and (ii) the 
increased ghost signals complicate the analysis of the 7-ray events (see 


2.5.1). 


A model of the particle backgrounds for the region outside the SAA was compiled 


before launch as documented in Atwood et al. (2009): above the geomagnetic cutoff rigidity 


(~ 10 GV) the most abundant species is primary galactic CR protons, at lower energies 
particle fluxes are dominated by trapped protons and electrons, in addition to the 7 rays 
due to interaction of CRs with the atmosphere of the Earth. Since launch, the particle model 
has been updated to include new measurements, (see, e.g., 
electron population). 


Ackermann et al. 2010a, for the 


2.3. Observing Strategy and Paths of Sources Across the LAT Field-of-View 

Fermi has spent over 95% of the mission in sky survey mode, with most of the remainder 
split between pointed observations and calibrations. Furthermore, in sky survey mode, the 
azimuthal orientation of the LAT is constrained by the need to keep the spacecraft solar 
panels pointed toward the Sun and the radiators away from the Sun. Specifically, in sky 
survey mode Fermi completes one full rotation in 0 every orbit. 

Therefore, during an orbital precession period the LAT boresight will cross a range 
of declinations ±25.6° relative to the rocking angle, while the LAT x-axis will track 
the direction toward the Sun. As a result, during the sky survey mode, each point 
in the sky traces a complicated path across the LAT FoV, the details of which depend on 
the declination and ecliptic latitude. Figure [9] shows examples of these paths for two of 
the brightest LAT sources. Furthermore, in sky survey mode operations the path 
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Fig. 7. — Rates at several stages of the data acquisition and reduction process on a typical 
day (2011 August 17). Starting from the highest, the curves shown are for the rates: (i) at the 
input of the hardware trigger process (trigger request), (ii) at output of the hardware trigger 
(trigger accept), (iii) at the output of the on-board filter, (iv) after the loose P7TRANSIENT 
y-ray selection, (v) after the tighter P7S0URCE y-ray selection, and (vi) the P7S0URCE y-ray 
selection with an additional cut on the zenith angle (0 Z < 100°). See § [3] for more details 
about the event selection stages. 

that the LAT boresight traces across the sky during any two orbit period is only 
slightly different than during the two previous or subsequent orbits. 


2.4. Ground-Based Data Processing 

Reconstructing the signals in the individual detector channels into a coherent picture 
of a particle interaction with the LAT for each of the several hundred events collected every 
second is a formidable task. We will defer detailed discussion of the event reconstruction 
and classification to § |3j here we describe just the steps to give a sense of the constraints. 


1. Digitization-, we decompress the data and convert the information about signals in indi- 
vidual channels from the schema used in the electronics readout to more physically motivated 
schema — such as grouping signals in the ACD by tile, rather than by readout module. 


2. Reconstruction : we apply pattern recognition and fitting algorithms commonly used in 
high-energy particle physics experiments to reconstruct the event in terms of individual TKR 
tracks and energy clusters in the CAL and to associate those objects with signals in the ACD 
(see § |3.2h . 


3. Event analysis: we evaluate quantities that can be used as figures of merit for the event 



Geographic longitude [°] Geographic longitude [°] Geographic longitude [°] 


Fig. 8. — Orbital background environment as a function of geographic location: (a) vertical 
geomagnetic cutoff rigidity; (b) mean rate of events at the input of the hardware trigger 
process (trigger requests, see § 3.1.1), which decreases roughly as the square root of the 
geomagnetic cutoff rigidity; (c) mean rate of events in P7S0URCE event class (see § 3.4.3) 
with the additional cut on the zenith angle (6 Z < 100°). The anti-correlation between the 
P7S0URCE event rate and the trigger request rate is the result of the loss of efficiency due 
to ghost signals (see § 2.1.4 and § 5.2.2). The black lines and points surrounding the white 


area represent the sides and vertices of the SAA polygon as defined for LAT operations. The 
LAT does not acquire science data while inside the SAA polygon. 
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Fig. 9. — Live time maps in instrument coordinates for the Vela pulsar (a) and the Crab 
(b). These plots are made in a zenith equal area projection with the LAT boresight at the 
center of the image, and the color scale shows the total live time each source was at 9 Z < 85°. 
Recall that 0 = 0 corresponds to the +x axis of the LAT (Fig. [Tj) . These plots cover the 
same two-year time range that defines the standard calibration samples (see § 3.6). 


from the collections of tracks, clusters and associated ACD information. Once we have ex- 
tracted this information we apply multi- variate analysis techniques to extract measurements 
of the energy and direction of the event and to construct estimators that the event is in fact 
a y-ray interaction (see 


3.3). 


4. Event selection: we apply the selection criteria for the various y-ray event classes (§ 3.4). 


In addition to these procedures, the processing pipeline is responsible for verifying the data 
integrity at each step and for producing and making available all the ancillary data products 
related to calibration and performance monitoring of the LAT. 

On the whole, the LAT data processing pipeline utilizes approximately 130 CPU years 
and 300 TB of storage each year. 

The ISOC can reprocess data with updated algorithms from any stage in the process. 
However, given the quantities of data involved, reprocessing presents challenges both to the 
available pool of computing power and storage space. Table [2] gives a rough idea of the 
resources used by each stage in the process. Given the available resources, reprocessing from 
the Reconstruction step is practical only once a year or less often, and reprocessing from the 
Digitization step is feasible only every few years. 


30 


Stage 

CPU Time 
[years/year] 

Disk Usage 
[TB/year] 

Digitization 

1.2 

26 

Reconstruction 

85 

167 

Event analysis 

5 

10 

Event selection 

1.3 

60 

Data monitoring 

25 

12 

Calibration 

5 

72 


Tabic 2: Various stages of the data processing pipeline. The numbers in the CPU time and 
disk usage columns refer to a typical year of data taking. 


2.5. Simulated Data 


In order to develop the filtering, reconstruction and event selection algorithms described 
in § [3] we used a detailed simulation of particle interactions with Fermi written in the 
Geant4 (Agostinelli et al. 2003) framework. This simulation includes particle interactions 


with a detailed material model of Fermi , as well as simulations of the uncalibrated signals 
produced in the various sensors in the three subsystems. 

The fact that the simulation is detailed enough to produce uncalibrated signals for each 
channel allows the simulations to be truly end-to-end: 


1. we maintain fidelity to the analysis of the flight data by processing simulated data with 
the same reconstruction and analysis algorithms as the flight data; 


2. on a related but slightly different note, we simulate the data as seen by the trigger and 
the on-board software, and process the data with a simulation of the hardware trigger and 
exactly the same on-board filter algorithms as used on the LAT; 


3. we can merge the signals from two events into a single event. As described in § 2.5.1 
rely on this feature to add an unbiased sample of ghost signals to simulated events. 


we 


2.5.1. Ghosts and Overlaid Events 


The presence of ghost signals in the LAT data reduced the efficiency of the event se- 
lections tuned on simulated data lacking this effect. So we developed a strategy to account 
for the ghosts. One of the triggers implemented in the LAT data-acquisition system, the 
PERIODIC trigger (see § 3.1.1), provides us with a sample of ghost events: twice per second 
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the LAT subsystems are read out independent of any other trigger condition. For the CAL 
and the ACD, all the channels, even those below threshold, are recorded (this is not possible 
for the TKR). Since these triggers occur asynchronously with the particle triggers, the ghost 
signals collected are a representative sample of the background present in the LAT. 


We merge the signals channel-by-channel from a randomly chosen periodic trigger event 
into each simulated event. In more detail, we analyze the signals in each periodic trigger, 
converting the instrument signals to physical units; analog-to-digital converter channels in 
the ACD and CAL, and time over threshold signals in the TKR, are converted into deposited 
energy. Since the intensity and spectrum of cosmic rays seen by the LAT depend on the 
position of the Fermi satellite in its orbit, the ghost events are sorted by the value of 
Mcllwain L (Walt;!2005) at the point in the orbit where the event was recorded. Then, during 
the simulation process, for each simulated event we randomly select a periodic trigger event 
with similar Mcllwain L, and add the energy deposits in this “overlay event” to those of the 
original simulated event, after which the combined event is digitized and reconstructed in 
the usual way. 


We have used the improved simulations made with this overlay technique to more ac- 
curately evaluate the effective area of the LAT, and to better optimize the event selection 
criteria when developing Pass 7. Furthermore, we are studying ways to identify and tag 


the ghost signals and remove them from consideration in the event analysis (Rochester et al. 

2010 ). 


2.5.2. Specific Source 'y-ray Signal Simulation 


We have developed interfaces between the software libraries that implement the flux 
generation and coordinate transformation used by the ScienceTools and our detailed Geant4- 
based detector simulation. This allows us to re-use the same source models that we use with 
the ScienceTools within our detailed Geant4-based detector simulation, insuring consistent 
treatment of Fermi pointing history and the source modeling. One application for the present 
work was simulating atmospheric 7 rays above 10 GeV for comparison with our Earth limb 
calibration sample (see § 3.6.3). 


2.5.3. Generic y-ray Signal Simulation 

In addition to simulating individual sources, we also simulate uniform 7 -ray fields that 
we can use to explore the instrument response across the entire FoV and energy range of the 
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LAT. These will henceforth be referred to as allGamma simulations, and they have these 
features: 

1. the 7 rays are generated with an A 1-1 number spectrum (where E is the 7-ray energy), 

1. e., the same number of 7 rays are generated in each logarithmic energy bin; 

2. the individual 7 rays are randomly generated on a sphere with 6 m 2 cross sectional area 
(i.e., large enough to contain the entire LAT and most of the spacecraft) centered at the 
origin of the LAT reference frame, i.e., the center of the TKR/CAL interface plane; 

3. the directions of the 7-rays are sampled randomly across 27T of downward-going (in the 
LAT reference frame) directions, so as to represent a semi-isotropic incident flux. 

For defining the IRFs we simulate 200 million 7 rays over the range log 10 (E/1 MeV) G 
[1.25,5.75] and a further 200 million 7 rays over the range log l0 (E/l MeV) € [1.25,2.75]. 
The net result is to produce a data set that adequately samples the LAT energy band and 
FoV and that far exceeds the statistics of any single point source. 


2.5.4 ■ Simulation of Particle Backgrounds 

In order to develop our event classification algorithms, and to quantify the amount of 
residual CR background that remains in each 7-ray event class, we require accurate models 
of both the CR background intensities and the interactions of those particles with the LAT. 
Accordingly, we use the same Geant4-based detector simulation described in the previous 
section to simulate fluxes of CR backgrounds. 

There are three components to the CR-induced background in the LAT energy range. 

• Primary CRs : Protons, electrons, and heavier nuclei form the bulk of charged CRs. Due 
to the shielding by the magnetic field of the Earth, only particles with rigidities above a few 
GV (spacecraft position dependent) can reach the LAT orbit. Therefore, the background 
from primary CRs is relevant above a few GeV. The LAT event classification provides a 
very effective rejection of charged particles, up to an overall suppression factor of 10 6 for 
CR protons (with a prominent contribution from the ACD). However, due to the intense 
flux of primary protons, a significant number of misclassified CRs enter even the P7S0URCE, 
P7CLEAN, and P7ULTRACLEAN 7-ray classes. Minimum-ionizing primary protons above a few 
GeV produce a background around 100-300 MeV in the LAT. However, protons that in- 
elastically scatter in the passive material surrounding the ACD (e.g., the micro-meteoroid 
shield) without entering the ACD can produce secondary 7 rays at lower energies which are 
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then indistinguishable from cosmic 7 rays, and which we refer to as irreducible background 
(see also § 4.4). 


• Charged secondaries from CR interactions'. CRs entering the atmosphere produce a num- 
ber of secondaries in interactions with the atmosphere itself. Charged particles produced in 
these interactions can spiral back out, constrained by the magnetic held of the Earth, and 
eventually re-enter the atmosphere. These particles are predominantly protons, electrons, 
and positrons and are an important background below the geomagnetic cutoff rigidity. Due 
to the high efficiency of the ACD in rejecting charged particles, the dominant contribution 
to this background in the P7CLEAN and P7ULTRACLEAN event classes is positron annihilations 
in the passive material around the ACD. The resulting 7 rays are again indistinguishable 
from cosmic 7 rays. 


• Neutral secondaries from CR interactions'. 7 rays and neutrons are also produced in local 
CR interactions. Unaffected by the magnetic held, they approach the LAT from the direction 
of the Earth with intensities peaking at the Earth limb. Neutrons are easily suppressed by 
the event classification scheme and do not contribute significantly to the background. The 
7-ray background is suppressed by pointing away from the Earth and excluding events that 
originate from near or below the Earth limb. Due to the long tails of the LAT PSF (see §[6]), 
however, a small fraction of the events originating from the Earth limb are reconstructed 


with directions outside the exclusion regions. Since the Earth limb is extremely bright (Abdo 


et al. 2009b), even this small fraction is an important residual background for celestial 7 rays. 
Since the PSF tails broaden with decreasing energy the main contribution is at energies below 
a few hundred MeV (see § 4.5). Similarly, a small fraction of 7 rays entering the LAT with 
incidence angles greater than 90° are mischaracterized as downward-going (~ 5% to ~ 0.01%, 
depending on the energy, incidence angle, and the event sample). Because the Earth limb is 
extremely bright and some part of it is almost always behind the LAT, 7 rays from the limb 
are the dominant component of “back-entering” background contamination. Furthermore, 
since we mis-estimate the directions of these back-entering 7 rays by > 90°, they are often 
reconstructed outside the Earth exclusion region. 


Figure 10 shows the average CR-induced particle intensities at the orbit of the LAT 
in the model that we use. For comparison the intensity of the extragalactic diffuse 7-ray 
emission measured by the LAT (Abdo et al. 2010f) is overlaid to demonstrate the many 
orders of background suppression necessary to distinguish it from particle background. The 
model was developed prior to launch based on data from satellites in similar orbits and 
balloon experiments (Mizuno et al. 2004). 


As the particle rates are strongly dependent on location in geomagnetic coordinates, the 
details of the orbit model are also important. For tuning the event analysis, or for estimating 
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Fig. 10. — Model of the LAT orbital position and particle direction-averaged CR-induced 
particle intensities (Mizuno et al. 2004) sampled from a 64 s live time background simulation 
run. The intensity of the extragalactic diffuse background emission measured by the LAT 
(Abdo et al. 2010f) is shown for comparison. Note that the event energy is reconstructed 
under the hypothesis of a downward-going 7 ray and in general does not represent the actual 
energy for hadrons. 


the background rates for typical integration times of months or years, the simulated time 
interval must be at least equal to the precession period of the Fermi orbit (53.4 days). 
Simulating these high particle rates for such a long time interval is quite impractical, in 
terms of both CPU capacity and disk storage requirements. For studies of background 
rejection we usually simulate an entire precession period to ensure a proper sampling of the 
geomagnetic history, but to limit the particle counts we generate events for only a few seconds 
of simulated time every several minutes, e.g., a typical configuration requires event generation 
for 4 seconds every 4 minutes of time in orbit. This partial sampling is a compromise between 
the limited CPU and disk usage, and the requirement of having good statistics. Considering 
the LAT background rejection power, in order to have sizable statistics after even the first 
stages of the event analysis are performed, we must start with a simulated background data 
set of over 10 9 CRs. 
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3. EVENT TRIGGERING, FILTERING, ANALYSIS AND 

CLASSIFICATION 


In this section we describe the analysis steps that determine which events make it into 
our final 7 -ray data sample, starting with the triggering and filtering performed by the on- 
board data acquisition system 


3.1 ), moving on to the reconstruction of particle interactions 


in the event 


3.3) and finally the definition 


3.2), the analysis of the event as a whole 
of the 7 -ray classes (§ |3T). The overall logical structure of this process is illustrated in 
Figure 11 


Trigger and filter 
§ 3.1 


Event reconstruction 
§ 3.2 


Event-level analysis 
§ 3.3 



Definition of standard photon classes 
S 3.4 



Fig. 11. — Logical structure of the analysis steps determining which events are selected for 
a given class. The references to section numbers are meant to help the reader navigate the 
content of § [3j 


The event analysis requires knowledge of the LAT, the physics of particle interactions 
within its volumes, and of the particle backgrounds in the Fermi orbit. As described in 


2.5, we use large MC samples of 7 rays and of CRs to devise the best procedures to extract 


estimates of energies and incident directions, and to classify events as either 7 rays or charged 
particle backgrounds. 


Finally, in § 3.5 we describe the publicly available LAT event samples, while in § 3.6 
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we describe the calibration sources, event samples and methods we use to validate and 
characterize the performance of the LAT using flight data. 


3.1. Trigger and On-Board Filter 

In this section we review the event triggering, the readout of the LAT and the filtering 
performed on-board in order to reduce the data volume downlinked to ground. 


3.1.1. Triggering the LAT Readout 

Each subsystem provides one or more trigger primitives (or trigger requests ) as detailed 
in the following list. 


• TKR (also known as “three- in-a-row” ) : issued when three consecutive x-y silicon layer 
pairs — corresponding to six consecutive silicon planes — have a signal above threshold (nom- 
inally 0.25 MIPs). This signals the potential presence of a track in a tower. Since many 
tracks cross between towers and/or have more than 3 x-y layers within a tower, the TKR 
trigger request is very efficient. 

• CAL_L0: issued when the signal in any of the CAL crystal ends crosses the low-energy 
trigger threshold (nominally 100 MeV). 

• CAL_HI: issued when the signal in any of the CAL crystal ends crosses the high-energy 
trigger threshold (nominally 1 GeV). 

• VETO: issued when the signal in any of the ACD tiles is above the veto threshold (nom- 
inally 0.45 MIP). It signals a charged particle crossing the tile. The trigger system has 
a programmable list of ACD tiles associated with each TKR tower; if ACD shadowing is 
enabled (which is the case for the nominal configuration) an additional ROI primitive is as- 
sembled when a TKR primitive in a tower happens in coincidence with a VETO primitive in 
the ACD tiles associated to the same tower. 


• CNO: issued when the signal in any of the ACD tiles is above the CNO threshold (nominally 
25 MIPs). It indicates the passage of a heavily ionizing nucleus (CNO stands for “Carbon, 
Nitrogen, Oxygem ) and it is primarily used for the calibration of the CAL (see § 7.3.2). 


The LAT has the ability to generate three other trigger primitives. We will refer to 
these as special primitives. Two of these are unused in flight and will not be discussed here. 
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The third, the PERIODIC trigger, runs at a nominal frequency of 2 Hz during all science data 
taking and is used for diagnostic and calibration purposes. 

We emphasize that although the trigger primitives provided by the TKR and the CAL 
are tower-based, a trigger initiates a readout of the whole LAT. 

Trigger primitives are collected by the Central Trigger Unit. All 256 possible combina- 
tions of the eight trigger primitives are mapped into so-called trigger engines. In brief, some 
trigger requests are allowed to open a trigger window of fixed duration (nominally 700 ns) 
and the set of primitives collected in this time interval is compared to a table of allowed 
trigger conditions. In case a trigger condition is satisfied, a global trigger (or trigger accept) 
is issued and event acquisition commences. The LAT readout mode (i.e., enabling the CAL 
and ACD zero suppression and selecting the CAL one-range or four-range readout) can be 
individually set for each trigger engine. In addition, trigger engines are scalable: for each 
trigger condition a prescale is specified, corresponding to the number of valid trigger requests 
necessary to issue a single global trigger (obviously no prescale is applied to engines intended 
for y-ray collection). 


Engine 

PERIODIC 

CAL _HI 

CAL_L0 

TKR 

R0I 

CN0 

Prescale 

Average rate [Hz] 

3 

1 

X 

X 

X 

X 

X 

0 

2 

4 

0 

X 

1 

1 

1 

1 

0 

200 

5 

0 

X 

X 

X 

X 

1 

250 

5 

6 

0 

1 

X 

X 

X 

0 

0 

100 

7 

0 

0 

X 

1 

0 

0 

0 

1500 

8 

0 

0 

1 

0 

0 

0 

0 

400 a 

9 

0 

0 

1 

1 

1 

0 

0 

700 

10 

0 

0 

0 

1 

1 

0 

50 

100 


Table 3: Definition of the standard trigger engines: primitives used (1: required, 0: excluded, 
x: either), prescale factors and typical rates. In this short-hand representation, engines are 
defined with highest precedence first: each combination of trigger primitives is mapped to 
the engine corresponding to the first condition it matches. Trigger engines 0, 1 and 2 are 
used by the LAT software for bookkeeping and to catch conditions that should not happen 
on orbit. 


In the nominal configuration for science data taking the CAL_L0 condition is inhibited from opening a trigger 
window, and therefore engine 8 is effectively disabled. 
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The standard trigger engine definitions are listed in Table [3j Trigger engine 7 is partic- 
ularly important for 7-ray events: it requires no special primitives, a three-in-a-row signal in 
the TKR (TKR), absence of ROI shadowing and CNO from the ACD, and no CAL_HI. Engine 6 
handles CAL_HI events and engine 9 handles events with enough energy deposition in the 
CAL to potentially cause ROI shadowing. 

Trigger engine 4 is typical for calibration data (e.g., heavy ion collection): it requires no 
special primitives, three-in-a-row in the TKR (TKR), a highly ionizing passage in the ACD 
(CNO) in close correspondence (ROI) and at least a moderate energy deposition in the CAL 
(CAL_L0). Engines 5 and 10 process very few 7 rays and are used primarily for calibration and 
monitoring, so they are prescaled to limit the readout rate and/or to match the downlink 
bandwidth allocation. 


3.1.2. Event Filtering 

As mentioned before, since limited telemetry bandwidth is available for data downlink, 
some event filtering is required. The on-board filter allows the coexistence of different filtering 
algorithms and in fact, in the nominal science data taking configuration, all the events are 
presented to three different filters: 

• the gamma filter, designed to accept 7 rays (with an output average rate of approximately 
350 Hz); 

• the HIP filter, designed to select heavy ion event candidates, primarily for CAL calibration 
(with an output average rate of approximately 10 Hz); 

• the DIAGNOSTIC filter, designed to enrich the downlinked data sample in events useful to 
monitor sensor performance and selection biases: specifically ~ 2 Hz of periodic triggers and 
an unbiased sample of all trigger types prescaled by a factor 250 (with an output average 
rate of approximately 20 Hz). 

The gamma filter consists of a hierarchical sequence of veto tests, with the least CPU- 
intensive tests performed first. If an event fails a test, it is marked for rejection and not passed 
on for further processing (however, a small subset of events that are marked for rejection 
are downlinked through the DIAGNOSTIC filter). Some of the tests utilize rudimentary, two- 
dimensional tracks found by the on-board processing. The gamma filter has several steps 
(listed in order of processing): 

1. reject events that have patterns of ACD tile hits that are consistent with CRs and do 
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not have the CAL_L0 trigger primitive asserted, making it unlikely that the ACD hits were 
caused by backsplash; 

2. accept all events for which the total energy deposited in the CAL is greater than a 
programmable threshold, currently set to 20 GeV; 

3. reject events that have ACD hit tile patterns that are spatially associated with the TKR 
towers that caused the trigger, provided that the energy deposited in the CAL is less that a 
programmable threshold (currently set to 350 MeV); 

4. reject events for which a significant energy deposition in the CAL (typically > 100 MeV) 
is present but the pattern of hits in the TKR is unlikely to produce any track; 

5. reject events for which rudimentary tracks match with individual ACD tiles that were hit, 
provided the energy deposited in the CAL is less than some programmable amount (typically 
5 GeV); 

6. reject events that do not have at least one rudimentary track. 

Although it may seem strange to apply the requirement that there be any tracks after cutting 
on matches between tracks and the ACD, recall that the on-board filter software is highly 
optimized for speed, and terminates processing of each event as soon as it is possible to reach 
a decision. Thus the testing for track matches is performed during the track-finding stage, 
at the time the tracks are constructed. 


3.2. Reconstruction Algorithms 


Event reconstruction translates the raw event information from the LAT subsystems 
into a high-level event description under the assumption of a 7 ray impinging on the LAT 


volumes within 90° of the boresight (see Fig. 12 for an illustrative event display). 


Here we will briefly summarize the event reconstruction algorithms underlying both 
Pass 6 and Pass 7 flavors of the event analysis; a more detailed description is given in 


Atwood et al. (2009). We want to stress that Pass 6 and Pass 7 use exactly the same 


reconstruction algorithms with the exception of the energy un-biasing (described in § 3.2.1), 
which was only applied to Pass 7 data. 


3.2.1. Calorimeter Reconstruction 

The starting point for the energy evaluation is the measured energy depositions in the 
crystals. The centroid of the energy deposition is determined and the principal axes of the 
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Fig. 12. — Event display of a simulated 27 GeV 7 ray (a) and zoom over the CAL (b) 
and TKR (c) portions of the event. The small crosses represent the clusters in the TKR, 
while the variable-size squares indicate the reconstructed location and magnitude of the 
energy deposition for every hit crystal in the CAL. The dotted line represents the true 7 -ray 
direction, the solid line is the CAL axis (§ 3.2.1) and the dashed lines are the reconstructed 


TKR tracks (§ 3.2.1 ). The backsplash from the CAL generates tens of hits in the TKR, with 


two spurious tracks reconstructed in addition to the two associated with the 7 ray (note 
that they extrapolate away from the CAL centroid and do not match the CAL direction). 
It also generates a few hits in the ACD, which, however, are away from the vertex direction 
extrapolation and therefore do not compromise our ability to correctly classify the event as 
a 7 ray. 
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shower are evaluated by means of a principal moment analysis. In the Pass 6 and Pass 7 
event reconstruction procedure, the energy deposition is treated as a single quantity, with 
no attempt to identify contamination from ghost signals. Work to develop an algorithm to 
separate the CAL energy deposition into multiple clusters and to disentangle ghost signals 
is ongoing (Rochester et al. 2010). The amount of energy deposited in the TKR is evaluated 
by treating the tungsten-silicon detector as a sampling calorimeter; this contribution is an 
important correction at low energies. 


We apply three algorithms to estimate the actual energy of an event: a Parametric 
Correction (PC), a fit of the Shower Profile (SP) and a Maximum Likelihood (LH) approach. 
The energy assigned to any given event is the energy from one or the other of these algorithms. 
The three methods were designed to provide the best performance in different parts of the 
energy and incidence angle phase space (in fact, the LH algorithm was only tuned up to 
300 GeV, while the SP algorithm does not work well for events below ~ 1 GeV). Accordingly, 
they provide different energy resolutions and their distributions have slightly different biases 
(i.e., the most probable values are slightly above or below the true energy) for different 
energies and incidence angles; more details can be found in Atwood et al. (2009). 


In fact, the only significant change in the Pass 7 event reconstruction relative to Pass 6 
is to apply separate corrections for the biases of each energy estimation algorithms. We used 
MC simulations to characterize the deviations of the most probable value of the energy 
dispersion from the true energy across the entire LAT phase space for the three methods 
separately. Such deviations were found to be typically of the order of few percent (with a 
maximum value of ~ 10%) and always significantly smaller than the energy resolution — with 
LH displaying a negative bias and PC and SP displaying a positive bias in most of the phase 
space. 


We generated correction maps (as functions of y-ray energy and zenith angle) and in 
Pass 7 the residual average bias for all the inputs of the final energy assignment (discussed 


m 


3.3.2) is less than 1% in the entire LAT phase space. 


3.2.2. Tracker Reconstruction 

For the TKR we merge and assemble clusters of adjacent hit strips into track candi- 
dates by combinatorial analysis. We have developed two methods: for CAL-Seeded Pattern 
Recognition (CSPR) the trajectory of the original 7 ray is assumed to point at the centroid 
of the energy released in the CAL; the Blind Search Pattern Recognition (BSPR) can be 
used when there is little or no energy deposit in the CAL sensitive volumes. 
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Both the CSPR and BSPR algorithms start by considering nearby pairs of TKR clusters 
in adjacent layers as candidate tracks (the CSPR algorithm limits the candidate tracks to 
those pointing toward the CAL energy centroid). Both algorithms then proceed by using 
a Kalman filtering technique (see Kalman 1960 Fruhwirth 1987, for details about Kalman 
filtering and its use in particle tracking) which tests the hypotheses that additional TKR 
clusters were generated by the same incident particle and should be associated to each 
of the candidate tracks, and adds any such clusters to the appropriate candidate tracks. 
Furthermore, as each cluster is added to candidate tracks the Kalman filter updates the 
estimated direction and associated covariance matrix of those tracks. Both the CSPR and 
BSPR algorithms are weighted to consider best candidate track to be the one that is both 
pointing toward the CAL centroid and is the longest and straightest. In the CSPR, the 
main axis of the CAL energy deposition is also considered, candidate tracks for which the 
TKR and CAL estimated directions differ significantly are disfavored starting at ~ 1 GeV, 
and increasingly so at higher energies. At the completion of the CSPR algorithm the best 
candidate track is selected and confirmed as a track, and the clusters in it are flagged as 
used. We iterate the CSPR algorithm until no further tracks can be assembled from the 
unused TKR clusters, then proceed with the BSPR. 


If more than one track is found in a given event, we apply a vertexing algorithm that 
attempts to compute the most likely common origination point of the two highest quality 
(i.e. , longest and straightest) tracks, and, more importantly, to use that point as a constraint 
in combining the momenta of the two tracks to obtain a better estimate of the direction of 
the incoming 7 ray. 


3.2.3. ACD Reconstruction 

The ACD phase of the event reconstruction starts by estimating the energy deposited 
in each of the tiles and ribbons. Subsequently, these energy depositions are associated to 
each of the tracks found in the TKR. More specifically, each track is projected to the ACD, 
and we calculate whether that track intersects each ACD tile or ribbon with non-zero energy 
deposition. Furthermore, if the track projection does not actually cross an ACD element 
with non-zero energy deposition, we calculate the distance of closest approach between the 
track projection and the nearest such ACD element. Finally, we use the distance calculation 
and the amount of energy deposited to sort all the possible ACD-TKR associations by how 
likely they are to represent a charged particle passing through the ACD and into the TKR. 
This ranking is used when considering whether the event should be rejected in later analysis 
stages. 


43 


3.3. Event Analysis 


The first step of the event-level analysis procedure is to extract simpler representations 
of the event characteristics from the complex structures assembled during the reconstruction 
phase. These high-level quantities are then used for the final assignment of the event energy 
and direction (among the outputs of the reconstruction algorithms described in § 3.2) and for 


the background rejection. The final product is an event-by-event array of simple quantities 
relevant for scientific analysis: energy, direction and estimates of the probability a given 
event is a 7 ray (P a n). 


3.3.1. Extraction of “Figure- of- Merit” Quantities 

For each event we reduce the output of the TKR, CAL and ACD reconstruction to a 
set of a few hundred figure-of- merit quantities whose analyzing power has been studied and 
optimized with MC simulations. 

It is important to note that the best track (and to a lesser extent the second best track) 
plays a particularly important role in later analyses. Specifically, although we calculate 
figure-of- merit variables — such as the agreement between the track direction and the axis of 
the CAL shower, or the distance between the track extrapolation and the nearest ACD energy 
deposition — for all the tracks in the event, many of the multi- variate analyses described in the 
rest of this section consider only those figure-of-merit variables associated with the two best 
tracks in the event. Furthermore, the figure-of-merit variables associated with the best track 
tend to carry significantly more weight in the multi-variate analysis — which is of primary 
importance, as most ghost tracks come from protons and heavy ions and tend to be longer 
and straighter than tracks from e + and e~ from a 7 -ray conversion. 

Early iterations of the event analysis split the events by energy and then applied different 
selections in the various energy bands. We found that this approach led to large spectral 
features at the bin boundaries. Therefore, we chose instead to scale many of our figure-of- 
merit variables with energy so as to remove most of the energy dependence. This allowed us 
to have a single set of event selection criteria spanning the entire LAT energy range. 


3.3.2. Event Energy Analysis 


The second step in the event analysis is selecting one energy estimate among those 
available for the event (see § 3.2.1). We apply a classification tree (CT) analysis (Breiman 


et al. 1984 Breiman 1996) to select which of the energy reconstruction methods is most 


likely to provide the best estimate of the event energy. 

Because of the numerous edges and gaps between the logs in the CAL, and because of 
the huge energy range and large FoV of the LAT, the quality of the energy reconstruction 
can vary significantly from event to event. Therefore, for each event we also apply a second 
CT analysis to estimate the probability that the estimated energy is within the nominal core 
of the energy dispersion. Specifically we define a scaled energy deviation, as described in 
more detail in § 7.1.1 from which most of the energy and angular dependence is factored 


(i.e. , the energy dispersion in the scaled deviation is largely energy- and angle-independent). 
We then train CTs that provide probability estimates this event is less than 2cr ( P 2a ) or 
3cr (P 3(T ) away from the most probable value of the energy dispersion. Finally, we define a 
reconstruction quality estimator Pe by combining P 2(J and P?, a \ 

Pe = s/PzA- ( 4 ) 

Large values of Pe indicate that the event is likely to be in the core of the energy dispersion, 
and so have an accurate energy estimate. 

In the Pass 7 analysis we did not use the energy estimates from the LH algorithm: 
by construction it is a binned estimator, and energy assignments tended to pile up at the 
bin boundaries (see § 7.2 for a more detailed discussion of the effects of removing the LH 


algorithm from consideration in the energy assignment). The removal of the LH energies 
causes a somewhat degraded energy resolution for those events where it would have been 
selected; we compensate for the loss of this energy estimator by requiring a slightly more 
stringent cut on the energy reconstruction quality when defining the standard event classes 
(e.g., see § 3.4.3). 


3.3.3. Ajialysis of the Event Direction 

The third step selects the measured direction of the incoming 7 ray from the available 
options. Those options are the directions as estimated from (i) the best track (ii) the best 
vertex and (iii and iv) the same two options using the centroid of the energy deposition in the 
CAL as an additional constraint on the direction of the incident 7 ray. Again, we use a CT 
analysis to combine the information about the event and determine which of the methods is 
most likely to provide the best direction measurement. 

As with the energy analysis, the quality of the direction reconstruction can vary signif- 
icantly from one event to the next. I11 this case, we have the additional complication that 
7 rays can convert at different heights in the TKR, giving us anywhere between 6 and 36 
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track position measurements. Therefore we use an additional CT analysis to estimate the 
probability P core that the reconstructed direction falls within the nominal 68% containment 
angle as defined by the simplified analytical model: 


C 68 (E) = 


\ 


Co • 


E Y 


100 MeV 




+ cf 


( 5 ) 


(where the values of the coefficients for front- and back-converting events are listed in Ta- 
ble [2]). The reader will notice the similarity with the functional expression used for the PSF 
prescaling described in § |6.1 


Conversion type 

Co [ c 

>i /-> r c 
C 1 

1 P 

Front 

3.3 

0.1 

0.78 

Back 

6.6 

0.2 

0.78 


Table 4: Coefficients of the analytical model for the 68% containment angle for the PSF 
from Eqn. ([5]) used for the multi-variate analysis used to evaluate the quality Pcore of the 
direction reconstruction. 


3.3.4 ■ Differences in Event Energy and Direction Analyses Between Pass 6 and Pass 7 

There is a particular subtlety to the event reconstruction that merits a brief discussion. 
Specifically, the different event reconstruction algorithms we use provide a set of choices for 
the best energy and direction of each event. The stages of the event-level analysis described 
in the previous two subsections select the algorithms that are more likely to provide the 
best estimates. Therefore, although the event reconstruction was unchanged except for 
the energy un-biasing, changes in the event-level analysis can result in individual events 
being assigned slightly different directions and/or energy estimates in Pass 7 with respect 
to Pass 6. Figure [13] illustrates such differences for the P7S0URCE events above 100 MeV 
that were also included in the P6_DIFFUSE 7-ray class. 


3.3.5. Rejection of Charged Particles in the Field of View 

The next step of the event analysis starts the process of classifying events as 7 rays 
or particle backgrounds, specifically identifying events for which evidence clearly indicates 
that a charged particle entered the LAT from inside the FoV. We refer to this stage as the 
Charged Particle in the Ficld-of-view (CPF) analysis. To accomplish this we first apply a 
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Fig. 13. — Event by event differences in the reconstructed energy (a) and direction (b) 
between Pass 6 and Pass 7 for a sample of P7S0URCE events above 100 MeV. The two 
histograms are based on a week of data, with a zenith angle cut at 100° to remove the 
bright Earth limb emission, typically viewed at large angle in LAT coordinates. In order to 
factor out the energy dependence of the PSF, the angular separation is measured in units 
of nominal 68% containment at the energy of the event. Note that in the vast majority of 
the cases the same direction is assigned in Pass 6 and Pass 7. On the other hand, the 
systematic offset in energy is due to the energy un-biasing introduced in Pass 7 (§ 3.2.1). 


series of filters that classify as background those events for which the best track has hits 
all the way to the edge of the TKR and (i) points to an active region of the ACD that has 
significant deposited energy, or (ii) points directly to less sensitive areas in the ACD, such 
as the corners or gaps between ACD tiles. 

After applying these filters we attempt to account for cases where the best reconstructed 
track does not represent the incoming particle well. This can happen for a variety of reasons, 
for example when the bias of the track finding algorithm toward the longest, straightest track 
causes us to incorrectly assign a backsplash particle from the CAL as the best track, or when 
the incoming particle passes largely in the gaps between the TKR towers and we incorrectly 
select some secondary particle as the best track. We classify as background the events for 
which the total energy deposited in the ACD is too large to be accounted for as backsplash 
from the CAL. It is important to note that this particular requirement loses the benefit of the 
ACD segmentation, so we apply a conservative version to all events, and a tighter version to 
events that have a track that passes within an energy dependent distance of significant energy 
deposition in the ACD. During extremely bright SFs this total ACD energy requirement can 
cause us to misclassify many 7 rays as CRs because of energy deposition in the ACD from 
SF X-rays (see Appendix [A]) . 
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All the events that are classified as charged particles at this stage are removed from 
further consideration as 7-ray candidates and are passed to separate analyses to identify 
various species of charged particles; see, e.g., Abdo et al. (2009e), Ackermann et al. (2010b), 


Ackermann et al. 

(2012c 

) and 

Ackermann et al. 

(2012d) 


Finally, we perform another CT analysis to combine all the available information into 
an estimate of the probability that the event is from charged particle backgrounds (Pcpf)- 
Although the performance of this selection depends on energy and angle, roughly speaking 
more than 95% of the background in the telemetered data is removed by means of these 
cuts, and about 10% of the 7-ray sample is lost. 


We note in passing that electrons and positrons cause electromagnetic showers that look 
extremely similar to 7-ray interactions, and the remaining stages of the event analysis are 
based on the topology (e.g., the shape, density and smoothness) of the energy deposition in 
the TKR (§ 3.3.6) and CAL (§ 3.3.7), and have very little additional discriminating power 


against such backgrounds. 


3.3.6. TKR Topology Analysis 

In the next stage of the event analysis we use information from the TKR to identify 
CR backgrounds that were not identified as such by the CPF analysis. These events are not 
immediately removed from the 7-ray analysis but only flagged to allow for removal from the 
higher purity event classes. In this part of the analysis we flag events saturating the energy 
deposition in the TKR planes with a pattern expected for heavily ionizing particles as well 
as events with very high energy deposition in the first hit layers of the TKR, which is a 
signature of a particle coming from the CAL that ranges out in the middle of the TKR. 

After these two flags are applied, we divide the events among 5 branches for evaluating 
their topologies in the TKR: the first one collects events with a vertex solution in the TKR, 
the remaining ones separate events with one track and events with many tracks, treating sep- 
arately events converting in the thin and thick section of the TKR. This is done because each 
of these different topologies has significantly different ratios of 7-ray signal-to-background 
contamination, and presents accordingly different difficulty for selecting 7 rays. For each 
branch we apply a different initial selection to remove many CR background events and 
then pass the remaining events to a CT analysis, which estimates the probability (Tf kr) the 
event comes from 7-ray interactions as opposed to charged particle backgrounds (and more 
specifically hadronic charged particle backgrounds). 

The variables used in the CT analysis were designed to emphasize differences between the 
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characteristics of 7 -ray conversions and hadronic CRs. This includes distinguishing between 
electromagnetic and hadronic showers by counting extra hits near the track, quantifying the 
complexity of the event in the TKR and how collimated the overall track structure is, looking 
how deep into the fiducial volume the event starts, and requiring that the ionization along 
the track is consistent with e ± pairs. 


3.3.7. CAL Topology Analysis 


Next we apply a CAL topology analysis that is similar in design to the TKR topology 
analysis. The first part of this analysis consists of few general cuts flagging events coming 
from the bottom and sides of the LAT. As with the TKR topology analysis, we then split 
the events into five branches depending on the topology; for each branch we apply a cut 
and a CT stage. Again, CTs are trained to select 7 rays versus hadronic charged particle 
backgrounds, and provide an estimate of the probability that the event is a 7 ray(PcAL)- 


The variables used in the CAL topology analysis are not only CAL derived quantities, 
but also involve comparisons between what is observed in the CAL and what is seen in the 
TKR. Among the important discriminants are how well the track solution points toward 
the CAL energy centroid, how well the direction derived from the CAL (via the moments 
analysis, see § 3.2.1) agrees with the track direction, the root mean square (RMS) width of 


the CAL shower, and the ratio of the highest energy in a single crystal to the total energy 
in all crystals. 


3.3.8. Event Classification 


After the analyses related to the three LAT subsystems, the last stage of our event 
analysis applies a final CT analysis utilizing all available information, notably including the 
outputs of the CT analyses from previous phases of the event analysis. This second-order 
CT analysis is particularly important for defining event classes with high 7 -ray purity, as 
discussed in 


3.4 


At this point we have a number of specifiers of event reconstruction and classification 
quality on an event-by-event basis: 


• energy reconstruction quality Pe (see 

• direction reconstruction quality P cor e 


3.3.2); 


see 


3.3.3); 
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7 -ray probabilities from the CPF analysis (Pcpf, see § 3.3.5) and the TKR (Ptkr, see 


3.3.6) and CAL (Pcal, see § 3.3.7) topology analysis; 


overall y-ray probability P a u from the final classification step of the event analysis. 


These are the basic ingredients for defining the standard y-ray classes described in § 3.4 


3.4. Standard Event Classes for y rays 

Making use of the classification quantities described in the preceding subsections we 
define event classes optimized for a range of astrophysical source analyses. We note that 
these event classes are nested: each succeeding selection is a strict subset of the previous 
ones. The relative selection efficiencies for the event classes as well as previous stages of the 
event selection process are shown in Figure [l4| 



Fig. 14. — LAT A eff integrated over the FoV as a function of energy at successive stages of 
the event filtering as estimated with simulated data. Since we require a direction and energy 
to use a y ray for science analysis we consider only events with at least one track found and 
that pass the fiducial cuts (see § 3.4.1). 


3.4.1. Track-Finding and Fiducial Cuts 

For the definition of the standard y-ray classes we confine ourselves to events that have 
some chance to be useful for most standard science analysis, namely events for which we 
have enough information to reconstruct an energy and a direction. Therefore, we require 
that the event has a reconstructed track to allow for a direction estimation. Furthermore, 


50 


we require that the track points into the CAL and crosses at least 4 radiation lengths of 
material in the CAL and that at least 5 MeV of energy is deposited in the CAL. This second 
requirement reduces the fiducial volume of the detector by rejecting off-axis 7 rays that pass 
near the top of the TKR and miss the CAL. Events with no tracks in the TKR or with less 
energy deposited in the CAL are not considered further as candidate 7 rays. The remaining 
data set (i.e., the events that are passed along as potential candidates for the standard 7-ray 
classes) is still composed almost entirely of background events. 


We note that these cuts remove from consideration two classes of events that might be 
useful for specific, non-standard, analyses. The first class consists of events that deposit all 
their energy in the TKR, either because they range out before they reach the CAL or simply 
because they miss the CAL entirely. In general the energy resolution is much poorer for these 
events, though at low energies (< 100 MeV) the CAL does not improve the energy resolution 
significantly. These events have been used effectively in the analysis of both Gamma Ray 
Bursts (GRBs) (Pelassa et al. 2010) and SFs (Ackermann et al. 2012a). The second class 
consists of events that do not have reconstructed tracks, but have enough information in the 
CAL to derive an estimate of the event direction (though without the TKR information the 
angular resolution for these events is highly degraded). These events can occur because the 
7 ray entered the LAT at a large incidence angle and missed most of the TKR, or simply 
because the 7 ray passed through the TKR without converting. 


3.4.2. P7TRANSIENT Class Selection 

For the analysis of brief transient sources (e.g., GRBs) a high purity is not required as 
the time selection itself limits the amount of background counts in the ROI. Accordingly, we 
define a P7TRANSIENT event class with only a few cuts, with the aim of achieving a residual 
background rate of a few Hz while maintaining a large efficiency for 7 rays. The list of cuts 
is short: 

1. the event must pass minimal cuts on the quality estimators Pe and P CO re(-pE > 0. and 
Pcore > 0 .); 

2. the event energy after all corrections must be greater than 10 MeV; 

3. the event energy after all corrections must not amount to more than 5 times the energy 
deposited in the CAL; 

4. the CPF analysis must not flag the event as a charged particle (Pcpf > 0.1); 

5. the event must pass a relatively loose cut on P a n(P a ii > 0.2). 


51 


It is worth noting that although the cuts are optimized somewhat differently in Pass 6 : 
a similar series of cuts are applied to define the P 6 _TRANSIENT event class, which is roughly 
equivalent to the P7TRANSIENT event class. 


3.4.3. P7S0URCE Class Selection 


We have several additional considerations when defining the P7S0URCE event class, which 
is intended for the analysis of point sources. While the 7 rays from a point source are clustered 
on the sky near the true source position, residual CR background can be modeled as an 
isotropic component and be accounted for in the analysis of the spectrum or position of the 
source. Although the exact numbers depend on the details of the PSF and the spectrum of the 
source, in general we require a background rate of less than ~ 1 Hz in the LAT FoV to ensure 
that we maintain a high enough signal-to-background ratio so that this has little impact 
on source detection and characterization. Furthermore, in contrast to studies of transient 
sources, for sources with an integral energy flux above 100 MeV of ~ 10~ n erg cm -2 s ^ 1 
the precision of spectral studies is often limited by systematic uncertainties below 1 GeV 
, increasing the importance of the precision and accuracy of the event energy and 
direction reconstruction. 


(see § 8.3 


Additionally, we note that for the P7S0URCE event class and both higher-purity event 
classes (P7CLEAN and P7ULTRACLEAN) we developed and optimized the cuts using on-orbit 
data samples (particularly the clean and dirty samples described in § 3.6.4) as well as MC 


simulations. In particular, we performed a comparison of many event quality parameters 
and reconstruction variables between the clean and dirty samples to identify characteristics 
of the CR background, which is greatly enhanced in the dirty sample, and devised selection 
criteria to remove it. Some examples of this procedure are described in more detail in § |4.3 


An important limitation of this technique is obviously that one needs a very low ratio 
of residual CR background to 7 rays in the clean sample to obtain large contrast between 
the clean and dirty samples. It can only be applied to data samples which already have 
a 7 -ray-to-residual-background ratio that is comparable to 1. Therefore, when optimizing 
these selections we only used events that are part of the P7TRANSIENT class (see§ 3.4.1) and 
additionally have a high value of P a n estimator introduced in § 3.3.8 


The P7S0URCE event class is defined starting from the P7TRANSIENT event class and 
includes tighter cuts on many of the same quantities. Specifically, the P7S0URCE event class 
selection requires that: 


1 . the event must not have been flagged as background based on topological analysis of the 


52 


CAL (§ 3.3.6) and TKR (§|3.3.6) reconstruction, and must pass a tighter cut on P E ; the cut 


on P E does reject a noticeable number of events, though by definition these events have less 
accurate energy estimates on average (P E > 0.1 or P E > 0.3 for events which originally used 
the LH energy estimate); 


2. the event must pass cuts to reject MIPs based on the event topology in the CAL and the 
energy deposited in the TKR; 


3. the event must pass a set of cuts on the agreement between the TKR and CAL direction 
reconstruction; 


4. the event must pass a tighter, energy dependent, cut on P core 


P core > max(0.025, 0.025 + 0.175(3. - log 10 (P/lMeV)); 


5. the event must pass a tighter cut on P a n 

P a ii > max(0.7, 0.996 — 1.4 x 10 _4 (max(5.4 — log 10 (P/lMeV), 0) 5 ' 1 2 3 )). 

As with the TRANSIENT event classes, it is worth noting that although the cuts are 
optimized somewhat differently in Pass 6 a similar series of cuts are applied to define the 
P6_DIFFUSE event class, which is roughly equivalent to the P7S0URCE selection. 


3. 4 - 4 - P7CLEAN Class Selection 

For the analysis of diffuse y-ray emission, we need to reduce the background contami- 
nation to a level of about ~ 0.1 Hz across the LAT FoV, to keep it below the extragalactic 
y-ray background at all energies. For comparison, the total Galactic diffuse contribution is 
~ 1 Hz, depending on where the LAT is pointing, though most of that is localized along the 
Galactic plane. 

The selection of P7CLEAN class events starts from the P7S0URCE event class and includes 
the following additional cuts: 


1. the event must pass a series of cuts designed to reject specific backgrounds such as CRs 
that passed through the mounting holes in the ACD, or the gaps along the corners of the 
ACD, with minimal costs to the y-ray efficiency 


see 


4.2); 


2. the event must pass cuts on the topology of the event in the CAL designed to remove 

hadronic CRs. 
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As with the two previous event classes, a similar series of cuts was applied to define the 
P6_DATACLEAN event class, which is roughly equivalent to the P7CLEAN event class. 


3.4.5. P7ULTRACLEAN Class Selection 

For the analysis of extragalactic diffuse 7 -ray emission we need to reduce the background 
contamination even further below the extragalactic 7 -ray background rate to avoid introduc- 
ing artificial spectral features. As illustrated in Figure [28] the residual contamination for the 
P7ULTRACLEAN class is ~ 40% lower than that of the P7CLEAN class around 100 MeV (the 
residual levels becoming more similar to each other as the energy increases and becoming 
the same above 10 GeV). 

The selection of P7ULTRACLEAN class events is relatively simple, consisting of a tighter, 
energy-dependent cut on P a n: 


P a ii > 0.996 — 0.0394(max(3.26 — logio(E/l MeV), 0 ) L ' 8 (Front) 

P a u > 0.996 — 0.006(max(4.0 — log w (E /I MeV), O) 3 0 (Back). 

3.5. Publicly Released Data 

At the time of this writing the LAT team has published results and released data 
for both the Pass 6 and Pass 7 event analyses, as well as several event classes and the 
associated IRFs for each iteration of the event analysis. Furthermore, as our understanding 
of the LAT has improved, we have updated the IRFs for particular event classes. Table [5] 
summarizes these data sets and associated IRFs. This paper will focus in particular on the 
performance of the P7S0URCE event class and the validation of the associated P7S0URCE_V6 
IRFs since this is the data set the LAT team currently recommends for most analyses. 
Of the other Pass 7 event classes, P7TRANSIENT is recommended for the analysis of short 
(< 1000 s) transient events such as GRBs, P7CLEAN is recommended for analyses requiring 
low CR background contamination, such as the study of large-scale diffuse emission, and 
P7ULTRACLEAN is recommended when CR background contamination must be minimized, 
even at the cost of some loss of effective area, such as when studying the extragalactic 
background. Accordingly, all plots, figures and tables will be made with the P7S0URCE event 
sample and the P7S0URCE_V6 set of IRFs unless stated otherwise. 

One other very important point is that the excellent stability of the LAT subsystems (see 
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2.1 ) means that changes in the instrument performance over time are not a consideration in 


defining the event analyses or IRFs. The small changes in performance at the subsystem level 


are easily addressed by calibrations applied during the event reconstruction procedure (Abdo 


et ah 2009f). Accordingly, to date the LAT team is able to produce IRFs that are valid for 


the entire mission. This in turn simplifies the data analysis task by removing the need to 
split the LAT data by time range. 


Event Class 

Pass 7 IRF Set 

Pass 6 Counterpart 

Pass 6 IRF Set 

P7TRANSIENT 

P7TRANSIENT_V6 

P6.TRANSIENT 

P6_V1 .TRANSIENT 
P6_V3_TRANSIENT 

P7S0URCE 

P7S0URCE_V6 a 

P6_DIFFUSE 

P6_V1_DIFFUSE 
P6_V3_DIFFUSE 
P6_V11 .DIFFUSE 

P7CLEAN 

P7CLEAN_V6 a 

P6_DATACLEAN 

P6.V3.DATACLEAN 

P7ULTRACLEAN 

P7ULTRACLEAN_V6 

- 

- 


Table 5: Publicly released event selections and IRFs. Note the slight change in naming 
conventions between Pass 6 and Pass 7. The Pass 7 naming convention emphasize the 
point that we may release multiple IRFs for the same event class as we improve the IRFs 
and background rejection. 


“We have also released P7S0URCE.V6MC and P7CLEAN.V6MC IRF sets, which feature a MC-based PSF that 
includes 9 dependence that we have used when we need to minimize the potential of instrument-induced 
variability (see §[6]). 


3.6. Calibration Sources and Background Subtraction Methods 

Because of the complexity of the LAT and of the physics simulations of particle inter- 
actions we cannot expect the MC to perfectly reproduce the flight data. For this reason we 
have developed validation data sets for the IRFs. Although no astrophysical source has per- 
fectly known properties, in practice there are several sources for which accurate background 
subtraction allows to extract a clean y-ray sample that we can use to validate the MC predic- 
tions. Table [6] summarizes these calibration sources and associated background subtraction 
techniques and MC samples. The remainder of the section discusses the particulars of these 
samples. 


Table 6. Summary table of calibration data samples. 
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e In the case of the Galactic ridge we cannot subtract background accurately enough for detailed quantitative validations; 
however, we can distinguish between regions of higher and lower y-ray-to-CR ratios. 
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3.6.1. The Vela Pulsar 


The Vela pulsar (PSR J0835— 4510) has the largest integral flux > 100 MeV of any 
7-ray source, and has been well studied in the LAT energy range (Abdo et al. 2009c, 2010h). 
Furthermore, the pulsed nature of high-energy y-ray emission gives us an independent control 
on the background. In fact, off-pulse y-ray emission is almost entirely absent. These factors 
combine to make the Vela pulsar an almost ideal calibration source. Unfortunately, the 
spectrum of Vela cuts off at ~ 3 GeV and Vela is nearly undetectable above 30 GeV. 

The selection criteria we use to define our Vela calibration samples are listed in Ta- 
ble § Specific calibration samples used for particular studies may include additional re- 
quirements. For example, the “P7TRANSIENT Vela calibration sample” includes all events in 
the P7TRANSIENT event class that pass the Vela calibration sample criteria. We used the 
TEMP02 packag^] (Hobbs et al. 2006) and a pulsar timing modej^] derived from data taken 
with the Parkes radio telescope ( Abdo et al.||2009c Weltevrede et al.||2010[ ) to assign a phase 
to each y ray. 



Fig. 15. — Phase-folded times for events in the P7TRANSIENT Vela calibration sample, which 
includes events between 30 MeV and 10 GeV. The signal and background regions, as defined 
in the background subtraction analysis, are highlighted. 

We can achieve excellent statistical background subtraction for any distribution (i.e. , 
spectrum, spatial distribution, any discriminating variable used in the event classification) 
by subtracting the distribution of off-pulse events (defined as phases in the range [0.7, 1.0]) 
from the distribution of on-peak events (defined as phases in [0.125, 0.175] U [0.5125, 0.6125]). 


^ http: / / www.atnf.csiro.au/research/pulsar /tempo2 / 
f ’ht.tp: / /fermi. gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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Figure [15] shows the distribution of phases of 7 rays in the P7TRANSIENT Vela calibration 
sample, including our standard on-peak and off-pulse regions. (Note that the PSF analysis 
6.2.1| uses slightly different definitions of on-peak ([0.12, 0.17] U [0.52, 0.57]) and off-pulse 
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[0.8, 1.0] regions; the slight difference in definition causes no significant change in the results). 

Admittedly, all bright pulsars are potential candidates; nonetheless for this analysis the 
consequent increase in event statistics does not warrant the procedural complication required 
to deal with a stacked sample of pulsars (see also § 6.2). 


3.6.2. Bright Active Galactic Nuclei 


At 1 GeV the 95% containment radius for the P7S0URCE event class is ~ 1.4° for front- 
converting events and ~ 2.5° for back-converting events (see § [6j) . Given the density of bright 
7-ray sources in the sky, above this energy the instrument PSF becomes narrow enough that 
we can use the angular distance a between a 7 ray and the nearest celestial 7-ray source as a 
good discriminator for background subtraction, particularly at high Galactic latitudes where 
there are fewer sources and the interstellar diffuse emission is less pronounced. Unfortunately, 
no single source is bright enough to provide adequate statistics to serve as a good calibrator. 
However, by considering 7 rays from a sample of bright and/or hard spectrum Active Galactic 
Nuclei (AGN) that are isolated from other hard sources we can create a good calibration 
sample for energies from 1 GeV to 100 GeV. 


Table [7] lists the AGN that we use here, Figure [16] shows their positions, and Figure [T7[ 
shows a comparison of their spectral index (r) and integral 7-ray flux between 1 GeV and 
100 GeV (P35) to other AGN in the second LAT source catalog (2FGL; Nolan et al. 2012). 
Note that different but overlapping sets of AGN are used for A e ff-related studies and PSF- 
related studies. More information about the source properties can be found in Abdoet al. 


(2010a 

) and 

Nolan et al. 

(2012) 


The selection criteria we use to define AGN calibration samples are listed in Table [6] 
As with the Vela calibration samples, specific calibration samples used for particular studies 
may include additional requirements (e.g., the P7TRANSIENT AGN calibration sample). 


To use this calibration sample to perform background subtraction we define signal and 
background regions in terms of the angular distance a between the 7 ray and the closest 
AGN. Specifically, we use a < 0.5° for the signal region and the annulus 3.87288° < a < 4° 
for the background region. These ranges are chosen so that the background region contains 
four times the solid angle of the signal region. We adopted the separation between the signal 
and background regions in order to minimize signal 7 rays from the tails of the PSF leaking 
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Source 

i n 

b n 

2FGL 

Used 

KUV 00311-1938 

94.17 

-81.21 

J0033. 5-1921 

PSF 

PKS 0118-272 

213.58 

-83.53 

J0120. 4-2700 

PSF 

B3 0133+388 

132.43 

-22.95 

J0136. 5+3905 

PSF 

OC 457 

130.79 

-14.32 

J0136. 9+4751 

A eff PSF 

PKS 0208-512 

276.12 

-61.76 

J0210. 7-5102 

^eff 

PKS 0215+015 

162.20 

-54.41 

J0217. 9+0143 

PSF 

S4 0218+35 

142.60 

-23.49 

J0221. 0+3555 

PSF 

3C 66A 

140.14 

-16.77 

J0222. 6+4302 

^eff 

AO 0235+164 

156.78 

-39.10 

J0238. 7+1637 

PSF 

PKS 0250-225 

209.72 

-62.10 

J0252. 7-2218 

^eff 

PKS 0301-243 

214.64 

-60.17 

J0303.4— 2407 

PSF 

NGC 1275 

150.59 

-13.25 

J0319. 8+4130 

PSF 

PMN J0334— 3725 

240.22 

-54.36 

J0334.3— 3728 

PSF 

PKS 0420-01 

195.28 

-33.15 

J0423. 2-0120 

PSF 

PKS 0426-380 

240.70 

-43.62 

J0428.6— 3756 

PSF 

MG2 J043337+2905 

170.52 

-12.62 

J0433. 5+2905 

PSF 

PKS 0440-00 

197.21 

-28.44 

J0442. 7-0017 

PSF 

PKS 0447-439 

248.81 

-39.91 

J0449.4— 4350 

Aeff PSF 

PKS 0454-234 

223.73 

-34.90 

J0457. 0—2325 

Aeff PSF 

1ES 0502+675 

143.80 

15.90 

J0508. 0+6737 

PSF 

TXS 0506+056 

195.40 

-19.62 

J0509. 4+0542 

PSF 

PKS 0537-441 

250.08 

-31.09 

J0538.8— 4405 

Aeff PSF 

TXS 0628-240 

232.68 

-15.00 

J0630.9— 2406 

PSF 

B3 0650+453 

171.20 

19.36 

J0654. 2+4514 

^eff 

PKS 0700-661 

276.77 

-23.78 

J0700.3— 6611 

PSF 

MG2 J071354+1934 

197.68 

13.61 

J0714. 0+1933 

^eff 

B2 0716+33 

185.06 

19.85 

J0719. 3+3306 

A e ff PSF 

S5 0716+71 

143.97 

28.02 

J0721. 9+7120 

^eff 

PKS 0727-11 

227.77 

3.13 

J0730. 2-1141 

PSF 

PKS 0735+17 

201.85 

18.06 

J0738. 0+1742 

PSF 

PKS 0805-07 

229.04 

13.16 

J0808. 2-0750 

PSF 

S4 0814+42 

178.21 

33.41 

J0818. 2+4223 

PSF 

PKS 0823-223 

243.97 

8.92 

J0825. 9-2229 

PSF 

S4 0917+44 

175.70 

44.81 

J0920. 9+4441 

PSF 

4C +55.17 

158.59 

47.94 

J0957. 7+5522 

PSF 

1H 1013+498 

165.53 

52.73 

J1015. 1+4925 

Aeff PSF 

4C +01.28 

251.50 

52.77 

J1058. 4+0133 

PSF 

TXS 1055+567 

149.57 

54.42 

J1058. 6+5628 

PSF 

Mkn 421 

179.82 

65.03 

J1104. 4+3812 

Aeff PSF 

Ton 599 

199.41 

78.37 

J1159. 5+2914 

A e ff PSF 

W Comae 

201.69 

83.28 

J1221. 4+2814 

PSF 

4C +21.35 

255.07 

81.66 

J1224. 9+2122 

PSF 

PKS 1244-255 

301.60 

37.08 

J1246. 7-2546 

Aeff PSF 

PG 1246+586 

123.74 

58.77 

J1248. 2+5820 

A e ff PSF 

S4 1250+53 

122.36 

64.08 

J1253. 1+5302 

PSF 

3C 279 

305.10 

57.06 

J1256. 1-0547 

PSF 

OP 313 

85.59 

83.29 

J1310. 6+3222 

-A-etf 

GB 1310+487 

113.32 

68.25 

J1312. 8+4828 

PSF 

PMN J1344— 1723 

320.48 

43.67 

J1344. 2-1723 

PSF 

PKS 1424+240 

29.48 

68.20 

J1427. 0+2347 

PSF 

PKS 1440-389 

325.64 

18.72 

J1443.9— 3908 

PSF 

PKS 1454-354 

329.89 

20.52 

J1457.4— 3540 

PSF 

PKS 1502+106 

11.37 

54.58 

J1504. 3+1029 

PSF 

AP Librae 

340.70 

27.58 

J1517.7— 2421 

PSF 

B2 1520+31 

50.18 

57.02 

J1522. 1+3144 

A e ff PSF 

GB6 J1542+6129 

95.38 

45.40 

J1542. 9+6129 

PSF 

PG 1553+113 

21.92 

43.95 

J1555. 7+1111 

A e ff PSF 

1H 1720+117 

34.11 

24.47 

J1725. 0+1151 

PSF 

B2 1732+38A 

64.04 

31.02 

J1734. 3+3858 

^eff 

PMN J1802— 3940 

352.44 

-8.42 

J1802. 6-3940 

PSF 

PKS 1830-211 

12.15 

-5.72 

J1833.6— 2104 

^eff 

S4 1849+67 

97.50 

25.03 

J1849. 4+6706 

^eff 

TXS 1902+556 

85.96 

20.51 

J1903. 3+5539 

PSF 

1H 1914-194 

18.24 

-14.30 

J1917. 6-1921 

PSF 

TXS 1920-211 

17.17 

-16.26 

J1923. 5-2105 

PSF 

1ES 1959+650 

98.02 

17.67 

J2000. 0+6509 

A e ff PSF 

PKS 2005-489 

350.37 

-32.61 

J2009. 5-4850 

PSF 

PKS 2023-07 

36.89 

-24.39 

J2025.6— 0736 

PSF 

PKS 2052-47 

352.58 

-40.38 

J2056. 2-4715 

PSF 

MH 2136-428 

358.29 

-48.32 

J2139. 3-4236 

PSF 

PKS 2155-304 

17.74 

-52.24 

J2158.8— 3013 

PSF 

BL Lacertae 

92.60 

-10.46 

J2202. 8+4216 

Aeff PSF 

PKS 2201+171 

75.68 

-29.63 

J2203. 4+1726 

PSF 

B2 2234+28A 

90.12 

-25.66 

J2236. 4+2828 

A e ff 

3C 454.3 

86.12 

-38.18 

J2253. 9+1609 

PSF 

PKS 2326-502 

332.00 

-62.30 

J2329. 2-4956 

PSF 


Tabic 7: List of the Active Galactic Nuclei used in the AGN calibration samples for A e g and 
PSF studies. 
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Fig. 16. — Positions of the AGN in the calibration samples, shown in a Hammer- Aitoff 
projection in Galactic coordinates. Circles mark the AGN used for A eS studies and crosses 
those used in PSF studies. 



r 

Fig. 17. — Integral 7 -ray fluxes between 1 GeV and 100 GeV (F 35 ) and spectral indices (r) 
of the AGN in the calibration samples. Circles mark the AGN used for A e g- studies and 
crosses those used in PSF studies. The dots mark all the other sources associated with AGN 
in the 2FGL catalog. 

into the background region. Figure [18] shows the squares of the angular separations between 
7 rays and the nearest bright AGN for all 7 rays in the P7TRANSIENT event class in 6° regions 
around the 25 AGN listed as the A e g calibration sample in Table [7] including our definitions 
of source and background regions. 


3.6.3. The Earth Limb 

The Earth’s atmosphere is a very bright 7 -ray source. Furthermore, at energies above 
a few GeV the 7 -ray flux seen by the LAT is dominated by 7 rays from the interactions 
of primary CR protons with the upper atmosphere. This consideration, together with the 
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Fig. 18. — Square of the angular separation between reconstructed 7-ray directions and 
the AGN for events in the P7TRANSIENT AGN calibration sample, which includes 7 rays in 
the energy range 1 GeV to 100 GeV. The signal and background regions, as defined in the 
background subtraction analysis, are highlighted. 


narrowness of the PSF at energies > 10 GeV, causes the Earth limb to appear as a very 
bright and sharp feature, which provides an excellent calibration source. Furthermore, we 
have selected data from 200 orbits during which the LAT was pointed near the Earth limb 
as the basis of the Earth limb calibration sample. The selection criteria we use to define 
the Earth limb calibration samples are listed in Table [6j When using the Earth limb as 
a calibration source we generally limit the energy range to energies > 10 GeV, primarily 
because at lower energies orbital variations in the geomagnetic field significantly affect the 
7-ray fluxes (however see § 7.5 for a counter-example). 


For this calibration source, we define our signal region as 111.1002° < 9 Z < 112.9545° 
and background regions just above and below the limb: 108.6629° < 9 Z < 109.5725° and 
114.5193° < 9 Z < 115.4675°. Note that these ranges are defined to give the same solid angle 
(0.067T sr) in the signal and background regions. 


Figure 19 shows the zenith angle distribution for all 7 rays in the P7TRANSIENT event 
class for the Earth limb calibration sample. 


3.6.4. The Galactic Ridge 

At energies above ~ 30 GeV no single source provides enough 7 rays for a good com- 
parison between flight data and MG. However, the combination of bright Galactic sources 
and Galactic diffuse backgrounds means that there is a very large excess of 7 rays coming 


61 



e z [°] 


Fig. 19. — Angle with respect to the zenith for events in the P7TRANSIENT Earth limb calibra- 
tion sample, which includes 7 rays with energies above 10 GeV. The signal and background 
regions, as defined in the background subtraction analysis, are highlighted. 


from the Galactic plane relative to high Galactic latitudes. 

The intensity of the 7-ray emission at low latitudes in the inner Galaxy is more than 
an order of magnitude greater than at high latitudes in the outer Galaxy. In contrast, the 
intensity of the CR background is approximately isotropic for observation periods longer 
than the 53.4-day orbital precession period. 

Unfortunately, since the Galaxy extends over much of the sky, and since the data set 
consists of several thousand orbits it is not practical to disentangle the variations of exposure 
from the spatial variations in Galactic diffuse emission without relying on detailed modeling 
of the Galactic diffuse emission. Accordingly, we use the Galactic ridge primarily when we are 
developing our event selections, rather than for precise calibration of the LAT performance. 


Specifically, we developed the event classes that require a high 7-ray purity, i.e. , the 
event classes used in the analysis of celestial point sources and diffuse emission (see § 3.4), in 


part by tuning our selection criteria to maximize the contrast between regions in the bright 
Galactic plane and at high Galactic latitudes. This helped to mitigate the risk that insuffi- 
cient statistics of the MC training sample or limited accuracy of the MC description of the 
geometry of the detector and the particle interactions in the LAT limited the discriminating 
power and accuracy of the event classification analysis. 


The selection criteria we use to define the Galactic ridge calibrations samples are listed 
in Table [6] In particular, we use the region (|6| < 1.5° , —40° < l < 50°, which was selected 
to maximize the total 7-ray flux) to define a clean data sample and the region (\b\ > 50°, 
90° < l < 270°) to define a dirty data sample. The ratio of 7 rays to CR background in 
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the clean region is more than an order of magnitude higher than the same ratio in the dirty 
region. Figure [20] shows the count maps for the P7S0URCE samples for both regions. 

Furthermore, to give a sense of the statistics of these samples at high energies, Figure [21] 
shows the Galactic latitude distribution for all 7 rays in the P7TRANSIENT event class with 
energies above 17783 MeV (i.e. , log 10 (E/1 MeV) > 4.25). 



ur 


O 
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Fig. 20. — Definitions of the clean and dirty regions, showing the counts in both regions 
in a Hannner-Aitoff projection. This figure uses the data in the P7S0URCE Galactic ridge 
calibration sample. 



Fig. 21. — Sine of Galactic latitude for events above 17783 MeV in the P7TRANSIENT Galactic 
ridge calibration sample. The Galactic latitude selections for the clean and dirty regions are 
highlighted. Note that the definition of the clean and dirty regions also include selections on 
Galactic longitude. 


3.6.5. Summary of Astrophysical Calibration Sources 

As we will see in the next sections, the IRFs depend heavily on 6 (and, to a lesser extent, 
on 0). Therefore, any detailed comparison between flight data and MG simulations must 
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account for the distribution of observing profile, particularly t ohs (9). flow best to account 
for the observing profile depends on the particulars of the calibration samples. 

For any point source, the observing profile is determined by the position of the source, 
the rocking angle of the LAT and the amount of time spent in survey mode relative to 
pointed observations. Figure [22] shows the observing profile for Vela for the first two years of 
the mission. Rather than produce a dedicated large statistics MC sample for Vela, we re-use 
our allGamma MC sample, re-weighting the events in that sample so as to match the Vela 
observing profile. 



Fig. 22. — Vela observing profile: starting from the top, the curves show the accumulated 
time as a function of cos 6 for the first two years of the mission, the accumulated observing 
time (accounting for the dead time when the LAT triggers), the time during which Vela was 
less than 100° from the zenith, and the time during which the entire 15° ROl around Vela 
was less than 100° from the zenith (or, equivalently, that Vela was less than 85° from the 
zenith). The shaded region corresponds to the area outside the LAT FoV. 


Similarly, we re-weight the allGamma MC to match the summed observing profiles of 

Unfortunately, since AGN are 


all of the AGN our sample, which is shown in Figure 23 


intrinsically variable, and since the AGN in this sample span a range of fluxes, this re- 
weighting technique will not work as well with this sample. On the other hand, by taking a 
large set of AGN, we reduce the bias due to the variability of any one particular source. In 
broad terms, our re-weighted MC sample reproduced the 9 distribution of the AGN sample 
to better than 2% (see § 5.5). Finally, we note that since the PSF is narrower above 1 GeV, 


and the ROI around each AGN is only 6° we do not apply the ROI-based 9 Z cut when 
building the AGN calibration samples. 


Since the Earth limb is a spatially extended source, we cannot apply the re-weighting 
technique we used for the Vela and AGN samples to account for the observing profile. On 
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Fig. 23. — Sum of the observing profiles for the AGN sample: starting from the top, the 
curves show the sum of the accumulated time as a function of cos 6 for the first two years of 
the mission, the sum of the accumulated live time, and the sum of the accumulated live times 
during which each AGN was less than 100° from the zenith. The shaded region corresponds 
to the area outside the LAT FoV. 


the other hand, since the data set consists of only 200 orbits, and the Earth limb emission is 
well understood above 10 GeV we can produce a MC simulation of the Earth limb emission 
for those orbits and compare it with the flight data (see 


2.5.2). 


Finally, Figure [24] shows the statistics available for each of the samples. This shows that 
the calibration sources span most of the LAT energy range, certainly from 30 MeV up to at 
least 100 GeV. 



Fig. 24. — Excess counts in the signal regions as a function of energy for Vela, stacked AGN, 
and Earth limb calibration samples for the P7TRANSIENT event class. The horizontal error 
bars indicate the energy binning used with the different samples. 
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BACKGROUND CONTAMINATION 


In this section we discuss the residual particle backgrounds, methods to estimate the 
contamination of LAT 7-ray samples by these backgrounds and how to treat such contami- 
nation in high-level analysis. In this context we define particle backgrounds as all events that 
are classified as 7 rays in a given LAT event class but originate from CRs or the interactions 
of CRs in the Earth’s atmosphere. Therefore, the particle backgrounds include both charged 
and neutral particles — including secondary 7 rays. 


The LAT background model is described in detail in § 2.5.4 We focus here on describing 
the particle background contamination in the high-purity event classes, i.e. , the ones used 
for single source, source population and diffuse emission analysis (P7S0URCE and above). We 
also focus on average background contamination for long (few months or longer) observation 
periods. The CR-induced particle background is extremely variable throughout the orbit 
of the LAT; therefore estimates of particle backgrounds for brief transient sources must be 
derived from a dedicated analysis of the periods of interest (e.g., as done in Abdo et al. 


2009d). 


4.1. Residual Background Contamination in Monte Carlo Simulations 

We can estimate the residual background in the various event classes by propagating the 
LAT background model (see § 2.5.4) through the full Geant4-based LAT detector simulation 
chain and applying the event classification analysis on the simulated data (see 
and 


2.2, 8 2.4 


2.5). 


In comparison to the pre-launch particle background model shown in Figure 10 we 
have implemented substantial improvements in our model of the primary CR protons and 
electrons. The effects of the geomagnetic cutoff on the directional and energy dependence 
of the primary CR flux in the pre-launch were based on a dipole approximation of the 
geomagnetic field. Currently, we simulate an isotropic flux of CRs, and trace the trajectory 
of each particle backward through the geomagnetic magnetic field, eliminating any particles 
that intersect the Earth or lower atmosphere. We are using the current version of the 
International Geomagnetic Reference Field (IGRF-11), a high order multipole model of the 


geomagnetic field ( 

Finlay et al. 

2010 

), and the publicly available trajectory tracing code 

described in 

Smart & Shea 

(2005 

. We oversample the primary proton and electron spectrum 


at high energies to obtain sufficient statistics up to ~ 600 GeV in reconstructed energy, and 
obtain rate predictions by appropriately deweighting events at the oversampled energies. 


A total of 2.2 x 10 12 primary protons and 1.6 x 10 8 primary electrons were generated. 


In 
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addition, the equivalent of 80 ks (~ 1 day) of instrument live time of background events from 
CR secondaries produced in the Earth atmosphere were simulated according to the spatial 
and spectral distributions in the pre-launch LAT particle background model (i.e., with no 
trajectory tracing). This intensive simulation effort was used to determine the CR-induced 
background between 100 MeV and ~ 600 GeV in the P7S0URCE, P7CLEAN, and P7ULTRACLEAN 
event classes. 


However, these simulations still have important shortcomings. For example, the model- 
ing of inelastic interactions of alpha particles and heavier ions does not match our observa- 
tions (see further discussion in § 4.2, especially Fig. 25 b). The particle background model for 


secondaries produced in CR interactions in the atmosphere is derived from measurements 
from satellites with different orbits and during different parts of the Solar cycle, and has 
been projected to the LAT orbit. Furthermore, due to our required background suppression 
factors of up to 10 6 , even small inaccuracies in our simulation of particle interactions with 
the LAT can potentially lead to large discrepancies between predicted and actual charged 
particle background rates. 

As a consequence we take several measures to account for any such shortcomings and 
minimize possible discrepancies. First, we compare a region on the sky with a high 7 ray 
to residual background ratio to a region with a low 7 ray to residual background ratio to 
isolate those contributions caused by the accumulation of background in the latter data set; 
and we compare those contributions to the predictions from the background simulations (see 


4.2). Second, we compare several key event parameters between data and simulation for 


high-purity 7-ray samples to crosscheck the agreement between data and simulation, and we 
slightly adjust the intensities of our CR background model based on a fit of the shapes of 
two key event classification variables with signal and background component (see § 4.3). 


4.2. Estimating and Reducing Residual Background Using Test Samples 


We already described the usage of data samples with different 7-ray-to-background ratios 
to define high-purity event classes in § |3.6.4| and we detailed the resulting event classification 
cuts in § 3.4.3 § 3.4.4 and § 3.4.5 A key part of the development of the event classes was 


using the same clean and dirty samples to search for and eliminate residual background that 
is either not simulated with sufficient accuracy or has passed the multi-variate event selection 
that was trained on a limited-statistics sample of simulated events. 


Two specific examples of such improvements to the event classification are described 
in this section. The first is designed to mitigate ACD inefficiencies around the mounting 
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Fig. 25. — (a) Distribution of minimal distance of the intersection point of the extrapolation 
of the best track in the event from the edge of the corresponding ACD tile; (b) distribution 
of energy deposited in the ACD tile that intersects the extrapolation of the best track 
in the event for a subset of P7TRANSIENT events — with some P7S0URCE cuts applied; (c) 
same as (b) but for the P7S0URCE event class. The energy ranges are 30-100 GeV for (a) 
and 100-300 GeV for (b) and (c). Only data from high Galactic latitudes \b\ > 50° and 
longitudes 90° < l < 270°, which features a low y-ray-to-background ratio, are shown (all 
the background studies described in this section use ~ 24 months of data, from August 2008 
to end of July 2010). The relative normalizations of the two Monte Carlo components (CR 
and 7) are adjusted to fit the data. 


holes of the ACD tiles. Figure 25 (a) shows the distribution of the closest distance of the 
extrapolation of the best track in the event from the edge of an ACD tile. Only events that 
are classified as 7 rays in the P7S0URCE class are shown. A peak starting 39 mm from the edge 
is clearly visible, corresponding to the closest distance of many mounting holes in the ACD 
from the edge of a tile, where charged particles entering the LAT often leave very small or 
undetectable signals due to inefficiencies in the ACD response near these mounting holes. A 
matching peak is visible in the simulation, showing a good example of the detailed description 
of the LAT detector model entering the MC simulation of the LAT. This particular source 
of residual background has not been removed by the event classification scheme but can be 
easily eliminated. For P7CLEAN and P7UTLTRACLEAN classes we remove all events where the 
best track extrapolates to a range between 35 and 45 mm from the closest edge of an ACD 
tile that additionally produces a signal in the first TKR layer (see § 3.4.4). 


The second example demonstrates the removal of effects from poorly simulated interac- 
tions like the inelastic interactions of alpha particles and heavier nuclei. Figure [25] (b) shows 
the distribution of the energy deposited in the ACD tile closest to the extrapolation of the 
best track in the event onto the ACD plane. The events included pass the P7TRANSIENT 


selection and additionally some of the cuts used to define the P7S0URCE class (omitting cuts 
that are effective for reducing heavy nuclei contamination). A peak at a few MeV is clearly 
visible in both data and simulation corresponding to residual protons traversing the ACD. 
The second peak above 10 MeV, corresponding to residual helium, is almost completely miss- 
ing in the simulation. Figure 25 (c) shows the same distribution after all selection criteria 


for the P7S0URCE class have been applied. The residual helium peak has been removed and 
data and simulation show good agreement. 


4.3. Estimating Residual Background from Distributions of Control Quantities 


There are limitations to the agreement achievable between simulated and experimental 
data with methods such as those shown in the previous section. In particular, the primary 
and secondary CR fluxes, which are important inputs for our simulation, are uncertainly 
known. Furthermore, efficiencies of the trigger and on-board filter might be under- or over- 
estimated in the simulation. Therefore we adjust the normalization of the total residual CR 
background independently in 15 energy bins between 100 MeV and ~ 600 GeV to better 
describe the counts observed in the calibration data sample with a low y-ray-to-background 
ratio defined in § |3.6.4| The scaling is based on the events in the P7S0URCE class. For 


these events, the distributions of two event properties, the P all estimator (see § 3.3.8) and 
the transverse size of the particle shower in the CAL, have different shapes for 7 rays and 
CRs. The shapes of the distributions are sufficiently distinct for extracting the contribution 
of both components by fitting a superposition of simulated 7 rays and CRs to the on-orbit 
data. Although the difference between the shapes of the distributions for 7 rays and CR 
background decreases with increasing energy for the P a u estimator, it does increase for the 
transverse shower size. Therefore the fit is performed on P a n for energies E < 3 GeV, and on 
the transverse shower size for E > 3 GeV. The CR background correction factors obtained 
by these fits are then used to adjust the residual background predicted by the simulations. 


Figure 26 compares the distribution of the P a n estimator and the transverse shower 
size between simulated and experimentally observed events after the normalizations of the 
predicted CR background and of the 7-ray simulation have been adjusted. Each plot refers 
to a representative energy band in which the fit was performed for the respective variable. 
Figure [27] shows the adjustment factors obtained in the fit as a function of energy. An ad- 
justment factor of 1 corresponds to the CR background intensity predicted in the simulation. 
The adjustment factors vary between 0.7 and 1.6, depending on energy. 


We use the predictions of the residual CR background from the MC simulation for the 
P7S0URCE, P7CLEAN. and P7ULTRACLEAN event classes multiplied by the adjustment factors 
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in Figure [27] as our best estimate of the residual background. We use the largest adjustment 
factor (1.59) an indicator of the relative uncertainty of our determination of the residual 
background. This uncertainty is found to be ~ 35%, i.e., (1.59 — 1.00) / (1.59) = 0.37. 




Fig. 26. — (a) Distribution of the P a n estimator in the energy range 200-400 MeV for 7 
rays and residual CRs; (b) distribution of the transverse shower size in the energy range 32- 
64 GeV. The ranges are representative of the energy intervals in which the fit was performed 
for each of the variables. Only P7S0URCE events from high Galactic latitudes \b\ > 50° and 
longitudes 90° < |/| < 270°, which feature a low 7 -ray-to-background ratio (i.e., the “dirty” 
calibration sample), are shown. The normalizations of the CR background simulation and 
the 7 -ray simulation are adjusted to fit the data. 

Figure [28] summarizes our best estimate of the differential particle background rates in 
the three high-purity event classes for the energy range between 100 MeV and 600 GeV. 
Our background model is likely inaccurate below 100 MeV, and therefore the background 
contamination cannot be reliably determined by means of Monte Carlo simulations in that 
energy range. 


4.4. Estimating Irreducible Backgrounds 


The term “irreducible” was introduced in § |2.5.4| for CR-induced background with a well- 
reconstructed 7 ray inside the LAT as its only signature and is therefore indistinguishable 
from 7 rays of cosmic origin. As the two main classes of irreducible backgrounds we listed 
the CR positrons that annihilate outside the detector, and the CR protons that inelastically 
scatter in the passive material surrounding the ACD. This irreducible background is mostly 
limited to energies below a few GeV due to its production mechanisms. We do not consider 
the 7 -ray emission from the Earth’s atmosphere as irreducible, as only 7 rays with large 
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Reconstructed energy [MeV] 


Fig. 27 . — Correction factors for the CR background simulation as a function of energy, 
determined from a fit to the experimental data. The filled band shows the 35 % systematic 
uncertainty. Note that the energies shown here are based on reconstruction under the hy- 
pothesis that the event is a 7 ray and most high-energy protons deposit only a small fraction 
of their total energy in the LAT. 
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Fig. 28 . — Best estimates of differential rates of residual particle backgrounds for the 
P7S0URCE (a), P7CLEAN (b), and P7ULTRACLEAN (c) event classes. Individual contributions 
from primary CR protons, primary CR electrons and the secondaries from CR interactions 
are shown; the corresponding count rates for the extragalactic 7-ray background measured 
by Fermi (Abdo et al. 2010f) are also overlaid for comparison. 


errors in reconstructed directions enter the samples usually chosen for high-level analysis. 
The contamination can be reduced by both a stricter selection on event quality to reject 
badly reconstructed 7 rays and by larger exclusion regions around the Earth. 

An estimate of the fraction of irreducible background in the cleanest event class (i.e. , 
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P7ULTRACLEAN) is informative as it represents a lower limit to the achievable background 
rejection. Since the irreducible background cannot be separated from the data we can deduce 
its amount only based on available MC information. 


To determine the irreducible background from positron annihilations, one can compare 
the relative fractions of electrons and positrons surviving at different stages of the event 
selection. The secondary CR leptons between 100 MeV and 3 GeV passing the on-board 
filter are composed of f~ u ~ 0.28 of electrons and f^ { ~ 0.72 of positrons. Technically, 
the positron component is the sum of a reducible and an irreducible part; however, at this 
stage, the data set is still overwhelmingly dominated by reducible charged CRs, so that the 
irreducible contribution is effectively negligible: 

f obf ^obf,red -/obgirr ~ ./obf, red' (®) 

The secondary CR leptons passing the P7ULTRACLEAN selection are composed of f~ c « 0.10 
of electrons and /+. = /+ red + /+ irr ~ 0.90 of positrons. Since the reducible electron and 
positron components are indistinguishable in the LAT, they scale identically 


/, 


f obf, red f uc 


/obf/i 


obf J uc . 


uc,red 


/« 


obf 


therefore we have: 


f+ . = 1 - f~ - f 

J ucarr J uc J \ 


f obf 

f obf f u 


uc,red 


/, 


0.64. 


( 7 ) 


( 8 ) 


obf 


(i.e., ~ 64% of the secondary leptons in the P7ULTRACLEAN event class are irreducible back- 
ground events from positron annihilations). 

The amount of irreducible background below 1 GeV from inelastic scatters of protons 
can be estimated by evaluating the fraction of the residual simulated CR protons that does 
not enter the volume surrounded by the ACD. This is the case for about 95% of the simulated 
CR protons passing the P7ULTRACLEAN selection. 

For geometric reasons these scatters predominantly occur at the edges of the LAT (about 
75% of the residual CR protons, while the remaining 25% scatter in the spacecraft body). 


Figure 29 shows the positions of the projected intersections with the top ACD plane for 
simulated CR proton events surviving the P7ULTRACLEAN selection with reconstructed ener- 
gies below 1 GeV. An enhancement of tracks from the edges of the LAT is clearly visible, 
but smeared out due to the finite accuracy of the direction reconstruction. This feature 
in fact suggests the possibility to suppress this type of irreducible background by rejecting 
events from the edges of the LAT. However, due to the size of the PSF at low energies 
tracks intersecting large regions of the LAT surface would have to be vetoed, resulting in an 
unacceptable loss of effective area. Therefore, such a veto has not been implemented in any 
of the event classes. 
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Fig. 29. — Positions of the projected intersections with the top ACD plane for simulated 
CR proton events surviving the P7ULTRACLEAN selection. The dashed line marks the edges 
of the area covered by the ACD. 


4.5. Treatment of Particle Backgrounds 


After the event selection, even the purest class will contain residual background from 
misclassihed CRs. This contribution needs to be accounted for in the modeling used for spec- 
tral and spatial analysis of celestial y-ray sources. In particular, in the maximum likelihood 
model fitting framework of the ScienceTools, the misclassihed CRs must either be modeled 
separately or subsumed into one of the diffuse y-ray components. 

The simplest approach is to add a new source to the model to account for the residual 
background due to particle leakage. In most cases, since CR rates are related to geomagnetic 
(i.e., Earth) coordinates, for time intervals greater than a few months residual background 
events become approximately isotropically distributed in sky coordinates. 


The possibility of deriving an effective background template rests mainly on the as- 
sumption that the incidence angle dependence of A e g is the same for CRs and y rays. In 
many cases we also use a simplifying assumption that the CR contamination rates in front- 
and back-converting event scale with the relative acceptances. However, as we will discuss 
below and in 


4.6, that is not always the case. 


The heuristic nature of such a source is evident if one considers how it changes for 
different event selections. When we analyze a real y-ray source using different event class 
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selections we expect to find the same spectral distribution within the known systematics; on 
the other hand, this template depends on the amount of residual background and therefore 
on the CR rejection efficiency of the event class, so the templates we derive for the various 
event selections are dramatically different. Each event class requires a dedicated background 
template. 

Under the above assumptions any isotropic 7-ray component (e.g., the contribution of 
unresolved extragalactic sources) is not separable from the background leakage by means of a 
standard high-level source analysis (likelihood analysis) without additional knowledge (e.g., 
of spectral differences), so the two components are collected into a single term, simply called 
the isotropic templat To derive an isotropic template for a given event selection we perform 
a likelihood analysis of the high- latitude (|6| > 30) sky, including all resolved individual 
sources and a model of the Galactic interstellar emissioi^J fitting the spectrum of the isotropic 
component. It follows that the derived isotropic template depends on the assumed model 
for the Galactic interstellar emission, notably on the inverse-Compton component, which is 
smooth and far from negligible even at high Galactic latitudes, since the interstellar radiation 
field and CR electrons are broadly distributed about the Galactic plane. Therefore each 
Galactic interstellar emission model requires a different isotropic template. 


Between ~ 400 MeV and 10 GeV the assumptions mentioned above are rather good. 
In the 2FGL catalog analysis (Nolan et al. 2012) a single isotropic template was used and 
no significant systematics were observed above 400 MeV. Outside this energy range the rate 
of residual background events in the back section is appreciably greater than for the front 
section: the use of a single isotropic template (describing a weighted average of front and back 
background contamination) leads to a small hardening of measured spectra of point sources 
(Nolan et al. 2012). The effect is maximum for low-significance, soft sources: on average the 
spectral indices of power-law spectra are hardened by less than half of the typical uncertainty 
in the measured spectral index. It is preferable to derive separate isotropic templates for 
front and back and use them in a combined likelihood approach^ if front and back events 
are kept separate in the analysis, but the magnitude of this effect does not warrant such a 
complication for many analyses. 


To derive the true isotropic 7-ray component from the measured isotropic component 
it is necessary to separately estimate and subtract the amount of residual background con- 


7 E.g., isO-p7v6source.txt , available at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 

8 E.g., gal-2yearp7v6jv0.fits , available at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 

9 For example, isotropic-iem-front-v02.txt and isotropic-iem-back-v02.txt, 
http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html. 


available at 
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tamination (see § 4.3). 


The strategy to account for the CR background by means of an isotropic template, 
however, fails in the case of the 7-ray emission from the Earth limb. Residual background 
events in the FoV due to limb emission reconstructed in the tails of the PSF will produce 
a distinct pattern on the sky; its shape will depend on the pointing history and the time 
and energy ranges under consideration, and will be different for front- and back-converting 
events. 


A more stringent cut on 9 Z will reduce the contamination at the expense of exposure 
in certain regions of the sky. In particular, below 100 MeV exclusion regions to effectively 
eliminate the residual background become prohibitively large and significant Earth limb 
emission remains in the data sample for the commonly-used zenith angle limit of 100°. 


For the analysis leading to the 2FGL catalog, 7 rays with 9 Z > 100° were rejected. The 
remaining Earth limb emission was characterized by a template derived from the residual 
emission visible in the 50 MeV to 68 MeV energy band and which extended up to about 
400 MeVP^j (for details see Nolan et al. 2012). However, this template should not be used for 
periods much shorter than two years. 


Finally, residual background associated with mischaracterized “back-entering’' 


(see 


7 rays 

2.5.3), are another specific background that does not follow the 7 ray acceptance. The 
back-entering” 


probability to accept “back-entering” 7 rays into the P7S0URCE event selection is ~ 1000 
times smaller than for “front-entering” 7 rays and they are assigned directions roughly 180° 
away from the true directions. We consider the effect of this background in § 4.6 and find 


that treating it as part of the isotropic background does not introduce significant errors into 
analyses of point sources. 


4.6. Propagating Systematic Uncertainties to High Level Science Analysis 


As discussed in the previous section, residual CR background is treated as an isotropic 
fictitious 7-ray source in high-level science analysis of astrophysical 7-ray sources. As this 
approximation becomes less than perfect, significant systematic uncertainties can arise. 


In addition, a slight inconsistency between front and back A e g (see § 5.5) complicates the 
issue further, causing additional uncertainties when deriving separately the isotropic emission 
for the two selections. In general, we can quantify the resulting systematic uncertainties by 


10 Available at http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog. 
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comparing estimates obtained from the front-only selection with the full data set. 


The isotropic templates derived for the 2FGL catalog analysis for P7S0URCE_V6 events 
(and released via the FSSC) can be used for analyses spanning timescales of many months. 
On short timescales, especially less than the ~ 53.4 day precession period of the Fermi 
orbit, changes in the distribution of geomagnetic latitudes through which the LAT passes 
cause the residual background rates to be strongly dependent on the exact orbital history 
of the spacecraft and on the CR spectra at different geomagnetic locations. Analyses based 
on short time selections could do better either by using a dedicated estimate of the CR 
and isotropic backgrounds, e.g., by using a nearby control region, or take particular care to 
assess the impact of a possibly incorrect spectrum and spatial distribution of the background 
counts. Often this is done by allowing the isotropic component some freedom in the fitting 
procedure; see, for example, § |8.4 


As discussed in § |4.1[ § |4.3| and § |4.4| P7CLEAN and P7ULTRACLEAN event classes have much 


lower levels of background contamination than P7S0URCE. Accordingly, these samples can be 
used to study the dependence of any particular analysis on the level of particle background 
contamination in P7S0URCE analysis. 


Finally, we have studied the distribution of residual CR backgrounds in the P7S0URCE 
event sample b s (E,p) by comparing the observed counts distribution n s (E,p) in that sample 
with the predicted distribution n s (E,p), which we obtain by scaling the distribution of the 
P7ULTRACLEAN n u (E,p) sample by the ratio of the exposure calculated with the P7S0URCE_V6 
IRFs S s (E,p) to the exposure calculated with the P7ULTRACLEAN_V6 IRFs S u (E,p). Specifi- 
cally, 

n s (E,p) =n u {E,p) £ *^'P) 

£u{E,p) 

b s (E,p ) =n s (E,p) - n s (E,p ) (9) 


We studied the correlation between residual background and exposure as a function 
of energy. The detailed results are beyond the scope of this paper, but in general b s is 
not strictly proportional to £ s . This implies that the effective acceptance for residual CR 
backgrounds in the P7S0URCE event sample is not the same as for 7 rays. 

Although the spatial distribution of the residual CR-background could impact studies 
of large-scale diffuse emission, the variation across a typical ~ 20° ROI used when analyzing 
point sources is less that the variation in the exposure (2 to 5%, depending on the energy). 
Furthermore, for bright sources with sufficient statistics to make high-precision measure- 
ments, the correlation factor between the source parameters and the normalization of the 
isotropic component typically has a very small magnitude (< 0.03). Accordingly we neglect 
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the spatial variations in the CR-backgronnd contamination when performing point-source 
analyses. 
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5. EFFECTIVE AREA 


In order to correctly evaluate the spectra of astrophysical 7 -ray sources we need to know 
the effective collecting area of the LAT. In fact, A e g depends on the geometrical cross section 
of the LAT as well as the efficiency for converting and correctly identifying incident 7 rays. 
Because of the complexity of determining these we use high statistics MC simulations to 
evaluate A e g. We then quantify any discrepancies between simulations and flight data, and 
if needed, correct the MC-based A e g accordingly. 

As mentioned in § [lj we express the effective area as a function of the incident 7 -ray 
energy and direction in the LAT instrument frame. Therefore, the exposure (£) at a given 
energy for any point in the sky depends on the effective area and the observing profile (see 
Eqn. Q). 

In practice, the observing profile depends on the direction in the sky, and is accurately 
known. Therefore, the uncertainties on A e g(E ,9, (p) are the dominant source of instrument- 
related systematic error. Of course, we must also consider the uncertainties on our measure- 
ments of the 7 ray direction (i.e., the PSF) and the 7 ray energy (i.e., the energy dispersion). 
However, as we will show in the next three sections, in many cases the uncertainty of A e g is 
more important to the analyses than those of the PSF and the energy dispersion. 


In § 5.1 and § 5.2 we will describe how we generate tables of A e g as a function of energy 


and incidence angle (for the front and back sections of the LAT separately), and how we 
apply small corrections to those tables to account for variations of A e g with orbital position 
and azimuthal direction of the incoming 7 ray. Then in § |5.3| we will describe how we have 


validated the MC predictions of the 7 -ray selection efficiency for all the stages of the analysis 
using calibration samples within the flight data set, while in § 5.4 we will describe corrections 


to the A eS tables motivated by disagreement between the measured and predicted efficiency 
in one step of the selection process. Finally, in § 5.5, § 5.6, and § 5.7 we will evaluate 
the systematic uncertainties on A e g and show how we propagate these uncertainties into 
estimated systematic errors on measured astrophysical quantities such as fluxes and spectral 
indices. 


5.1. Effective Area Studies with Monte Carlo Simulations 


The starting point of the A e g evaluation is a dedicated allGamma sample (§ 2.5.3). 
Since the 7 rays are generated uniformly in log(E) and solid angle, the effective area in any 
of the bins in which the parameter space is partitioned can be expressed in terms of the total 
number of generated events N gen and the number of events passing the 7 -ray selection 
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criteria within the specific bin centered at E = Ei, 9 = 9; and 0 = 0*,: 


A e s(Ei, 9j, (j>k) = (6 nr 2 


N, 


gen 


2tt 


logic E max - log 10 E ] 


10 J -'min 


A Qj'kJ \log 10 E maX:i - log 10 E ] 


min,i 


( 10 ) 


where A is the solid angle subtended by the bin j, k in 9 and 0, E rnm and F max give 
the energy range of the allGamma sample and E min i and E maxi are the boundaries of the 
i th energy bin. (See § 2.5.3 for more details about the numerical factors and about the 
allGamma simulations in general). In practice, since the effective area is routinely averaged 

:e and r 

logm E n 


over 0 in scientific analysis, we factor out the 0 dependence and rewrite Eqn. ( 10 ) as 

logic E raaz 
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logio E } 
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logio En 


x R^Qjifa) (11) 


where R(E, 9, 0) is a small (of the order of 10%) correction factor whose average over 0 is 1 
by construction for any E and 9 (see § 5.2.3 for more details). Examples of effective area 
tables averaged over 0 are shown in Figure [30} 



io 2 io 3 io 4 io 5 

Energy [MeV] 



10 2 1 0 3 1 0 4 1 0 5 

Energy [MeV] 


Fig. 30. — Graphical representation of the effective area tables for the P7S0URCE_V6 class, 
front (a) and back (b) sections of the LAT. 


When describing the instrument performance, we more commonly show the effective area 
at normal incidence as a function of the energy or the angular dependence of the effective 


area for a given energy (usually 10 GeV), as shown in Figure 31 


The integral of the effective area over the solid angle, called the acceptance 


A{E) = / A eS (E,9,(f))dn = 


A eS (E, 9 , 0) sin 9 d9 dtp, 


( 12 ) 


'0 Jo 


is another widely used performance measure and is shown as a function of energy in Figure [32} 
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Fig. 31. — On-axis effective area as a function of the energy (a) and angular dependence (b) 
of the effective area at 10 GeV for the P7S0URCE class. 




Fig. 32. — Acceptance as a function of energy for the P7S0URCE class (a) and for the other 
standard y-ray classes (b). 


Formally, the FoV is defined, at any given energy, as the ratio between the acceptance 
and the on-axis effective area: 


F0V < £ > ^ AMt^O) <13) 

Figure [33] shows that the peak FoV of the LAT for the P7S0URCE event class is of the order 
of 2.7 sr between 1 and 10 GeV. At lower energies the FoV decreases with energy, as 7 
rays converting in the TKR at large angles pass through comparatively more material and 
therefore are less likely to trigger the instrument. A similar (smaller) effect is observed 
at very high energy, where, due to backsplash from the CAL, it becomes more difficult to 
reconstruct events at large angles. Finally we note that, for geometrical reasons (we require 
events in the standard classes to intersect the CAL) the FoV for the back section is typically 
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larger than that for the front section. 



Fig. 33. — LAT FoV as a function of the energy for the P7S0URCE class. Front- and back- 
converting events are shown separately. 


5.2. Corrections to the Effective Area Derived from Simulations 

In this section we describe three refinements we made to the effective area characteri- 
zation based on experience with flight data. In all three cases we have simulated the effects 
that we had previously ignored or averaged out, but that we discovered could significantly 
impact particular scientific analyses. 


5.2.1. Correction for Ghost Events 


As explained in more detail in § 2.1.4, after the start of LAT operations, it became 


apparent that ghost signals led to a significant decrease in effective area with respect to the 
pre-launch estimates, for which this effect was not considered. The overlay procedure used 
to account for this effect, first introduced in the P6_V3 set of IRFs, is described in detail 
in 


2.5.1 and its impact on the effective area is shown in Figure [34 


5.2.2. Live Time Dependence 

The effect of ghost signals is corrected on average as described in the previous section. 
A smaller correction is necessary to account for the detailed dependence of A e g- on the CR 


81 



Fig. 34. — Effective area at normal incidence for the P6_DIFFUSE class, in the pre-launch 
version (P6_V1JDIFFUSE, not corrected for ghost effects) and in the updated version includ- 


ing a modeling of ghost effects (P6_V3JDIFFUSE). Adapted from Rando & the Fermi LAT 


Collaboration (2009). 


rates. To account for this we need an estimator for the rate of CRs entering the LAT; the 
obvious one is the trigger rate, but technical issues make this choice impractical. A variable 
that can be easily obtained from the pointing history hies and which is linearly correlated 
with trigger rate is the live time fraction Fi , the fraction of the total observing time in which 
the LAT is triggerable and not busy reading out a previous event. The average value of F t 
is ~ 90% and varies between 82% and 92% over the Fermi orbit. 


We bin events from a sample of periodic triggers according to the corresponding live 
time fraction and for each bin we produce a dedicated allGamma simulation (§ 2.5.3); for 


each of the resulting overlay data sets we derive A eS . We have found the dependence on 
the incidence angle to be small and so we choose to neglect it when studying this effect. 
Furthermore, we have found that at a given energy A e g varies linearly with the live time 
fraction. We perform a linear interpolation in each energy bin in accord with 


A eS (E,Fi) = A eS (E ) • (c 0 (E)F l + Cl {E)) (14) 


separately for front and back events. 

In Figure [35] we plot Cq and C\ as a function of energy for front-converting events in the 
P7S0URCE event class. As shown by the solid lines we use a simple piecewise linear fit to 
describe the energy dependence and the fit parameters are stored in the A e ff tables. As we 
mentioned, the effective area derived from allGamma simulations that have overlaid periodic 
trigger events is effectively corrected for the average effect of this live time dependence, so 
we treat this additional correction (which can be positive or negative) as a modulation of 
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Fig. 35. — Energy dependence of the linear fit parameters Co (a) and C\ (b) for the A e ff 
live time dependence given in Eqn. (14) for front-converting events in the P7S0URCE: : FRONT 
event class. 


the tabulated A e g. Correction parameters are read in from the A eS hies of each set of IRFs, 
and used to correct the calculated exposures. 


The resulting corrections to the average A e s can reach —30% at 100 MeV, decreasing 
at higher energies to < 5% above ~ 10 GeV. The uncertainties in the corrections are much 
smaller; studies using flight data confirm the Monte Carlo-based predictions to better than 
2%. Note that over a 53.4 day orbital precession period this effect will tend to the overall 
average correction described in 


5.2.1 with less than 1% variation across the sky. 


Since the correction to the effective area is based on live time fraction, which 
is very strongly correlated with the CR intensity, it avoids any direct biases from 
long term changes in the CR intensity associated with the influence of solar ac- 
tivity on the geomagnetic field. However, the correction does not address the 
possibility that the CR population changes during the solar cycle in such a way 
as to change the effective area dependence on the live time fraction. Given the 
small change in the daily averaged LAT trigger accept rate observed in the mis- 
sion to date (< 5%), we neglect this effect. 


5.2.3. </> Dependence 

The tabulated values of A e g are averaged over the azimuthal angle of incidence and 
shown in Figure [36} Much of the azimuthal dependence of the effective area is geometrical, 
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due to the square shape of the LAT and the alignment of the gaps along the x and y axes. 

The RMS variation of the effective area as a function of 0 is typically of the 
order of 5% and exceeds 10% only at low energies (< 100 MeV) or far off-axis 
( 9 > 60°) where the effective area is small, and at very high energies (> 100 GeV) 
where the event rate is small. 



Fig. 36. — Total effective area at 10 GeV as a function of the incidence angle 9 and the 
azimuthal angle 0 for the P7S0URCE event class. The plot is shown in a zenith equal area 
projection with the LAT boresight at the center of the image; the concentric rings correspond 
to 0.2 increments in cos(0). 


In order to parametrize the azimuthal dependence of A e g we fold the azimuthal angle 
into the [0°,45°] interval and remap it into [0, 1] using the following transformation: 


£ = - 
7T 


7T 


mod — ) 


7 r 


(15) 


(the transformation maps 0° to 1, 45° to 0, 90° to 1 and so on). The allGamma events are 
binned in energy and 9 and in each bin a histogram of £ is fitted (see Fig. 37 for an example). 
Front- and back-converting events are treated separately. The fitting function is: 


m = 1 + 90 ^- ( 16 ) 

The absolute scale is not important: we normalize the correction to result in an average 
multiplicative factor of 1, so that the average A e g- is tabulated. The fitted parameters g 0 an d 
gi are stored in the A e g tables. 

We note in passing that the default for high-level analyses is to disregard the azimuthal 
variations when calculating exposures because they average out when intervals of months 
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Fig. 37. — Example of A e s azimuthal dependence fit. The plot refers to the bin centered 
at 7.5 GeV and 30° for the P7S0URCE class, front section — a similar fit is performed in each 
(. E,8 ) bin. The folded azimuthal angle £ is defined in Eqn. (15) and the fit function in 
Eqn. (16). Note that this plots shows A e s relative to A e g for £ = 0. 


or longer are considerecf^j Although the combined 9 and 0 dependence of the ob- 
serving profile averages out only on year-long timescales, the 8-fold symmetry 
of the LAT combined with the rotation of the x-axis to track the Sun results 
in effective averaging over 0 on short time scales. In fact, we have found that 
ignoring the 0 dependence of the effective area results in only a small variation 
of the exposure on 12-hour time scales (< 3% RMS at all energies). 


5.3. Step by Step Performance of Cuts and Filters 

Before describing the studies we performed to validate our event selections it is worth 
recalling that we require background rejection of ~ 10 6 while retaining high efficiency for 7 
rays. To achieve this we must select events based on many different criteria that are applied in 
several stages (see § [2]and § [3]) , which complicates the task of measuring the overall efficiency. 
For validation purposes we specifically examine the agreement of the selection efficiencies ( 77 ) 
between flight data and MG simulation. Therefore, the most relevant quantity is 

t-) ^/data 
it = , 

Vnic 



11 Information on how to include the (/) dependence of the effective area in exposure calcula- 
tions can be found at http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/binned_likelihood_tutorial.html 
and http: / / fermi.gsfc.nasa.gov / ssc / data/analysis / scitools/help / gtltcube.txt . 
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which is in general a function of direction, energy and conversion point of the incident 7 
ray. The approach described in this section can be used both to quantify the systematic 
uncertainties in A e g and possibly to correct the IRFs derived from the MC simulations, in 
cases of severe discrepancies with the flight data (e.g., see § 5.4). 


5.3.1. Background Subtraction 

To measure the efficiency of a cut on flight data, we can perform background subtraction 
before and after the cut and compare the excess in the signal region in the two cases. To 
the extent that the background subtraction is correct, the cut efficiency is simply the ratio 
of the number of background subtracted events before and after the cut we are testing: 

n S A - rn b 1 

V = , 

Tls, 0 ^^ 6,0 

where the subscripts s and b indicate the signal and background regions, the subscripts 0 
and 1 indicate the samples before and after the cut under test, and r is the ratio of the size 
of the background region to the size of the signal region. Note that the sample after the cut 
is a subset of the sample before the cut. Therefore, for a reasonable sample size (n s > 10), 
the statistical uncertainty of the efficiency is: 

r _ f [(w B|0 - n S| i) + r 2 {n b) 0 - n M )]?y 2 + (n B| 1 + r 2 n M )( 1 - ?y) 2 
V [ (n s ,o - rn bfi ) 2 

An example application of this technique using data from the Vela calibration data set is 
shown in Figure |38| 

In many cases our cuts have significant overlap — two cuts may reject many of the same 
events. I 11 such cases, measuring the efficiency of the second cut after the first cut has 
been applied would give a quite different result than measuring the efficiency of the second 
cut without the first cut. Therefore, whenever possible, we use the background subtraction 
technique described above to evaluate the efficiency of each step of our event selection both 
independently of the other steps as well after all the other cuts have been applied. 




5.3.2. Track-Finding and Fiducial Cut Efficiency 

Although the track-finding and fiducial cuts are applied midway through our event 
selection process we choose to discuss them first for two reasons: (i) we require an event 
direction to be able to perform a background subtraction on flight data and (ii) we express 
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Fig. 38. — Measurement of the efficiency of a cut using the background subtraction technique: 
(a) count spectra in the signal and background regions for the flight data set — note that the 
background counts have been scaled by 1/2 to account for the different phase ranges of the 
two regions; (b) excess in the signal region before and after the application of the test cut 
for the flight data sets — note that, for both (a) and (b), the corresponding plots for the MC 
data sets are not shown; (c) efficiency of the test cut rj for the flight and MC data sets; (d) 
ratio R of the efficiencies as measured with flight data to the MC predictions. In this case 
the baseline sample is the P7S0URCE event selection on the Vela calibration sample and the 
cut being tested is the P7CLEAN event selection. 


A e ff as a function of E and 6. Therefore, for performance studies we need to apply some 
minimum event quality and fiducial cuts before considering events for analysis. The standard 
cuts require at least one track found, with at least 5 MeV of energy in the CAL and that 
the track extrapolates through at least 4 radiation lengths in the CAL (see § 3.4.1). 


Most events that fail these cuts have either poorly reconstructed directions, poorly 
reconstructed energies, or both — which makes it difficult to study the performance as a 
function of energy and direction. Therefore, we choose to study the efficiency of the track- 
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finding and fiducial cuts by selecting events that almost fail these cuts in the P7S0URCE 
calibration samples. Specifically, we study the fraction of events that are very close to the 
cut thresholds and verify that the flight data agree with our MC simulations, see Figure |39| 



Fig. 39. — Validation of the track-finding and fiducial cuts. The data points show the fraction 
of P7S0URCE events in the Vela (< 10 GeV) and Earth limb (> 10 GeV) calibration samples 
that pass the fiducial cuts with some margin: having more than the absolute minimum 
number of hits on a track (black points), or having more than 20 MeV of energy deposition 
and crossing more than 5 radiation lengths in the CAL (gray points). The curves show the 
MC predictions for comparison. 


5.3.3. Trigger Conditions and Trigger Request Efficiency 


The LAT hardware trigger, trigger configuration and on-board filter are described in 
§ 3.1| For our purpose here we first consider the fractions of y-ray events that have one of the 
five physics trigger conditions (TKR, ROI, CAL_L0, CAL_HI, CNO) asserted. Furthermore, only 
two of them (TKR and CAL_HI) effectively serve to initiate a trigger request. Of the others, 
ROI does not exist without TKR, CAL_L0 alone does not open the trigger window and from the 
point of view of selecting y-ray events CNO is primarily a veto rather than a trigger. Since 
the CAL_HI requires at least 1 GeV in a CAL channel, and our fiducial cuts and quality cuts 
require at least one track, we are effectively using the TKR condition as the primary trigger 
for y rays up to such energies that the CAL_HI is very likely to be asserted. 


Using the DIAGNOSTIC filter events described in § |3.1 we can measure the fractions of 
all trigger requests that have individual primitives asserted. However, because of the high 
particle background rates this does not really probe the trigger stage of the y-ray selection 
process. On the other hand, for each of the five relevant trigger primitives we measure 


this efficiency as a function of energy for the allGamma sample as well as the P7TRANSIENT 
selection. For the P7TRANSIENT selection we also measure the fractions of events with each 
trigger primitive asserted and compare these to the MC predictions. These are shown in 
Figure 40 The only notable discrepancy is that the MC over-predicts the fraction of events 
having ROI asserted at very high energies; this is likely related to imperfect simulation of the 
backsplash from the CAL, and since the CAL_L0 and CAL_HI are typically asserted for these 
events, it does not affect the trigger readout decision (see Table |3|. 




Fig. 40. — (a) Fractions of all triggered events with various trigger conditions asserted as a 
function of the energy of the simulated event and (b) fraction of all P7TRANSIENT events that 
have the same conditions asserted as a function of the reconstructed energy of the event, for 
MC simulations (lines) and flight data (points). “Multiple TKR” refers to the events with 
more than one three-in-a-row combinations asserted. 

• The CAL4L0 trigger condition requires 100 MeV of energy deposited in any CAL channel. 
This condition starts to be asserted for 7 rays at ~ 300 MeV and reaches full efficiency for 
events contained within the CAL at ~ 1 GeV. 

• The CAL_HI trigger condition requires 1 GeV of energy deposited in any CAL channel. 
This condition starts to be asserted for 7 rays at ~ 3 GeV and reaches full efficiency at 
~ 15 GeV. 

• The CN0 condition requires any ACD tile to have very large signal, consistent with the 
passage of a heavy ion. The CN0 condition actually serves more to veto than to select an 
event for readout, but only becomes active at very high energies (> 100 GeV), where CAL_L0 
and CAL_HI are already active. 

• The ROI condition requires any ACD tile in a predefined region of interest associated with 
a TKR tower to have a signal above 0.45 MIP. The ROI condition actually serves more to 
veto than to select an event for readout and becomes active at a few GeV. 
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• Since the TKR condition serves as the primary trigger it is very difficult to measure the 
efficiency of the condition. However, we can estimate how well the MC simulates this effi- 
ciency by studying how well it models cases where the events almost fail the trigger request 
conditions. This is possible by calculating the fraction of events with exactly one combi- 
nation of three hit layers in a row in the TKR. Most high-energy 7 rays have many such 
combinations and would have triggered the LAT even if one hit had been lost. However, for 
back-converting low-energy events (< 100 MeV) the fraction of single combination events 
becomes significant. 

Since the fiducial cuts for our standard event selection require that a track is found, and 
that the track extrapolates into the CAL, the interplay between the trigger and the standard 
event selections is actually quite simple. 

• At very low energies, the LAT trigger starts to become efficient at ~ 10 MeV and follows 
the efficiency of the TKR condition, which becomes fully efficient by ~ 100 MeV. 

• At around 1 GeV the CAL_L0 condition becomes active. By design, this is considerably 
lower than the 10 GeV where the R0I starts to be asserted because of backsplash. 

• Above ~ 10 GeV the CAL_HI condition becomes active and we no longer rely on the TKR 
condition as the primary driver of the trigger. 

Furthermore, events that are rejected because R0I is asserted are extremely unlikely to 
pass standard event class selections. Taken together, this means that the only part of the 
LAT energy band where the trigger has strong influence on A e s is below ~ 100 MeV. 


5.3-4- On-Board Filter Efficiency 


Although the GAMMA filter has many different cuts (see § 3.1.2), most events that are 


rejected by the gamma filter would be either rejected by the fiducial cuts (see § 3.4.1) or 
by the P7TRANSIENT event selection. Accordingly, we choose to study the efficiency of the 
GAMMA filter as a whole, and only for those events which pass all the other selection criteria 
for the P7TRANSIENT class event sample. 


Since we downlink a small fraction of the events that fail the GAMMA filter (see § 3.1), 
we can check that this is indeed the case for this diagnostic sample. However, the large 
prescale factor (250) for the DIAGNOSTIC filter and high level of background rejection in the 
P7TRANSIENT selection severely limit our statistics for this study, and we can do little more 
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Fig. 41. — Fraction of events passing the P7TRANSIENT selection that also pass the GAMMA fil- 
ter for the DIAGNOSTIC filter flight data, compared to the MC predictions. Note that point 
at 42 GeV is well above the 20 GeV high-energy pass criteria where the gamma filter be- 
comes fully efficient, so the low measured value of the efficiency is likely due to a statistical 
fluctuation in the subtracted background. 


than confirm that the GAMMA filter is highly efficient for events in the P7TRANSIENT sample 
(see Fig. |4ll). 


5.3.5. P7TRANSIENT Class Efficiency 

Measuring the efficiency of the selection for the P7TRANSIENT event class is the most 
technically challenging part of the A e g validation for two reasons: (i) at the output of the 
GAMMA filter the background rates are still great enough to overwhelm almost all traces of 
y-ray signals and (ii) the event analysis directs events that are tagged as likely due to CRs 
away from the remainder of the y-ray selection criteria, which means that many of the CT 
analyses are never applied, and cannot be used to construct a cleaner sample on which we 
can measure the efficiency of any of these cuts independently of the other cuts. 

Figure |42|(a) shows the efficiency of each part of the P7TRANSIENT selection on events 
from the allGamma sample that have passed the gamma filter. As stated above, the high 
levels of CR background make it infeasible to use flight data to obtain stringent constraints 
on the efficiency of these selection criteria. Therefore, similarly to what we did for the 
track-finding and fiducial cuts, we also study the events that almost fail these cuts in the 
P7TRANSIENT calibration samples: these comparisons are shown in Figure |42|(b) . 

Because of the difficulty in validating the efficiency of the P7TRANSIENT event selection 
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criteria, we have chosen to use the consistency checks described in 
systematic uncertainty in our A e g representation. 


531 to estimate the 




Fig. 42. — (a) Fractions of events in the allGamma sample passing the gamma filter and 
track-finding and fiducial cuts that also pass each part of the P7TRANSIENT (§ 3.4.2) event 


selection, (b) Fractions of events in the P7TRANSIENT calibration samples that pass similar 
cuts with some margin for MC simulations (lines) and flight data (points). The numbers in 
the legends refer to the list of cuts in § 3.4.2| 


5. 3. 6. P7S0URCE, P7CLEAN and P7UL TRA CLEAN Class Efficiencies 

Although the P7TRANSIENT event class is dominated by residual background across the 
entire LAT energy range, the background levels are at least reduced to the point where the 
y-ray signals in the calibration samples described in § 3.6 are clearly detectable. This makes 


the validation of the effective area from this point on much easier. We can compare the 
efficiency of each cut, as measured on flight data, to the MC prediction with the method 
described in 


5.3.1 These comparisons are shown in Figure [43 


5.4. In-Flight Effective Area 


We observed that the efficiency for one part of our Pass 6 event selection was sys- 
tematically lower near 10 GeV for flight data than for the allGamma we used to evaluate 
the effective area, and we attempted to correct the effective area tables to provide more 
accurate flux measurements for y-ray sources. To be more specific, the offending cut is the 
Pass 6-equivalent of the P7S0URCE cut on the quality of direction reconstruction described 
3.4.3| (item [4] in the numbered list). It is important to note that for Pass 7 we chose 


m 
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Fig. 43. — Fractions of events in the P7TRANSIENT calibration samples that pass each part 
of the (a) P7S0URCE event selection (§ |3.4.3[ ) and (b) P7CLEAN (§ |3.4.4[ ) and P7ULTRACLEAN 


3.4.5) event selections for MC simulations (lines) and flight data (points). The numbers 


in the legends refer to the list of cuts in § 3.4.3 and § 3.4.4 


instead to make this particular cut less stringent to avoid the need to make such a correction, 
so that the in-flight corrections discussed here do not apply to the effective area tables for 
the Pass 7 standard y-ray classes. Furthermore, we traced the discrepancy to limitations 
in pre-launch calibration algorithm of CAL light asymmetry that resulted in degraded po- 
sition and direction resolution in the CAL: above ~ 1 GeV consistency between the TKR 
and CAL position and direction measurements is a strong indicator of accurate direction 
reconstruction. 


We measure the ratio of cut efficiency between flight data and MC as a function of 
energy and incidence angle (see Figure [44]) — as described in § 5.3.1 — and we use it to correct 
the MC based A e g. Because most of our calibration sources have limited statistics relative to 
the allGamma samples we are forced to use fewer bins when calculating the efficiency ratios. 
Furthermore, to avoid inducing sharp spectral features in measurements of y-ray sources we 
smooth the energy and angle dependence of the efficiency ratio. Specifically, we split the 
data into two cos^) bins: [0.2, 0.7] and [0.7, 1.0]. For values between the bin centers (i.e. , 
cos(0) € [0.45,0.85]) we perform a linear interpolation in cos(6*). For values outside that 
range we use the correction factor from the appropriate cos((9) bin. 

This procedure yields a correction map that we use to convert from the MC A e g to 
our best estimate of the true in-flight A e g for the P6_V11_DIFFUSE set of IRFs, as shown 
in Figure [45} Note that these IRFs should only be used together with the corresponding 
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Fig. 44. — Ratio of the P6.V11 .DIFFUSE to P6.V3 .DIFFUSE A e g for front-converting (a) and 
back-converting (b) events. The P6_V11 .DIFFUSE A e g tables have a correction factor relative 
to P6_V3_DIFFUSE that is based on the ratio of the efficiencies between flight and simulated 
data for the selection cut on the direction reconstruction quality (see item [4] in § 3.4.3). 
The underflow bins (white areas) had no A e g in P6.V3.DIFFUSE, As we did not interpolate 
the correction factor along the energy axis, fluctuations in the correction factor lead to the 
vertical bands visible in this figure. 


rescaled models for Galactic and isotropic diffuse emission 


provided by the FSSCQ 


5.5. Consistency Checks 


In this section we will describe several consistency checks we performed to estimate 
how well we understand A e g. Each of these tests consists of splitting a specific y-ray event 
selection into two subsets and comparing the fraction in each of the subsets from the flight 
data to the predictions from MG simulations. 


As a simple check of our method, we split the event sample into subsets of events 
converting on the +x and — x sides of the LAT. As noted in § |5.2|the </> dependence of A e g is 


strongest between directions toward the corners of the LAT relative to directions toward the 
sides of the LAT; accordingly, we split the data into events coming from the sides (£ < 0.5) 
or corners (£ > 0.5), based on Eqn. (15). Finally, we tested the 6 dependence of the A e g by 
splitting the data into on-axis (cos 9 > 0.7) and off-axis (cos 9 < 0.7) subsets. The results of 
these tests are shown in Figure 46 In each of these examples, we used the P7S0URCE event 


12 Specifically, glLiem_v02-P6-VlFDIFFUSE.fit and isotropic-iem-v02-P6-VlFDIFFUSE.txt, available at 

http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html. 
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Fig. 45. — Comparison of acceptance between P6_V3_DIFFUSE and 
The only difference between the IRFs is the application of the 
P6_V11 -DIFFUSE described in 8 |R4 


P6_V11 .DIFFUSE IRFs. 

corrections to A e g for 


sample as the starting point. 


We also compare the fluxes we measure with different event classes. By doing so we can 
check the accuracy of our measured efficiency loss for each of the selection cuts to go from 
one event class to the next. Technically we do this by asking what fraction of events in one 
event class also remain into an event class with tighter selection. Figure 46 also shows the 
results of comparing the P7CLEAN selection to the P7S0URCE selection. 


We also compared the fluxes we measure with front-converting events relative to back- 
converting events (Fig. [47|. Since the primary difference between these two parts of the LAT 
is in the distribution of conversion material, this test is especially sensitive to issues with 
our Geant4 simulation and probes our modeling of the trigger and track-finding efficiency. 


In each case we find that the fraction of events in each subset for the flight data are 
consistent with MC predictions to better than 15%. In fact, for most of the cases the 
agreement is far better than that, closer to the 2-3% level. The most significant discrepancies 
we see are between front-converting and back-converting events (Fig. 47). 


5.6. Uncertainties on the Effective Area 

5.6.1. Overall Uncertainty of the Effective Area 


From the consistency checks described in the previous section we arrive at a rough 
overall estimate of the uncertainty of A e g, which is shown in Figure 48. Note that this 
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Fig. 46. — Ratio of the fraction of events as seen in flight data relative to the prediction 
from MC simulations: (a) 7 rays converting on the +x side of the LAT; (b) 7 rays pointing 
in 45° ranges of azimuth toward the corners of the LAT (i.e., £ > 0.5) ; (c) 7 rays from near 
the LAT boresight, i.e., for which cosd > 0.7; (d) 7 rays passing the P7CLEAN selection. 

estimate is assigned simply to account for the largest observed inconsistency, namely the 
mismatch between the front-converting and back-converting events. Roughly speaking, for 
the P7S0URCEJV6 and P7CLEANJV6 event classes these uncertainties may be be quoted as 10% 
at 100 MeV, decreasing to 5% at 560 MeV and increasing to 10% at 10 GeV and above. It is 
important to note that these uncertainties are statements about overall uncertainty of A e g 
at various energies, and do not include any statement about what types of deviations we 
might expect within the stated uncertainty bands, nor about the point-to-point correlations 
in any systematic biases of A e g. Those questions are addressed in the next sections. 
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Fig. 47. — Ratio of the fraction of front-converting events as seen in flight data relative to 
the prediction from MC simulations: 


5.6.2. Point-to-Point Correlations of the Effective Area 


Since our selection criteria are generally scaled with energy (see § 3.3.1), we expect any 
biases of A e g to be highly correlated from one energy band to the next. This point is very 
important when estimating the size of potential instrumental spectral artifacts. 

We have studied the point-to-point correlation of the effective area through the consis- 
tency checks described in § 5.5 where it is evident that the deviations from unity are not 
independent for neighboring energy bins. In order to quantify this correlation we first scale 
the values r t of the data-to-Monte Carlo ratio R in each of the N energy bins (indexed by 
z), turning them into normalized deviations 


dj, 


(r* - m) 


byEii (c-m) 


( 20 ) 


(N~ 1 ) 

and m is the arithmetic mean of the rA and then we construct the metrics 


~ d i-nf 
(N - n) 


Tn = 


( 21 ) 


The quantity in Eqn. (21) is related to a reduced XN-n( 2a;); ^ is small for highly positively 
correlated deviations (in which case the differences between neighboring bins are generally 
small), while the expectation value is 2, for all values of n, for normal uncorrelated errors. 

Figure [49] shows this metric, for different values of n, calculated on the front vs. back 


consistency check shown in Figure 47 (Only the Vela calibration data set is considered 


here). The extremely small value of T\ indicates that neighboring logarithmic energy bins 
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Fig. 48. — Systematic uncertainty band on A e g as a function of energy. The interplay 
between the steeply decreasing effective area and the degrading energy resolution below 
100 MeV and the resulting impact on spectral analysis will be thoroughly discussed in § [TJ 


are highly positively correlated (as can be naively inferred from the plot in Fig. 47), while 
on the scale of half a decade in energy (n = 4) there is little evidence of a correlation. This 
implies that the systematic uncertainties on the effective area are not likely to introduce 
significant spectral features over scales much smaller than half a decade in energy (which is 
much larger than the LAT energy resolution). The results for all the consistency checks are 
summarized in Table [8] 




Fig. 49. — (a) Values of r n (n — 1,2 and 4) defined in Eqn. (21 ) for the front/back consistency 
check with the Vela data set shown in Figure [47} The parent distribution of r n for random 
normal uncorrelated deviations is shown in gray, while the black histogram represents the 
values of r n for 10,000 random permutations of the original data points. (Note that in both 
cases the distributions are independent of n .) (b) Cumulative probability distribution for 
random normal uncorrelated deviations. 
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Consistency check 

Tl 

P(t < Tl) 

+x vs. —x 

1.24 

4.9 x 10" 2 

Azimuthal dependence 

2.21 

0.69 

9 dependence 

1.05 

1.6 x 10~ 2 

P7CLEAN selection 

0.95 

8.4 x 10~ 3 

Front vs. back 

0.66 

5.0 x 10” 4 


Tabic 8: Summary of the point-to-point correlations for the consistency checks described 


m 


5.5. Most of them indicate a strong positive correlation of the systematic biases between 


adjacent bins on the effective area (the expectation value for random normal uncorrelated 


deviations is r n = 2 for all values of n) 


For many analyses, especially those that are limited by statistics, it is enough to con- 
sider the overall uncertainty and allow for worst case deviations within the stated uncertainty 
bands. However, doing so will result in very conservative systematic error estimates. We will 
discuss this in more detail in § |5.7| when we describe techniques to propagate the estimates 


of the uncertainty of A e g to uncertainties on quantities such as fluxes and spectral indices. 
Furthermore, we will come back to the issues of point-to-point correlations and the poten- 
tial induced spectral features when we discuss the uncertainties associated with the energy 
reconstruction in § [TJ 


5.6.3. Variability Induced by Errors in the Effective Area 

As a source moves across the LAT FoV and A e g changes with the viewing angle, any 
errors in the A e g parametrization as a function of 9 potentially could induce artificial vari- 
ability. We have searched for such induced variability with Vela. We split the data set into 
12 hour periods (indexed by i) and compared the number of 7 rays observed (rif) during each 
period with the number of 7 rays we predict (hi) based on the fraction of the total exposure 
for Vela that we integrated during that 12 hour period. 

On average, our Vela calibration sample contains 230 (176) on-peak (off-pulse) P7S0URCE 
class events in the 100 MeV-10 GeV energy band every 12 hours. Since the off-pulse region 
is twice the size of the on-peak region, background subtraction yields an average on-peak 
excess of n t = 142 7 rays with an average statistical uncertainty of a* = 16 7 rays in each 
time interval. 

The exposure calculation requires several inputs: 
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1. the spacecraft pointing and live time history, which are binned in 30 s intervals; 

2. the P7S0URCE_V6 A e g parametrization, which we use when deriving A e s(E, t), the effective 
area for 7 rays for Vela for each 30 s interval; 


3. the P7S0URCE_V6MC PSF (see § 6.1) parametrization, which allows us to calculate the 
energy dependent containment C(E,t, 15°) within the 15° ROI for each 30 s interval; 


4. a parametrization of the spectrum of Vela F(E) so that we may correctly integrate 
Aeg(E,t) and C(E,t, 15°) over the energy range. 


In order to minimize dependence on the modeled flux, we calculate the exposure in- 
dependently for 24 energy bins (which we index by j). The exposure in a single time and 
energy bin is 

£ij= f ' J ' A s (E, t)C(E , t, 15°) -J^— dtdE. (22) 

We can then express the expected number of 7 rays in each time and energy bin (fijj) as a 
fraction of the total number of 7 rays in that energy bin (rij), 


n i,j ~ 


£■ ■ 


n 


j ■ 


(23) 


Then we sum riij and ruj across energy bins to hnd n* and n t . 

We have performed this analysis, dividing the first 700 days of the Vela data sample 
into 1400 12-hour time intervals and using the phase-averaged flux model 

F(E) oc E~ r e~^o (24) 


with T = 1.38 and E 0 = 1.360 GeV 
variation, 

the normalized residuals, 


ru 


o et al. 2009c). Figure 50 shows the fractional 


(25) 


ru 


ru 


ru 


ru 


<Ji 


(26) 


and the distribution of the normalized residuals for each of the time-intervals. The normal- 
ized residuals are very nearly normally distributed. Furthermore, the Fourier transform of the 
time-series (Fig. 51 ) shows only a small peak corresponding to the orbital precession period 
and is otherwise consistent with Poisson noise. Note that unlike more complicated analyses 
that involve fitting the flux of a point source, this analysis is testing only the accuracy of 
the Agfj representation (and to a much lesser extent, the 15° containment of the PSF). We 
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attribute the peak in the Fourier spectrum to small, incidence angle-dependent 
errors in the effective area. As the orbit precesses, the range of incidence angles 
sampled, and hence the bias in the calculated exposure, varies slightly. 

Although we performed this analysis with 12-hour time intervals, as noted 
in § 1 2 . 3 1 , the LAT boresight follows very similar paths across the sky during 
successive 2-orbit periods. Therefore, the level of instrument-induced variability 
observed with 12-hour time intervals is likely to be indicative of the systematic 
uncertainties for variability analyses down to 3-hour time scales. 




Fig. 50. — Check for variability-induced by errors in the A e s tables P7S0URCE Vela calibra- 
tion sample and the P7S0URCE_V6MC IRFs. Panel (a) shows the fractional difference and 
normalized residuals between the observed counts in the 100 MeV to 10 GeV band and the 
prediction based on the fraction of the total exposure accumulated during each of 1400 12- 
hour time intervals. Panel (b) shows the normalized residuals, which are very well fit with a 
Gaussian with unit width and zero mean. 


Although the estimate used in § |5.6.1| that the systematic uncertainty of A e g- is less 


than the disagreement between the extreme cases is quite conservative for long time scale 
observations, it is somewhat less conservative for shorter observations. For example, in 
observations less than the Fermi orbital precession period of ~ 53.4 days, a particular 
region of the sky might be preferentially observed at incidence angles where the bias of A e s 
is particularly large, or during parts of the orbit in which Fermi is exposed to particularly 
high CR background rates and the correction described in § |5.2.2| leaves some residual bias 


in the calculated exposure. Finally, we have observed that ignoring the </> dependence of 
the effective area (§ 5.2.3) can induce artificial quarter-yearly periodicity in the fluxes from 


directions near extremely bright sources, in particular the Vela pulsar. 
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Fig. 51. — The discrete Fourier transform of the time series. The only peak visible above 
the noise floor corresponds to the 53-day orbital precession period. Note that the figure is 
normalized and the vertical scale is expressed in units of the statistical uncertainty. 

5.7. Propagating Uncertainties on the Effective Area to High Level Science 

Analysis 

As we hinted in the previous section, translating uncertainties on A e g into systematic 
errors on quantities such as fluxes and spectral indices depends on the particular analysis 
and requires assumptions about the variation of A e g within the uncertainty bands. 



5.7.1. Using Custom Made IRFs to Generate an Error Envelope 


A somewhat brute force approach to this problem is to generate IRFs that represent 
worst case scenarios for measuring specific quantities like fluxes or spectral parameters and 
use these bracketing IRFs to repeat the analysis and extract the variation in the measured 
quantities. Of course, the nature of the variations between the IRFs depends on the quantity 
in question. 


We address this in a generic way by scaling A e g by the product of the relative systematic 
uncertainty 

SA,„(E) 


((E) = 


(27) 


A M (E) 

and arbitrary bracketing functions B(E) taking values in the [—1, 1] interval. Specifically we 
define modified A e g as 


K«(E,0) = A, S (E,9) ■ (1 + e(E)B(E)) . 


(28) 
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The simplest bracketing functions, B(E ) = ±1 clearly minimize and maximize A e g 
within the uncertainty band. On the other hand, to maximize the effect on the spectral 
index in a power-law fit, we choose a functional form that changes sign at the pivot or 
decorrelation energy Eq (i.e., the energy at which the fitted differential flux and spectral 
index are uncorrelated): 


B{E) 


± tanh ( - log (E/E 0 ) 


(29) 


The parameter k controls the slope of the transition near E 0 \ in practice we use k = 0.13, 
which corresponds to smoothing over twice the LAT energy resolution of A E/E ~ 0.15. The 
bracketing IRFs used for effective area studies are listed in Table [9] Figure 52 shows the 


Name 

B(E) 

c -flux Jo 

+ 1 

C-flux-hi 

-1 

cJndexsoft 

+ tanh(| log(E/ Efl) 

cJndexJiard 

— tanh(| log(E/ E 0 )) 


Table 9: Bracketing A e g and corresponding energy-dependent scaling functions used to create 
them. 


bracketing functions for cJndexsoft and cJndexJiard and their effects on the on-axis A&. 




Fig. 52. — Bracketing IRFs for A e g designed to estimate possible systematic biases of 
the measurement of the spectral index, (a) energy dependence of the scaling parameter 
e(E)B(E)‘, (b) acceptance of the resulting scaled A e g for P7S0URCE. For this example we 
used a pivot energy E 0 = 1 GeV. 

We have studied two sources from the AGN sample to obtain estimates of the effects of 
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instrumental uncertainties on the measured fluxes and spectral parameters, namely: 

• PG 1553+113 — associated with 2FGL J1555. 7+1111, which has spectral index of 1.66 ± 0.02, 
making it one of the hardest bright AGN; 

• B2 1520+31 — associated with 2FGL J1522. 1+3144, which has a spectral index of 2.37 ± 0.02 
(when fit with a power-law), making it one of the softest bright AGN. 

In each case we used the ScienceTools (version v9r25p2) to perform a series of binned max- 
imum likelihood fits to a 20° x 20° region centered at the source position over the energy 
range 100 MeV-100 GeV. For each individual fit we followed the same procedure: 

1. used the y-ray and time interval selection criteria as for the Vela calibration sample (see 

§H; 

2. included all 2FGL sources within 20° in our likelihood model, with the same spectral 
parametrizations as were used in the 2FGL catalog; 

3. included models of the Galactic diffuse emission ( ring-2year.P76-vO.fits ) and the isotropic 
diffuse emission (isotrop-2year-P76source-vl.txt) rescaled by the inverse of the function used 
to rescale the effective area (so as to ensure the same distribution of expected counts from 
these diffuse sources); 

4. freed the spectral parameters for all point sources within 8° from the center of the region, 
as well as the overall normalizations of both diffuse components. 


Following what was done for the 2FGL catalog, we used a simple power-law 


dN 

~dE 


N 0 



-r 


to model the differential flux from PG 1553+113 and a “log-parabola” 


dN 

~dE 


N 0 



-(a+pin(E/E 0 )) 


(30) 


(31) 


to model that of B2 1520+31. (All of this will be also relevant for the tests with bracketing 
PSFs described in § 6.5.1 and with energy dispersion included in the likelihood fit described 


m 


7.4). 


Table [TO] and [FT] shows the fits results for PG 1553+113 and B2 1520+31, respectively, 
using these A eff bracketing functions, as well as the integral counts (F25) and energy (S25) 
fluxes between 100 MeV and 100 GeV. 
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Bracketing 7l e ff 

N 0 

[MeV -1 cm -2 s -1 ] 

r 

E 2 5 

[cm -2 s -1 ] 

S 25 

[MeV cm -2 s -1 ] 

Nominal 

2.54 x 10 -12 

1.68 

6.91 x 10 -8 

1.19 x 10 -4 

C-flux-hi 

2.75 x 10 -12 

1.67 

7.32 x 10 -8 

1.31 x 10 -4 

C-fluxJo 

2.37 x 10 -12 

1.69 

6.54 x 10 -8 

1.09 x 10 -4 

cJndex-hard 

2.57 x 10 -12 

1.64 

6.45 x 10 -8 

1.30 x 10 -4 

cJndexsoft 

2.53 x 10 -12 

1.73 

7.44 x 10 -8 

1.11 x 10" 4 


Table 10: Fit parameters and integral fluxes obtained using the A e g bracketing IRFs for 
PG 1553+113. Note that the quoted precision is roughly equivalent to the fit uncertainties 
and the pivot energy is E 0 = 2240 MeV for this source. 


Bracketing A e g 

No 

[MeV -1 cm -2 s -1 ] 

a 


F 25 

[cm -2 s -1 ] 

<5*25 

[MeV cm -2 s -1 ] 

Nominal 

5.23 x 10 -10 

2.24 

0.08 

4.09 x 10 -7 

1.33 x 10 -4 

C-flux-hi 

5.60 x 10 -10 

2.26 

0.07 

4.44 x 10 -7 

1.44 x 10 -4 

c_ flux Jo 

4.90 x 10 -10 

2.22 

0.08 

3.79 x 10 -7 

1.24 x 10 -4 

cJiidex-hard 

5.20 x 10 -10 

2.15 

0.10 

3.93 x 10 -7 

1.35 x 10 -4 

C-indexsoft 

5.27 x 10 -10 

2.33 

0.05 

4.29 x 10 -7 

1.32 x 10 -4 


Table 11: Fit parameters and integral fluxes obtained using the A e s bracketing IRFs for 
B2 1520+31. Note that the quoted precision is roughly equivalent to the fit uncertainties 
and the pivot energy is E 0 = 281 MeV for this source. 
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The ranges of the fit values indicate propagated uncertainties from the uncertainty in 
A e ff (Table [l2|. It is important to note that the systematic error estimates resulting from this 
technique represent conservative estimates within the instrumental uncertainties, rather than 
random variations. Furthermore, many of the bracketing IRFs are mutually exclusive, so 
when considering relative variations between 7-ray sources it is more appropriate to compare 
how the relative values change for each set of bracketing IRFs. 


Parameter 

B2 1520+31 
(soft) 

PG 1553+113 
(hard) 

SN 0 /N 0 

+7.2% 

-6.3% 

+8.0% -6.9% 

hr (ha) 

+0.09 

-0.09 

+0.05 -0.05 

5/3 

+0.02 

-0.02 

- 

8 F 25 /F 2 5 

+8.5% 

-7.2% 

+7.7% -6.6% 

ss 25 /s 25 

+8.1% 

-6.9% 

+10.0% -8.3% 


Table 12: Systematic variations arising from uncertainties in the effective area. For the 
spectral index (T or a) and spectral curvature (/3) we give the absolute variation with respect 
to the nominal value (e.g., hr). For the flux prefactor and the integral fluxes we give the 
relative variation with respect to the nominal value (e.g., 5N 0 /N 0 ). 


5.7.2. Using a Bootstrap Method to Generate an Error Envelope 


Alternatively, given a family of plausible A eff curves, we can use a weighted bootstrap 

for more details) for propagating the systematic 


approach (see Efron & Tibshirani 1993 


uncertainties on A e g. The weighed bootstrap approach is closely related to the bracketing 


IRFs method described in the previous section and to the methods discussed in Lee et al. 


(2011) in the context of the analysis of Chandra X-ray data. 


The basic idea is that, for each trial, the event data are bootstrap resampled using a 
weighting based on an effective area scaling function that is drawn from a family of plausible 
curves. The simplest of such families of A e ff curves (Fig. 53) can be constructed starting 


from Eqn. (28) and multiplying the scaling function e(E)B(E) by a normally distributed 
random number £ with zero mean — which effectively becomes the parameter controlling the 
family itself: 

A' eS (E, 9 , 0 = A eS (E, 9)-( 1 + £e(E)B(E)) (32) 


We have found the results of the weighted bootstrap with the family of scaling functions 
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Energy [MeV] 


Fig. 53. — A family of effective area scaling functions obtained for different values of the 
normal random variable £ in Eqn. (32). 


defined in Eqn. (32) to be in good agreement with those of the bracketing IRFs approach 
described in the previous section. This should not be surprising as the bracketing IRFs 
use the ±lcr excursions of the bracketing function and the weighted bootstrap draws from 
a Gaussian distribution of function scalings that are also based on that same bracketing 
function. The real benefits of the weighted bootstrap arise when one has families of plausible 
effective area functions that have more complicated dependencies (e.g., on incidence angle 
as well as energy) such that exposures would need to be recalculated to apply the bracketing 
method. 


5.8. Comparison of Pass 6 and Pass 7 


The Pass 7 event classes were designed to meet the same high-level analysis needs as 
the Pass 6 classes. Table [5] summarizes the correspondence between the Pass 7 standard 
y-ray classes and their closest Pass 6 equivalents. The corresponding acceptances are shown 
in Figure [54} 


As discussed in § [3} the main technical improvement in Pass 7 was optimizing the event 
selections on simulations that included the ghost signals, and on flight data. The principal 
outcome for Pass 7 is a substantial increase of A eS below ~ 300 MeV, especially for the 
cleaner event classes. At 100 MeV the P7S0URCE_V6 and P7CLEAN_V6 y-ray classes feature 
an acceptance of ~ 0.3 m 2 sr, to be compared with the ~ 0.1 m 2 sr of the Pass 6 equivalent 
(Fig. 54). Pass 7 has opened a window on astronomy with the LAT data below 100 MeV, 
though the reader should bear in mind the specific caveats in 
spectral analyses at these energies. 


6.5 and § 7.4 when performing 
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Fig. 54. — Comparison of the acceptances for the Pass 7 and the Pass 6 standard 7 -ray 
classes. 

Additionally, in Pass 7, the energy dependence of A e g is smoother. This was accom- 
plished in conjunction with obtaining a better overall understanding of the effective area 
itself, owing to the extensive use of on-orbit data to verify the fidelity of our MC simulation 
at each step of the 7 -ray selection process. 
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6. POINT-SPREAD FUNCTION 


As discussed in S 12.1.1 


at low energies the PSF is determined by multiple scattering. 
For example, the calculated multiple scattering for a normally incident 100 MeV 7 ray 
converting in the middle of a thin (front) section tungsten foil is ~ 3.1° (space angle). The 
68% containment angle as measured from flight data for 100 MeV 7 rays near the LAT 
boresight averaged over all towers is ~ 3.3° and is in agreement with the MC simulations 


see 


6.2 ). The small difference is due to missing measurements when the trajectories happen 


to pass through regions without SSD coverage and the fact that the electron and positron 
from a conversion can undergo hard scattering processes such as bremsstrahlung. 


If multiple scattering were the only consideration, the PSF should become narrower 
as E~ 1 . The measured PSF however improves more slowly with energy, instead falling as 

. This slower improvement relative to that expected for pure multiple 
scattering is also due to missed measurements and hard scattering processes and is predicted 
by the MC calculations. 


- 1 - 0.78 


see 


3.3.3 


Above a few GeV the narrowing of the PSF with energy is limited by the finite hit 
resolution of the SSDs. The strip pitch of 228 /mi and the lever arm for the direction 
measurement result in a limiting precision for the average conversion of ~ 0.1° at normal 
incidence. The transition to this measurement precision-dominated regime as predicted by 
the MC should occur between ~ 3 GeV and ~ 20 GeV. Estimates of the limiting PSF from 
flight data however indicate a worse performance above ~ 3 GeV; the PSF 68% containment 
levels off at almost double the calculated value (i.e. , ~ 0.16°). This departure is one of 
the few instances where the MC results significantly differ from real data (see § 6.2). The 


LAT collaboration has identified limitations in the pre-launch calibration algorithms of the 
CAL light asymmetry (for more details see Abdo et al. 2009f) as the primary cause of 
these discrepancies, and is assessing the improvement of the flight data-derived PSF at high 
energies for data which were reprocessed with improved calibration constant 


13 


see 


3.3.3 


for details of how the CAL energy centroid is used in the event direction analysis). 


6.1. Point-Spread Function from Monte Carlo Simulations 

Equivalently to what we stated in § [lj the PSF is the likelihood to reconstruct a 7 ray 
with a given angular deviation 5v = \v' — v |. We write it as P. 


1 ^ http: // fermi.gsfc . nasa.gov /ssc / data / analysis /documentation / Cicerone /Cicerone_L AT JRF s /IRFJPSF . html 
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As for the effective area, events from a dedicated allGamma MC simulation that pass 
the selections for the event class in question are binned in true energy and incidence angle. 
Note again that we ignore any 0 dependence of the PSF. As discussed below, in almost 
all cases the 0 dependence of the PSF is much weaker than the 6 dependence, 
and in § 8.4 we show that ignoring the 6 dependence induces at most a ~ 4% 
RMS variation of the flux. To allow for some additional smoothing of the variations 
of the parameters of the PSF fitting function with energy and angle a procedure based on 
running averages is implemented: in each bin events belonging to nearby bins (by default 
±1 bin, where the bin sizes are 0.1 in cos 6 and 0.25 decades in energy) along the energy 
and angle axes are included. For every bin we build a histogram with the angular deviations 
of detected 7 rays. This distribution is fitted and the fit parameters are saved. Note that, 
although the PSF is parametrized in the LAT reference frame, the angular deviation is the 
same whether expressed in celestial or LAT reference frames. 

Since the PSF varies with 0, it is often useful to consider the PSF averaged over the 
observing profile 


2.3) for a source of interest: 


P(Sv] E) = 


f P(dv;E,0)A eff (E,0)t ohs (0)dn 
J A eff (E,0)t ohs (0)dn 


(33) 


6.1.1. Point-Spread Function: Scaling and Fitting 


In our parametrized description of the PSF most of the energy dependence is factored 
into a scaling term: 


S P (E) = 


U 


Co • 


E 


V 

100 MeV J 


+ 4 - 


(34) 


Despite a careful investigation, we did not find a simple satisfactory description of the 6 
dependence to be incorporated in the scaling function. 

When building our MC-based PSF we use a set of scaling function parameters based 
on pre-launch simulations and confirmed with analysis from beam tests with the Calibration 
Unit (see § 7.3.1 for more details). The values of these parameters are shown in Table 13 


Note that S'p(E) has same functional form as for C 68 in § 3.3.3, however we have updated 


the parameters slightly based on the scaling observed in our allGamma sample. 
We define the scaled angular deviation x as: 

8v 

Sp(E) ' 


x = 


(35) 
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Conversion type 

c 0 

[°] 

/-> r c 

c i 

)' 

P 

Front 

3.32 

0.022 

0.80 

Back 

5.50 

0.074 

0.80 


Table 13: Parameters of the angular deviation scaling function Sp for the PSF parametriza- 
tion. 


An example of scaled deviation is shown in Figure [55} The effect 


the profile almost independent of energy, in that the maximum is 
all energy bins while the PSF 68% containment varies by almost 
from 100 MeV to 100 GeV. 


of the scaling is to make 
always close to x — 1 for 
two orders of magnitude 


Before the fit is performed, each scaled deviation histogram is converted into a proba- 
bility density with respect to solid angle. The result is illustrated in Figure |55j 




Fig. 55. — (a) Scaled angular deviation histogram and (b) PSF estimate in the range E = 
[5.6, 10] GeV and 6 = [26,37]° for the P7S0URCEJV6 event class, front section. 


Note that, although the scaling removes most of the energy dependence, the simulation 
indicates significant variation of the PSF with 6. At larger incidence angle the tracks must 
cross more material in each TKR plane. At energies below ~ 1 GeV this degrades the PSF 
owing to the increased multiple scattering, while at higher energies (above ~ 1 GeV) the 
additional complication of hard scattering processes in the TKR and additional hits in the 
TKR from the nascent electromagnetic shower complicate the track finding and degrade the 
PSF. Figure [56] shows how the scaled containment radii evolve with energy and incidence an- 
gle. To test the <f> dependence of the PSF we have also measured the containment 
radii independently for events with £ > 0.33 and £ < 0.33, where £ is the folded 
azimuthal angle defined by Eqn. (15). The 68% and 95% containment radii for 
the two £ ranges differ by < 5% for all energies and angles except at high energies 
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(> 10 GeV) and low incidence angles (cos 6 > 0.7) for back-converting events, and 
at even higher energies (> 100 GeV) and large incidence angles (cos 6 < 0.7) for 
front-converting events. Even then the maximum difference in the 68% contain- 
ment radius for events is only 10% at 10 GeV and 25% at 100 GeV. In summary, 
the variations of the containment radii with 0 are many times smaller than the 
corresponding variations with 9 for all but the highest energies (>~ 100 GeV). 



Fig. 56. — Scaled 68% containment radius (top) and ratio of 95% to 68% containment radii 
(bottom) as a function of ( E,9 ), for front- and back-converting events in the P7S0URCE_V6 
IRFs. 


The base function for the PSF is the same as used by the XMM-Newton mission (Kirsch 


et al. 

2004; 

Read et al. 

2011 


K(x,a, 7) = 


2'ko 2 


1 

1 - - 

7 


1 x 2 
27 a 2 


-7 


(36) 


which Kirsch et al. (2004) refer to as a King function ( King 1962), and is isomorphic to the 
well-studied Student’s t- distribution (Student 1908). Note that this function is defined so as 
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to satisfy the normalization condition: 

K(x, a, 7) 2nx dx — 1; (37) 

the extra 2ttx comes from the integral over the solid angle dVt = sin(x) dxdcj) ~ 2ttx dx. 
However, at very low energies the PSF widens to the point that the small angle approximation 
fails by more than a few percent and K (x) sin(x) must be normalized numerically. 

To allow for more accurate descriptions of the tails of the distributions, we use the sum 
of two King functions to represent the dependence of the PSF on scaled deviation for a given 
incidence angle and energy: 

P{P) ,/core A ("H *^"core; Tcore) T (1 ./core ) A ('p 0"taib Ttail ) • (38) 



The a and 7 values are stored in tables of PSF parameters as SCORE, STAIL, GCORE and GTAIL 
respectively. Because of the arbitrary normalization used in fitting the PSF function, f core 
must be extracted from the NTAIL table parameter, in conjunction with SCORE and STAIL: 


fo 


1 

1 + NTAIL • STAIL 2 /SC0RE 2 ' 


(39) 


The fitting function has been revised several times since the development of the first 
preliminary response functions. The version described here is the one currently being used 
and is different from e.g., that used for P6_V3 IRFs. A description of the fit functions used 
in the past is given in § |6.1.2 


6.1.2. Legacy Point-Spread Function Parametrization 

The first set of publicly released IRFs, P6_V3, used a slightly different PSF parametriza- 
tion. Specifically, it allowed for only one a parameter and fixed the relative normalization 
of the two King functions by constraining the two to contribute equally at Xb = 2y/5 a. So 
for P6_V3 we used 


P(P) ,/core A (x . (X, 7core) (1 fc ore) A (x . (7. 7tail); 


(40) 


with 


fo 


1 


1 + I< (x b , a, 7core)/A' (x b , a, 7 ta ii) ’ 


(41) 
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6.2. Point-Spread Function from On-Orbit Data 


During the first year of the mission we observed that for energies greater than a few 
GeV the distributions of 7 rays around isolated point sources were systematically wider than 
the expectations based on the PSF estimated from the MC simulations. We observed the 
same discrepancies for pulsars and blazars. In order to obtain a more accurate description 
of the core of the PSF for sources that are observed at a typical range of incidence angles, 
starting from the then-current Pass 6 IRFs, we derived the PSF directly from flight data, 
by means of a stacking analysis of selected point sources. The details of this analysis are 
described in Ackermann, M. et al. (2012, in preparation). Here we summarize the procedure, 
the associated uncertainties, and the impact on high-level source analysis. 


6.2.1. Angular Containment from Pulsars 


In the 1-10 GeV energy range bright pulsars are excellent sources for evaluating the 
in-flight PSF: not only they are among the brightest 7 ray sources, providing abundant 
statistics, the pulsed emission is very stable and the angular distributions of 7 rays around 
the true positions can be estimated readily by phase selecting 7 rays. The Vela pulsar is 
the brightest pulsar and an analysis with adequate statistics can be based on Vela alone. 
In the remainder of this section we describe the procedure for Vela only; the extension 


to an arbitrary number of pulsars is straightforward (see Ackermann, M. et al. 2012, in 


preparation). 


We use the P7S0URCE Vela calibration sample, which we divided into front- and back- 
converting subsamples, and bin the data in 4 energy bins per decade. We then calculate the 
pulsar phase for each 7 ray. As mentioned in § 3.6.1 we define [0.12, 0.17] U [0.52, 0.57] as the 


on” interval and [0.8, 1.0] as “off” (see Fig. 15 for the phase histogram). Next we calculate 


the containment angles: the position of Vela is known with a precision that greatly exceeds 
the LAT angular resolution and which can be assumed as the true source position. We 
make histograms of the angular deviations from Vela for the on-peak and off-pulse intervals 
and normalize them for the relative phase ranges. To estimate the PSF from flight data we 
measure the containment radii from the difference between the histograms. 


6.2.2. Angular Containment from Active Galactic Nuclei 


Above ~ 10 GeV spectral cutoffs of pulsars leave AGN as the only attractive sources for 
studying the PSF. In Ackermann, M. et al. (2012, in preparation) we address the potential 
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contribution from pair halos around AGN and conclude that we see no indication that this 
phenomenon is the explanation for the PSF being broader than predicted. Thus, we treat 
AGN as point sources. Many 7-ray sources are considered to be only “associated” with AGN, 
as opposed to “firmly identified” (see Nolan et al. 2012, for a discussion of the distinction), 


because of the limited angular resolution of the LAT. In the present analysis we consider 
only sources with high-confidence associations. As for pulsars, the positions of AGN are 
known with high precision from other wavelengths, and the angular distances from the true 
directions can be calculated. To accumulate enough statistics we stacked several sources and 
performed a joint analysis. We selected AGN from among the LAT sources with the highest 
significances above 10 GeV outside the Galactic plane. This energy limit is set by the source 
density: below a few GeV the LAT PSF is broad enough that nearby sources frequently 
overlap. 


A significant difference with respect to the pulsar analysis is the necessity of modeling 
the background in evaluating the distribution of angular deviations. We assume that after the 
stacking of the sky regions far from the Galactic plane the background count distribution can 
be assumed to be isotropic. At each energy the background is modeled as a flat distribution 
normalized by the amplitude in an annulus centered on the stacked data set. The inner 
radius of this annulus was chosen to be significantly larger than the region containing 7 rays 
from the stacked AGN sample. The uncertainty of the containment radius in each energy 
bin was set to the RMS of a large sample of MC realizations for the signal and background 
distributions. 


6.2.3. Point- Spread Function Fitting 


We have developed a procedure (described in detail in Ackermann, M. et al. 2012, in 


preparation), to fit our PSF model to the measured containment radii for different energy 


ranges. Given the statistical limitations we use a single King function, Eqn. (36). For 


the same reason we do not measure the dependence of the PSF on the incidence angle, 
i.e., we calculate an acceptance- weighted average over the incidence angle. We first fit the 


experimental 68% and 95% containment radii (f?68 and R95) with Eqn. (36). Then we extract 
a new scaling relation. And finally, we use the fitted (rather that the measured) 68% and 
95% containment radii to obtain a new set of PSF parameters for each energy bin. By using 
the fitted containment radii, this procedure smooths out the statistical fluctuations across 
the energy bins. 

The 68% and 95% angular containment radii for the flight-based P7S0URCE_V6 PSF are 
shown in Figure |o7| 
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Fig. 57. — 68% and 95% containment angles as a function of energy for the P7S0URCE_V6 
event class. 


6.3. Uncertainties of the Point-Spread Function 


The uncertainty of the derived PSF was estimated by comparing the 68% and 95% PSF 
containment radii from a set of calibration point sources with the corresponding containment 
radii derived from the P7S0URCE7V6 PSF (Ackermann, M. et al. 2012, in preparation). The 


and 95% containment radii measure of the accuracy representation of the PSF in the 
core and tail, respectively. 


The analysis was performed as a function of energy with 4 energy bins per decade. To 
determine the accuracy of the PSF fit as a function of incidence angle, subsamples were also 
studied in which 7 rays were additionally split into three bins of cos 6 ([0.2, 0.5), [0.5,0.75), 
and [0.75,1.0]). 


Figure 58 shows the 68% and 95% containment radii as a function of energy for front- 
and back-converting 7 rays averaged over incidence angle. The smooth lines show the model 
predictions for the MC (P7S0URCE_V6MC) and in-flight (P7S0URCE_V6) IRFs. At energies 
below 3 GeV the containment radii match the MC PSF with fractional residuals no larger 
than 10%. Above 3 GeV the MC PSF begins to systematically underestimate the 68% 
containment radius by as much as 50% for both front- and back-converting 7 rays. As shown 
in Figure [58] the P7S0URCE_V6 PSF reproduces the flattening of the energy dependence of 
the PSF containment at high energies. However owing to the limitations of using a single 
King function to parametrize the PSF, this model over-predicts the PSF tails as represented 
by the 95% containment radii. 


At large incidence angles (cos# G [0.2, 0.5]) the LAT PSF broadens by approximately 
a factor of 1.5. Due to the sky survey observing mode of Fermi and the decreased effec- 
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Fig. 58. — 68% and 95% containment radii for 7 rays averaged over all incidence angles 
(cos 9 G [0.2, 1.0]) as a function of energy for front (top) and back (bottom). Data points 
with error bars show the measured containment radii derived from the Vela and AGN PSF 
calibration data sets. Solid lines show the P7S0URCE_V6 model predictions in each energy 
bin; dashed lines show the predictions from the MC simulations (P7S0URCE_V6MC). Residual 
plots indicate the fractional deviation with respect to P7S0URCE_V6. 


tive area at large incidence angles the fit of the in-flight PSF is dominated by 7 rays at 
smaller incidence angles (cos 9 € [0.5, 1.0]). The variation of the PSF with incidence angle is 
most relevant for the analysis of transient phenomena in which the time-scale of interest is 
comparable to or shorter than the orbital period of Fermi , such as GRBs and short-period 
time-series analyses. As shown in Figure [59| the agreement of in-flight and MC PSF models 
with the data appears worse at large incidence angles (cos# G [0.2, 0.5]) although the limited 
statistics limit a rigorous comparison. The in-flight PSF model, which does not incorporate 
6 dependence, significantly under-predicts the width of the 68% containment radius for both 
front- and back-converting 7 rays. The effect of these discrepancies on high-level analysis 
will be considered in 5 16.5 
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Fig. 59. — 68% and 95% containment radii for off-axis 7 rays (cos 6 1 G [0.2, 0.5]) as a function 
of energy for front (top) and back (bottom). Data points with error bars show the measured 
containment radii derived from the Vela and AGN PSF calibration data sets. Solid lines show 
the P7S0URCE_V6 model predictions in each energy bin; dashed lines show the predictions 
from the MC simulations (P7S0URCE_V6MC). Residual plots indicate the fractional deviation 
with respect to P7S0URCE_V6. 


6.4. The “Fisheye” Effect 

For sources observed only in a narrow range of incidence angle, particularly near the 
edges of the LAT FoV, we must consider an additional complication: particles that scatter 
toward the LAT boresight are more likely to trigger the LAT and be reconstructed than 
particles that scatter away from the boresight. Furthermore, since the event selections do 
not require the separate reconstruction of both tracks from the 7 — y e + e~ conversion and 
the reconstruction code will estimate an event direction using a single track if no vertex is 
found, for some events we base our direction estimate only on the particle that scattered 
more toward the LAT boresight. This effect increases as the energy decreases, since multiple 
scattering causes larger deviations. 
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Over the course of a year the position of the LAT boresight relative to any given direction 
in the sky is fairly uniformly distributed in azimuth; however, for shorter periods several 
factors can lead to large non-uniformities (see § 2.3). Furthermore, by construction the PSF 


is averaged over azimuth. Therefore, for long (> 1 year) integration periods the treatment of 
the PSF described in the previous sections is perfectly adequate and the fisheye effect simply 
results in a broadening of the PSF which is correctly described by the IRFs. However, for 
shorter periods the fisheye effect can result in systematic biases in localization, which should 
be accounted for. To do so we consider the polar and azimuthal components of the angular 
separation between the true and reconstructed y-ray directions in the allGamma sample. 
Taking v as the true y-ray direction (in the LAT frame), v' as the reconstructed direction, 
and z as the LAT boresight we can define the local polar and azimuthal directions: 


Z X V 

I z x LI 


e = 


0 X V 
10 X v\ 


Then we calculate the component of the misreconstruction along each: 
<50 = sin -1 ^0 • (v' — v ) j 59 = — sin -1 (§ ■ (v' — v ] 


(42) 


(43) 


(the extra negative sign in the above equation is applied so that the fisheye effect represents 
a bias toward positive values of 56). ft is also worth noting that simply considering the 
distributions of differences between the true and reconstructed 6 is complicated by the fact 
that the amount of solid angle varies with 6. 


Figure 60 shows how the mean of 59 varies as a function of energy and incidence angle 
for the P7S0URCE event selection. Although the bias can be very large at high incidence 
angles and low energies, it is important to recall that (i) the PSF is quite wide in those 


cases, and (ii) there is relatively little acceptance in that region. Figure 61 shows the ratio 
of the mean to the RMS of 59 for the same energies and angles, and we see that except for 
the furthest off-axis events and the lowest energies, the fisheye effect is a small contributor 
to the overall width of the PSF when considering persistent y-ray sources. 

We should recall that the P7TRANSIENT event selection does not include as tight con- 
straints on the quality of the event reconstruction as the P7S0URCE selection. Accordingly, 
the fisheye effect is more pronounced for this event selection, as can be seen in Figures [62] and 


63 This is of particular importance for GRBs, for which almost all of the exposure might 
be at a single incidence angle. For soft GRBs in the LAT, the bias for the P7TRANSIENT y 
rays near 100 MeV that contribute the most to the localization can be up to 6° at 50° off 
axis. Furthermore, since the orientation of the LAT boresight relative to a GRB 
might not change significantly during the GRB outburst, these y rays will tend 
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Fig. 60. — Mean of 59 for MC simulated data as a function of energy for several incidence 
angles for the P7S0URCE event selection for front-converting (a) and back-converting (b) 
events. There are not enough statistics for front-converting events at 6 — 75.5° to extract 
reliable values. 




Fig. 61. — Ratio of the mean to the RMS of 59 for MC simulated data as a function of 
energy for several incidence angles for the P7S0URCE event selection for front-converting (a) 
and back-converting (b) events. There are not enough statistics for front-converting events 
at 9 = 75.5° to extract reliable values. 

to be biased in the same direction, causing an overall bias in the localization of 
the GRB. Therefore, when statistics permit, more robust localizations can be 
obtained by using events with energies > 200 MeV, or by using P7S0URCE rather 
than P7TRANSIENT class events. 







120 




Fig. 62. — Mean of 66 for MC simulated data as a function of energy for several incidence 
angles for the P7TRANSIENT event selection for front-converting (a) and back-converting (b) 
events. There are not enough statistics for front-converting events at 6 = 75.5° to extract 
reliable values. 




Fig. 63. — Ratio of the mean to the RMS of 69 for MC simulated data as a function of energy 
for several incidence angles for the P7TRANSIENT event selection for front-converting (a) and 
back-converting (b) events. There are not enough statistics for front-converting events at 
9 = 75.5° to extract reliable values. 


6.5. Propagating Uncertainties of the Point-Spread Function to High Level 

Science Analysis 


Uncertainties in the PSF parametrization lead to an imperfect source model in the high- 
level source analysis (likelihood fitting), and therefore to systematic uncertainties in source 
spectra 


6.5.1), localizations and measurements of source extensions (§ 6.5.2) 


An additional source of systematic uncertainty comes from neglecting the 9 dependence 
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of the PSF in the derivation of the parametrization from flight data. In particular, this can 
lead to artificial variability as the PSF changes with the varying observing profile. 


6.5.1. Using Custom IRFs to Generate an Error Envelope 


We can create custom PSFs in a similar manner to that used for the effective area (see 


5.7), though with the additional complication that we have to explore variations in the 


shape of the PSF as a function of energy and incidence angle. Because the dependence of 
the PSF on 7 is not intuitive we choose to express the bracketing functions in terms of the 
observable quantities Rq 8 and r = R 95 /R 68 rather than in terms of a and 7. Specifically, in 
each bin of energy and incidence angle, we can define the bracketing values R' e8 and r' in 
terms of R e8 and r: 


-^68 — -^ 68(1 + e 6s(B)B 68 (E)) 

r' = v(l + e r (E)B r (E)). (44) 


We can then solve for the King function parameters <+ and 7' which would correspond to 
these values. 


We re-analyzed both the B2 1520+31 and the PG 1553+113 ROIs with each of the 
PSF bracketing functions listed in Table [14] using the procedure described in § |5.7| Based 
on the quality of the fits described in 


6.3 in particular the residuals on the 68% and 95% 


containment radii, we have assigned: 


^68 (E) — 10 % 

e r (E) = 50%. (45) 

Tables [l5l and [l6l shows the fits results for PG 1553+113 and B2 1520+31 as well as the 
integral counts and energy fluxes between 100 MeV and 100 GeV obtained using these PSF 
bracketing functions. The ranges of the fit values indicate propagated uncertainties from the 
uncertainty in the PSF (Table [17]). 

The greater influence of the uncertainty of the PSF on the flux and spectral measure- 
ments for the softer source (B2 1520+31) comes about because at lower energies the wider 
PSF makes resolving sources more difficult and results in greater correlation with the Galac- 
tic and isotropic diffuse components in the likelihood fit, which varied by up to ±4% and 
±5% respectively. 
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Name 

B§s(E) 

B r (E ) 

c_scalehi_t_nom 

+1 

0 

cscalelo_t_nom 

-1 

0 

C-pivothLLnom 

tanh(| log(E/ E 0 )) 

0 

cjpivotlojijnom 

— tanh(| log(E/ E 0 )) 

0 

cjnom_t_scaleh% 

0 

+1 

C-nom^Lscalelo 

0 

-1 

C-nom_Lpivothi 

0 

tanh(| log (E/E 0 )) 

cxnomjtxpivotlo 

0 

— tanh(| log(£y E 0 )) 


Tabic 
As in 


14: Bracketing PSFs and the energy-dependent scaling functions used to create them. 
5.7 we use k = 0.13, which corresponds to smoothing over A E/E ~ 0.30. 


Bracketing PSF 

N 0 

[MeV -1 cm -2 s -1 ] 

r 

F 25 

[cm -2 s -1 ] 

S 25 

[MeV cm -2 s -1 ] 

Nominal 

2.54 x 

10 -12 

1.68 

6.91 x 10 -8 

1.19 x 10 -4 

c_nom_t_scalelo 

2.48 x 

10 -12 

1.70 

6.94 x 10 -8 

1.13 x 10 -4 

c_nom_t_scalehi 

2.54 x 

10 -12 

1.66 

6.68 x 10 -8 

1.22 x 10 -4 

C-nom-t-pivotlo 

2.46 x 

10 -12 

1.68 

6.67 x 10 -8 

1.15 x 10 -4 

C-nomxt-pivothi 

2.58 x 

10 -12 

1.68 

6.99 x 10 -8 

1.21 x 10 -4 

cscalelo-t_nom 

2.46 x 

10 -12 

1.67 

6.53 x 10 -8 

1.18 x 10 -4 

C-sealehijijnom 

2.63 x 

10 -12 

1.70 

7.32 x 10 -8 

1.21 x 10 -4 

C-pivotlo-tjnom 

2.60 x 

10 -12 

1.70 

7.30 x 10 -8 

1.18 x 10 -4 

C-pivothi-t-nom 

2.49 x 

10 -12 

1.66 

6.50 x 10 -8 

1.21 x 10 -4 


Table 15: Fit parameters and integral fluxes obtained using the PSF bracketing IRFs for 
PG 1553+113. Note that the quoted precision is roughly equivalent to the fit uncertainties 
and the pivot energy is E 0 = 2240 MeV for this source. 


6.5.2. Effects on Source Extension 

For sufficiently long exposures, the LAT can spatially resolve a number of q-ray sources. 
In the 2FGL catalog, 12 spatially extended LAT sources were identified (Nolan et al. 2012) 
using 24 months of LAT data; and several additional extended sources have been recently 
resolved from these data by Lande et al. (2012) using special techniques for modeling the 
spatial extension. Understanding the possible spatial extension of LAT sources is important 
for identifying multiwavelength counterparts, and using a source model with the correct 
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Bracketing PSF 

[MeV" 1 

N 0 

cm - 

2 s- r 

a 

I 

P 

[cm 

^25 

~ 2 s' 1 ] 

S25 

[MeV cm" 2 

s' 1 ] 

Nominal 

5.23 

X 

10- 

-10 

2.24 

0.08 

4.09 

X 

10- 

-7 

1.33 

X 

10- 

-4 

C-norri-Lscalelo 

5.34 

X 

io- 

-10 

2.31 

0.06 

4.28 

X 

10- 

-7 

1.33 

X 

10 - 

-4 

C-nomJscalehi 

5.02 

X 

io- 

-10 

2.19 

0.09 

3.85 

X 

io- 

-7 

1.30 

X 

10 - 

-4 

c_nom-t_pivotlo 

5.00 

X 

io- 

-10 

2.22 

0.09 

3.86 

X 

io- 

-7 

1.26 

X 

10 - 

-4 

C-noni-t-pivothi 

5.43 

X 

io- 

-10 

2.29 

0.06 

4.34 

X 

io- 

-7 

1.38 

X 

10 - 

-4 

cscalelo-t_nom 

4.87 

X 

io- 

-10 

2.21 

0.08 

3.77 

X 

io- 

-7 

1.26 

X 

10 - 

-4 

cscalehLt_nom 

5.57 

X 

io- 

-10 

2.27 

0.07 

4.40 

X 

io- 

-7 

1.40 

X 

10 - 

-4 

cjpivotlo_t_nom 

5.27 

X 

io- 

-10 

2.31 

0.06 

4.26 

X 

io- 

-7 

1.34 

X 

10 - 

-4 

cjpivothijtjnom 

5.16 

X 

io- 

-10 

2.17 

0.10 

3.91 

X 

io- 

-7 

1.32 

X 

10 - 

-4 


Table 16: Fit parameters and integral fluxes obtained using the PSF bracketing IRFs for 
B2 1520+31. Note that the quoted precision is roughly equivalent to the fit uncertainties 
and the pivot energy is Eq = 281 MeV for this source. 


Parameter 

B2 1520+31 
(soft) 

PG 1553+113 
(hard) 

5 N 0 / N 0 

+6.6% 

-6.8% 

+3.5% -3.4% 

8 T ( 8 a) 

+0.07 

-0.07 

+0.02 -0.02 

5/3 

+0.02 

-0.02 

- 

8 F 25 /F 25 

+7.7% 

-7.7% 

+6.0% -5.9% 

8 S 25 /S 25 

+5.4% 

-5.5% 

+2.4% -4.9% 


Table 17: Systematic variations arising from uncertainties in the PSF. For the spectral index 
(r or a) and spectral curvature (/3) we give the absolute variation with respect to the nominal 
value (e.g., <7T). For the flux prefactor and the integral fluxes we give the relative variations 
with respect to the nominal value (e.g., 5N 0 /N 0 ). 


spatial extent produces more accurate spectral fits and avoids biases in the model parameters. 


In addition to using a correct spatial model, the accuracy of the PSF can also affect 
the analyses of extended sources; and using an incorrect PSF will result in biases both in 
the fitted model parameters and in the significance of any spatial extension that is found. 
For example, flight data indicate that the MC-based PSF in the P7S0URCE_V6MC IRFs is too 
narrow at energies > 3 GeV (§ 6.2). Fitting the extension of SNR IC 443 (Abdo et ah 2010e) 


with a uniform disk for data from the first two years of observations and using the flight- 
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determined PSF in P7S0URCE_V6, we find a best fit disk radius of IC 443 of a = 0.35° ± 0.01° 
(Lande et al. 2012). By contrast, fitting these same data, but using the MC-based PSF 
in P7S0URCE_V6MC, we find a best fit radius of cr = 0.39°. This corresponds to a ~ 10% 
systematic bias in the measurement of the extension of IC 443. Since IC 443 is a fairly hard 
source, with photon index T = 2.2, the bias found for softer sources would be somewhat 
smaller. 


Use of an imperfect model of the PSF can also affect the assessment of the statistical 
significance of the extension of the source. We determine the significance of a measured 
source extension using the test statistic 


TS ext — 2 log(L ext /L pt ), 


(46) 


This is twice the difference in log-likelihood found when fitting the source assuming it is 
spatially extended versus assuming it is a point source. The use and validity of this formula 
is described in (Lande et al. 2012). For IC 443, using the P7S0URCE_V6 IRFs, we find 
TSext = 640, corresponding to a formal statistical significance of ~ 25cr. Performing the 
same analysis using P7S0URCE_V6MC, we obtain TS ex t = 1300, corresponding to ~ 37o\ 


6.5.3. Effects on Variability 


Neglecting the dependence of the PSF on incidence angle (§ 6.3) may introduce time- 


dependent biases of the estimated flux of the source, due to the different distribution of the 
source position with respect to the LAT boresight in each time interval. To estimate the 
potential size of such biases we consider how great an effect on source fluxes we might see if 
we were to naively use the nominal containment radius for aperture photometry analysis (i.e., 
an analysis where we simply count the 7 rays within a given aperture radius). Specifically, 
we use the MC based PSF to extract the containment radii averaged over the first two years 
of the mission R(E) in the direction of the Vela pulsar. Then we split the data into much 
shorter time intervals (indexed by i) and compute the fraction of 7 rays Ci(E ; vela) that fall 
within R(E; vela) for each of those shorter intervals: 


Ci(E , vela) = 


I J 0 R<E ' vehl) A e g(E, 6; vela)t 0 bs,i(6 l ; vela) sin ddddr 
f A eB (E,6)t ohS: i(0;vela) sin Odd 


(47) 


Here we consider 68% and 95% containment radii. Figure 64 shows how 6*68,1 [E\ vela) 
and 695^ (E] vela) vary for each 12-hour time interval over the first 700 days of routine science 
operations. If there were no 6 dependence to the PSF this figure would show lines at 0.68 and 
0.95. However we clearly see that ignoring the 6 dependence will cause slight mis-estimates 
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in the fraction of 7 rays falling within the mission-averaged i? 6 8 and R 95 . In fact, the width 
of the bands (i.e. , the spread of that mis-estimation) indicate the errors in flux we would 
expect if we were making a flux estimate based on aperture photometry with aperture cuts at 
R. It is worth pointing out explicitly that the C^^{E) band is significantly narrower, simply 
because most of the 7 rays are already contained within the Rg^(E) and the derivative of 
the PSF is smaller. 

Understanding the effect of the PSF dependence on incidence angle will have on a 
full likelihood analysis is much more complicated, as it depends on the other sources in 
the likelihood model. However, it is reasonable to take the width of the Cqs(E) bands as 
indicative of the magnitude of the effect when analyzing sources in complex regions where 
source proximity is a issue and the width of the C^{E) bands as indicative when analyzing 
isolated sources. 



Energy [MeV] Energy [MeV] 


Fig. 64. — Fraction of the PSF contained at the mission averaged 68% and 95% containment 
radii in the direction of the Vela pulsar for the P7S0URCE_V6MC IRFs as a function of energy for 
front-converting (a) and back-converting (b) 7 rays for each of 1400 12-hour time intervals. 


6.6. Comparison of Pass 6 and Pass 7 


As a consequence of the issues discussed in § |5.4[ Pass 7 event classes have less system- 
atic uncertainty for A e g at the expense of a slightly broader PSF across the entire parameter 
space. 

In addition Pass 7 standard classes (P7S0URCE and cleaner) feature a PSF derived from 
data above 1 GeV; see § |6.2| Below 1 GeV, where the PSF is derived from MC simulations, 


the PSF parameters are recalculated in order to ensure a smooth variation of the containment 
levels as a function of energy (§ |6.2.3). 
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Figure 65 compares the 68% containment angles for Pass 6 and Pass 7 classes recom- 
mended for routine analyses of 7 -ray point sources. While the difference between P6_V3_DIFFUSE 
and P7S0URCE_V6 is due mostly to the data-derived PSF in the latter, there is also a clear 
difference across all energies between the P7S0URCEJV6 PSF and in-flight P6_V11_DIFFUSE 
PSF. 
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Fig. 65. — Comparison of the 68% containment radius of the PSF for P7S0URCE_V6 with 
respect to P6_V3_DIFFUSE (a) and P6_V11_DIFFUSE (b). The P6_V3 .DIFFUSE is plotted for 
normal incidence; P6_V11_DIFFUSE and P7S0URCE.V6 do not include dependence on the in- 
cidence angle. 


While the difference in the containment radii we described is typically within the PSF 
uncertainties this comparison was performed with the incidence angle dependence ignored, 
so the tabulated PSF is an average over the FoV weighted by the average exposure over long 
timescales. This may be an issue for the analysis of bright sources over short timescales, 
when the observation angle distribution in the LAT reference system is significantly different 
from the average. Future Pass 7 IRFs may address this by reintroducing the 9 dependence. 


Current Pass 7 IRFs parametrize the PSF with a simplified version of Eqn. (38) for 
which the tail term is omitted as described in 5 16.2.31 The less accurate model of the 


PSF tails in the Pass 7 classes with on-orbit PSF (to date P7S0URCE.V6, P7CLEANJV6 and 
P7ULTRACLEAN_V6) is probably not an issue for most analyses, but caution is recommended 
when comparing the dispersion of 7 rays around the true position of bright sources with the 
expected distribution derived from the IRFs. 
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ENERGY DISPERSION 


The dispersion of measured energies around the true values (i.e. , the energy dispersion, 
or redistribution function) is generally asymmetric, with the most prominent tail being to- 
ward lower energies. This feature, characteristic of thin electromagnetic calorimeters (such 
as the LAT CAL, as opposed to full-containment calorimeters) makes the energy redistri- 
bution difficult to parametrize. As a matter of fact, since most source spectra are steeply 
falling with energy, low-energy tails in the energy dispersion are relatively harmless, while 
overestimating the event energy can potentially lead to overestimating the hardness of the 
spectrum. In the event selections we specifically make an effort to suppress the high-energy 
tails, rejecting events for which we might overcompensate with the energy corrections (see 


3.2.1 and § 3.3.2) 


While the width and shape of the energy redistribution are well understood, by default 
the energy dispersion is not taken into account in the standard likelihood analysis, primarily 
due to computational limitations. Furthermore, as we will see in § 7.4, the effect of neglecting 


the energy redistribution is usually small enough that it can be ignored. When that is not 
the case, it is easy to check the magnitude of the bias induced with dedicated simulations. 
If needed, the energy dispersion can be accounted for in the spectral analysis either by 
specifically enabling the functionality in the binned likelihood ( ScienceTools version 09-26- 
00 or later) or by means of unfolding techniques, as we briefly describe in 


7.4 


It is worth stressing that an in-flight validation of the energy response is much more 
difficult than the corresponding validations of the effective area and the PSF described in 
the previous section, as there is no known astrophysical calibration source that provides a 
spectral line at a well defined energy that would play the role that point sources have for 
the PSF (see, however, the discussion in § 7.3.3). The energy reconstruction validation will 


be described in 5 7.3 


7.1. Energy Dispersion and Parametrization from Monte Carlo Simulations 

As for A e g and PSF, in order to derive the energy dispersion events from an allGamma 
simulation belonging to a specified event class are binned in true energy E and incidence 
angle 6 and the distribution of measured energy E' is fitted (and the parameters tabulated). 
A complex functional form is necessary and attention is required to avoid fit instability and 
overly rapid variation of the fit parameters. 
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7.1.1. Scaling 


6.1.1) prompted a similar approach 
The scaling function Sd cur- 


The success of scaling the PSF distributions (see 
to the histogramming and fitting of the energy deviations, 
rently employed, derived by fitting the 68% containment of the measured fractional energy 
distribution (E 1 — E)/E across the entire E-9 plane, is 


Sd(E,9) = c 0 (log 10 Ef + ci(cos6*) 2 + c 2 log 10 if + c 3 cosd + c 4 log 10 E cos 6 + c 5 . (48) 


Front- and back-converting events are treated separately. The parameters used since the 
Pass 6 analysis are listed in Table 18 These are included in the IRF hies for Pass 6 and 
Pass 7 event classes. 



c 0 

Cl 

c 2 

c 3 

c 4 

C5 

Front 

0.0210 

0.0580 

-0.207 

-0.213 

0.042 

0.564 

Back 

0.0215 

0.0507 

-0.220 

-0.243 

0.065 

0.584 


Table 18: Numerical values of the coefficients defining the energy resolution scaling function 
Sd in Eqn. (48) for the Pass 6 and Pass 7 event classes. 


The scaling function Sd allows us to define a scaled energy deviation x for which much 
of the energy and angular dependence is already accounted for: 

(E' - E) 

X ~ S d (E, 9)E 



This scaled value is calculated for each simulated event and histogrammed, as shown in 
Figure 66 Scaling the distribution achieves two main effects: the width of the core of the 
distribution is almost constant for all energies and incidence angles, and the distribution 
thus expressed is significantly simpler to parametrize than the un-scaled version. 


7.1.2. Fitting the Scaled Variable 


The scaled energy deviation, as illustrated by the example in Figure [66j has a well 
defined core accompanied by elongated tails. We found we can effectively fit this structure 
with several piecewise functions of the form: 


f I 

X — Xq 

V 2 

a 


R(x, x 0 , a, 7) = N exp 


7 ' 


(50) 
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Fig. 66. — Histogram of the scaled energy deviation, as defined in Eqn. (49), fitted with the 


function D(x) in Eqn. (51). The plot refers to the ( E,6 ) bin centered at 7.5 GeV and 30 c 
for the front section. 


Indeed we model the energy dispersion D(x) within each ( E , 6) bin as four piecewise func- 
tions: 


D{x) = < 


N l R(x, xq, a L , 7r) 

if 

(x — Xq) < —X 

N l R(x, x 0 , 07, 72) 

if 

(x - x 0 ) e [~x, 0] 

N r R(x, xo, a r , 7 r ) 

if 

(x - x 0 ) G [0, x] 

N R R(x,x 0 ,a R , 7 r) 

if 

(x — Xo) > x 


(51) 


The values of the split point x and of the four exponents 7 of the energy dispersion parametriza- 
tion in Eqn. (51) are fixed as specified in Table 19 Moreover, the relative normalizations 
are set by requiring continuity at x = xq and \x — xo| = x and therefore the fit is effectively 
performed with a total of six free parameters, which are stored in the IRF FITS hies: the 
overall normalization N r = Ni (NORM), the centroid position Xo (BIAS), the two core scales 
cr r (RSI) and 07 (LSI) and the two tail scales a R (RS2) and 07 (LS2). 


X 

7 l 

7 1 

7 r 

7 R 

1.5 

0.6 

1.6 

1.6 

0.6 


Table 19: Numerical values of the split point x and of the four exponents 7 in Eqn. (51 ) that 
are fixed when fitting the scaled energy deviations. 


7.1.3. Energy Resolution 

The energy resolution is a figure of merit which is customarily used to summarize in 
a single number the information contained in the energy dispersion parametrization. As 
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illustrated in Figure 67 we define the energy resolution as the half width of the energy 
window containing 34% + 34% (i.e. , 68%) of the energy dispersion on both sides of its most 
probable value, divided by the most probable value itself. We note that this prescription 
gives slightly larger values of energy resolution than using the smallest 68% containment 
window. 



Fig. 67. — Energy dispersion at 10 GeV for front-converting P7_S0URCE events 30° off-axis. 
The most probable value of the distribution and the 68% containment window are indicated. 


Figures [68] and [69] shows the energy resolution for P7S0URCE events as a function of energy 
and incidence angle. As mentioned in § |2.1.2 the energy resolution has a broad minimum 
between ~ 1 and ~ 10 GeV, degrading at lower energies due to the energy deposited in the 
TKR and at higher energies due to the leakage of the electromagnetic shower out the sides 
and the back of the CAL. Conversely the energy resolution tends to improve as the incidence 
angle increases. This is especially true at high energy, where a longer path length in the CAL 
implies less shower leakage — though front-converting events more than 55° off-axis tend to 
exit the sides of the CAL and have worse energy resolution. 


7.1.% Correlation Between Energy Dispersion and PSF 


As stated in §[lj we factorize the IRFs, effectively assuming that the energy and direction 
measurements are uncorrelated. We use our allGamma MC sample to test this hypothesis. 
We locate each event within the cumulative PSF and measured energy distributions: 

not 

Pp — I P(a;E,v)da 

Jo 

f E ' 

P D = D{E'-,E,v)dE' 

Jo 


( 52 ) 
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Fig. 68. — Energy resolution (68% containment half-width, as decribed in the text) as a 
function of energy and incidence angle for front (a) and back (b) conversions. 




Fig. 69. — Energy resolution as a function of energy on-axis (a) and incidence angle at 
10 GeV (b) for the P7S0URCEJV6 event class. 


If the IRFs perfectly described the MC sample, both of these distributions should be flat. 
Furthermore, if there are no correlations between the energy measurement and the direction 
estimate, the two-dimensional combined distribution should also be flat. In practice this is 
nearly the case; the correlation coefficient between Pp and Pp is small ( | Cp, p \ < 0.1) across 
most of the energy and incidence angle range, as shown in Figure 70 (a). However, in many 
bins of (. E , 6), we do observe a small excess of counts with Pp ~ 1 and Pp ~ 0. Those values 
correspond to events that are in the tails of the PSF and have a energy estimate that are 
significantly lower than the true y-ray energy. To quantify the magnitude of this effect, we 
considered the fraction of highly anti-correlated events f a with Pp > 0.98 and Pp < 0.02 as 
a function of (E,9). In the absence of correlations, this fraction should be f a ~ 4 x 10 -4 . 
In fact, for highly off-axis 7 rays in some energy ranges this fraction approaches f a ~ 0.02; 
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however, averaged over the FoV, it is in the 0.001 to 0.0035 range for both front- and back- 


converting events in at all energies, as shown in Figure 70 (b) 


Although we can not reproduce this analysis with flight data, we have studied the 
correlations between the energy and direction quality estimators (Pe, see 
see 


3.3.2 and P„ 


3.3.3) and found similarly small effects. 


o 

- 0.02 

-0.04 

-0.06 

-0.08 

- 0.1 

- 0.12 

-0.14 

-0.16 

-0.18 

- 0.2 


-T-r-n-mj . l l ' ' ' '""I 

(a) -a- P7SOURCEV6 Front 

P7SOURCEV6 Back 

- . q. • ■ - Q ' 

■ & "43-' 

-! " E> - ,---43-- 1 


I I II I I I I I I I I ll I I till 


10" 


10 J 


10 


10 5 

Energy [MeV] 



Fig. 70. — Correlation between the energy dispersion and the PSF in the allGamma (§ 2.5.3) 
sample: (a) correlation coefficient Cp t p between Pp and Pp as a function of energy and (b) 
fraction of events with Pp > 0.98 and Pp, < 0.02 (i.e. , fraction of events which are in the tail 
of the PSF and also in the low side tail of the measured energy distribution). 


In summary, averaged over several orbital precession periods any biases caused by the 
correlation between the energy dispersion and the PSF are negligible compared to other 
systematic uncertainties we consider in the paper. However, it is certainly a potential con- 
tributor to instrument-induced variability (see § 8.4). 


7.2. Spectral Effects Observed with Simulations 


the maximum likelihood (LH) energy correction algorithm described 
y in Pass 6) is by construction a binned energy estimator. We have 
observed that it introduces spectral artifacts corresponding to the bins used in creating the 
likelihood parametrization. In addition, due to the fact that the correction is not reliable 
above a few hundred GeV, the method is specifically designed not to return values above 
300 GeV. As a consequence, it tends to concentrate events into a relatively narrow feature 
just below this energy. 


As noted in 5 3.3.2 


m 


3.2.1 (and used onl 


Both aspects are illustrated in Figure [7TJ where count spectra from standard allGamma 
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simulations (§ 2.5.3) are shown for both Pass 6 and Pass 7. While the overall shapes of the 
spectra reflect the LAT acceptance for the corresponding event classes (P6_V3_DIFFUSE and 
P7_S0URCE_V6), a clear sawtooth structure is visible in the Pass 6 count spectrum (with a 
typical width comparable with the LAT energy resolution and peak-to-trough amplitude of 
the order of ~ 5%), along with a prominent feature at 300 GeV, above which energies from 
the LH estimator cannot be selected 




Fig. 71. — Finely binned (50 bins per energy decade) count spectra for events passing 
the P6_V3_DIFFUSE (a) and the P7J30URCEJV6 (b) event selections; (a) also shows the LH 
energy estimates for the subset of events where they are available (note that the LH energy 
estimate is not always available, nor is it always selected when it is available). The final 
energy measurement in Pass 6 reflects the artificial sawtooth structures of the underlying 
LH energy estimator, which was removed in Pass 7. For both panels the vertical hashing 
corresponds to the bin boundaries of the LH energy estimator. 


The binning of the count spectra in Figure [Tl] (50 bins per energy decade) is deliberately 
much finer than the instrumental resolution — and therefore much finer than the binning one 
would use for a real spectral analysis of Fermi - LAT data. As a matter of fact, the features of 
the LH energy estimator have little or no effect in most practical situations. One noticeable 
exception is the search for spectral lines, such as might occur from the annihilation or decay 
of a massive particle into a two-body final state including 7 rays (see Abdo et ah 2010c, for 
a more detailed discussion). 

As already mentioned, in Pass 7 we chose to remove the LH algorithm from consider- 


ation in the energy assignment — see Figure 71 (b) 
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7.3. Uncertainties in the Energy Resolution and Scale 


In this section we will briefly review the systematic uncertainties on the absolute energy 
scale and the energy dispersion. The energy measurement is a complex process involving 
several different steps (see § 3.2.1 and § 3.3.2). Verification of the process is correspondingly 


complex, involving tests with sea-level CR muons, tests at accelerators, and analysis of flight 
data. 


7.3.1. The Calibration Unit Beam Test Campaign 


Since a direct calibration of the LAT with a particle beam was impractical for schedule 
and cost reasons, the LAT Collaboration assembled for this purpose a dedicated Calibra- 
tion Unit (CU) composed of two complete flight spare towers, one additional CAL module 
and five ACD tiles. An intensive beam test campaign was performed on the CU, between 
2006 July and November, with the primary goal of validating the MC simulation used to 
define the event reconstruction and the background rejection. The CU was exposed to beams 
of bremsstrahlung 7 rays, protons, electrons, positrons and pions — with energies ranging be- 
tween 50 MeV and 280 GeV (Table 20) — produced by the CERN Proton Synchrotron (PS) 


and Super Proton Synchrotron (SPS). In addition, the CU was irradiated with 1-1.5 GeV/n 
12 C and 131 Xe beams in the Gesellschaft fiir Schwerionenforschung (GSI) facility with the 
purpose of studying the instrument response to heavy ions. 


Particle 

Line 

Energy 

Triggers 

7 

PS 

0-2.5 GeV 

12 M 

7 (tagged) 

PS 

0.02-1.5 GeV 

4 M 

e~ 

PS 

1-5 GeV 

6.4 M 

e + 

PS 

1 GeV 

2.5 M 

7l~ 

PS 

5 GeV 

0.6 M 

V 

PS 

6-10 GeV 

0.6 M 

e~ 

SPS 

10-280 GeV 

17.8 M 

7l~ 

SPS 

20 GeV 

1.6 M 

P 

SPS 

20-100 GeV 

0.8 M 

12 c 

GSI 

1-1.5 GeV/n 

3 M 

131 Xe 

GSI 

1-1.5 GeV/n 

1.5 M 


Table 20: Summary of the 2006 Calibration Unit beam test campaign, adapted from Baldini 
et ah (2007). 
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A complete review of the results of the beam test campaign is beyond the scope of this 


paper. The reader is referred to 

Baldini et al. 

(2007 

) for a description of the experimental 

setup and a review of the main results and to 

Ackermann et al. 

(2010b 

) for additional related 


material. 

From the standpoint of the energy measurement one of the primary goals of the beam 
test campaign was to validate the leakage correction algorithms and, ultimately, the ability 
of our MC simulation to reliably predict the energy resolution at high energy. The tests 
confirmed that we understand the overall shower development and the energy resolution 


to better than 10% up the maximum available electron energy, i.e. , 280 GeV (Ackermann 


et al. 2010b). As we will see in the next section, this implies that the effect of systematic 
uncertainties on the energy resolution itself, when propagated to the high-level spectral 
analysis, is essentially irrelevant in the vast majority of practical cases. 

The most significant discrepancy was observed in the raw energy deposited in the 
calorimeter, as measured values were on average ~ 9% higher when compared with simula- 
tions (with smaller fluctuations around this value, slightly dependent on the particle energy 
and incidence angle). The origin of this discrepancy is unknown, and possibly could be at- 
tributed to residual environmental effects not properly accounted for, although an imperfect 
calibration of the CU cannot be excluded. In §|7.3.3|we will see that we have indications from 


flight data that we do understand the absolute energy scale to a precision better than 


7.3.2. Crystal Calibrations with Cosmic-Ray Data 


The calibration of the CAL crystals is the starting point of the energy reconstruction 
chain and underlies all the subsequent steps (a detailed description of LAT on-orbit calibra- 
tions is given in Abdo et al. 2009f). As explained in § 2.1.2, the large dynamic range of the 


LAT CAL (2 MeV-70 GeV per crystal) is achieved by means of four independent chains of 


electronics per crystal, with different amplifications (or ranges ), as summarized in Table 21 


In the most common readout mode (the so called single-range) the highest gain range that 
is not saturated is selected as the best estimate and converted into a digital signal that is 
eventually used in the event reconstruction. 


What is relevant for the discussion here is the on-orbit calibration of the crystal response, 
as determined by the crystal light yield and the linearity of the electronics. The lower energy 
scales are calibrated using primary protons that do not undergo nuclear interactions in the 
LAT, which are selected through a dedicated event analysis. For each crystal the most 
probable value of the energy deposition, corrected for the path length, is compared with the 
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Range 

Diode 

Gain 

Energy range 

MeV /Bin 

lex8 

Small 

High 

2 MeV-100 MeV 

0.033 

LEXf 

Small 

Low 

2 MeV-1 GeV 

0.30 

hex8 

Large 

ffigh 

30 MeV-7 GeV 

2.3 

HEXf 

Large 

Low 

30 MeV-70 GeV 

20 


Table 21: Readout ranges and energy conversion factors for the CAL crystals, adapted from 


Abdo et al. (2009f). 


MC prediction and the conversion factor between digital signal and MeV is computed. This 
procedure was first tested and validated with sea-level CR muons. The high-energy ranges 
are calibrated in the energy range of overlap between LEXl and hex8 using events collected 
by a special on-board filter for selecting heavy ions (see Fig. 72), read out in full four-range 
mode (i.e., with all four energy ranges read out at each log end). The non-linearity of the 
electronics is characterized across the entire energy range using a dedicated internal charge 
injection system and is corrected for in the energy measurement. 



Fig. 72. — Distribution of the path length-corrected crystal energy deposition for a sample 
of events used to inter-calibrate the low-energy and high-energy CAL energy ranges. The 
peaks corresponding to the most abundant Galactic CR species are clearly visible, though 
the composition differs somewhat from the Galactic CRs because of secondary production 
in the ACD and TKR. 

Though we originally intended to use heavy primary nuclei for an independent calibra- 
tion of the high-energy ranges, the uncertainty in the scintillation efficiency of heavy nuclei 
in the CsI(Tl) crystals relative to electromagnetic showers makes the heavy ion peaks unsuit- 
able for an independent cross-check of the absolute scale. The magnitude of the effect was 
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, but we are un- 
certain how to scale the effect to the typical on-orbit energies (~ 5 GeV/nucleon or greater) 
with our desired accuracy of a few percent. 

We have used protons and Galactic CRs from Be to Fe at incidence angles ranging from 
on-axis to 60° off axis to demonstrate that the crystal scintillation signal is proportional to 
path length for each species — i.e. that the crystal response is piecewise linear over factors of 
two in signal. Because the peaks from Be to O are closely spaced and easily resolved, this 
method demonstrates that the CAL energy scale is linear over at least the range ~ 180 MeV 
to ~ 1500 MeV per crystal, so any residual error in absolute energy scale applies equally 
over that entire range. Unfortunately it is difficult to bridge the gap between protons at 60° 
(~ 22 MeV) and He on axis (~ 45 MeV), and between He at 60° (~ 90 MeV) and Be on 
axis (~ 180 MeV), so we cannot demonstrate linearity in this region with this technique. We 
are exploring alternative CR event selections and geometries to cover these energy regions. 
Nonetheless, the relative variations of the peak positions (for both non interacting protons 
and heavier nuclei) can be measured with high accuracy, which effectively allows us to 
monitor the time drift of the absolute energy scale and prompt the update of new calibration 
constants when necessary. The results in Figure [73] are in good agreement with the pre- 
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Fig. 73. — Relative variation of the absolute energy scale, as measured from the position of 
the proton peak position, throughout the first three years of the mission. 

launch estimates of the crystal light yield attenuation due to radiation damage, which results 
predominantly from trapped charged particles in the SAA. A fit with the function 

(t-t 0 ) 

s(t) = (po - pi) + pie^~ (53) 

yields a time constant r of the order of ~ 2 years, and predicts an overall shift (~ p 0 — pi) 
of the energy scale (which can be corrected for) on the order of 4.5% after 10 years. 
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measured with 600-1600 MeV/nucleon C and Ni beams (Lott et al. 2006 
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7.3.3. Absolute Measurement of the Energy Scale Using the Earth’s Geomagnetic Cutoff 


As mentioned at the beginning of this section, there are essentially no astronomical 
sources with spectral features that are sharp enough and whose absolute energies are known 
to such a level of accuracy that they can be effectively exploited for an on-orbit validation 
of the absolute energy scale. 


One exception, perhaps unique in the energy range of the LAT, is the narrow 67.5 MeV 
pion decay line predicted to originate from interactions of primary CRs with the surface of 
the Moon (Moskalenko & Porter:2007). Unfortunately this feature is located at the lower end 
of the LAT energy range. While it is not inconceivable that the LAT possibly could provide 
the first evidence for the existence of such a line, the limited energy and angular resolution 
at these energies, together with the brightness of the limb of the Moon in continuum 7 rays, 
make the lunar pion line impractical as an absolute energy calibrator. 



Fig. 74. — Examples of simulated CR electron count spectra in bins of Mcllwain L for the 
Fermi orbit. They are averaged over the FoV and folded with the energy resolution. Effec- 
tively these are the templates that the electron count spectra from real data are compared 
with in order to measure the absolute energy scale. The peaked shape is the result of the 
power-law spectrum of primary electrons being convolved with the screening effect of the 
Earth’s magnetic field, which, in each bin of Mcllwain L, is effective below a certain cutoff 
energy. The relative normalizations reflect the fact that Fermi spends the most live time at 
low Mcllwain L. 

A practical alternative is the geomagnetic rigidity cutoff in the CR electron spectrum. 
The LAT is effectively shielded from CRs by the Earth’s magnetic field at a low rigidity. 
The convolution of this shielding with the primary power-law spectrum results in a peaked 
spectral feature, whose shape and absolute energy can be predicted with good precision by 
taking advantage of accurate models of the geomagnetic field that are available (jFinlay et ah 
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2010) and particle tracing computer programs (Smart & Shea] 2005). See § 4.1 for additional 
details. 


Figure 74 shows an example of such predictions in different parts of the Fermi orbit. 
The count spectra, averaged over the LAT FoV and folded with the LAT energy resolution, 
are shown in linear, rather than logarithmic, scale in order to emphasize the peaked spectral 
shape. Comparison of the predicted and measured peak positions is an in-flight measurement 
of the systematic bias in the energy scale. Since the magnitude of the Earth’s magnetic field 
varies across the Fermi orbit, this approach has the potential to provide a series of calibration 
points, specifically between ~ 6 GeV and ~ 13 GeV. The fact that this is the energy 
range for which the energy resolution of the LAT is best is beneficial for the measurement 
itself. Moreover, since both electrons and 7 rays generate electromagnetic showers in the 
detector, they are effectively indistinguishable from the calorimeter standpoint — so that 
energy measurements for one species directly apply to the other. 

We have used this approach for an in-flight calibration of the absolute energy scale 
using one year of data. The details of the analysis are beyond the scope of this paper and 


are discussed in Ackermann et al. (2012c). The main conclusion is that the measured cutoff 
energies exceed the predictions by 2.6% ± 0.5% (stat) ± 2.5% (sys) in the range 6-13 GeV. 


7. 3. 4- Summary of the Uncertainties in the Energy Resolution and Scale 


Here we summarize the results related to the systematic uncertainties of energy mea- 
surements with the LAT; in the next section we will discuss how these impact the high-level 
spectral analysis. 


The results from the CU beam test campaign indicated that the energy resolution, as 
predicted by our MC simulations, is accurate to within 10% up to the maximum accessible 
beam energy (i.e., 280 GeV). As we will see in the next section, this implies that this 
particular systematic effect is negligible in most practical situations. 


Though we were not able to identify the reason for the - 
deposited in the CU calorimeter measured at the beam tests 


9% discrepancy in the energy 


7.3.1), we do have compelling 


7.3.3) that in 


indications from the measurement of the CR electron geomagnetic cutoff 
fact the systematic error on the absolute energy scale is smaller than that. We stress that 
this measurement involves the real LAT, with the real flight calibrations, in its on-orbit 
environment. Although the measurement derived from the CR electron cutoff applies only to 
a small portion of the LAT energy range, the ~ 9% effect measured at the beam test affected 
essentially the entire range of energy and incidence angle, with a weak dependence on both 
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see also the remarks in 8 7.3.2 about the additional evidence for the fractional error on the 


energy scale being constant over a wide energy range). This evidence is further supported 
by the excellent internal consistency of the analysis measuring the electron and positron 
spectra separately (Ackermann et al. 2012d). However, it is more difficult to constrain the 
energy scale at the low and high ends of the LAT energy range, where the energy resolution 
degrades by a factor of ~ 2. Based on the full body of information currently available we 
conclude that that the energy scale for the LAT is correct to +20/ — 50% of the energy 
resolution of the LAT at a given energy. This corresponds to an uncertainty of +2/ — 5% 
on energy scale over the range 1-100 GeV, and increases to +4/ — 10% below 100 MeV and 
above 300 GeV. 


Finally, the measured energies of the Galactic CR peaks are being monitored to gauge 
the time stability of the absolute energy scale, which we can control at the 1% level by 
applying calibration constants on a channel-by-channel basis. 


7.4. Propagating Systematic Uncertainties of the Energy Resolution and Scale 

to High Level Science Analyses 


Uncertainties in the energy resolution and scale introduce systematics in spectral anal- 
yses, the magnitude of which depends on energy and the spectral shape of the sources under 
study. Since these systematic effects are a consequence of event redistribution between energy 
bins, effects are also strongly coupled to the energy dependence of A eS . As the implications 
of the systematic uncertainties are so dependent on the analysis, the main purpose of this 
section is to illustrate the basic ideas and the tone of the discussion is deliberately general. 


The effect of the finite energy resolution is generally negligible for spectral analysis 
above ~ 100 MeV, which as already mentioned is one of the main reasons why it is not 
taken into account in the standard likelihood analysis. Even above ~ 100 GeV, where the 
energy resolution is significantly degraded by the leakage of the shower from the CAL, the 
effect is negligible compared with other sources of systematics — at least for sources with steep 
spectra, like primary CR electrons (Ackermann et al. 2010b). MG studies with simulated 
point sources, though, show that the finite resolution can induce a bias in the count spectra 
as high as 20-30% at 50 MeV for power-law spectra with T ~ 1.5. The effect is strongly 
dependent on the spectral index and is less severe for softer sources, as shown in Figure [75} 


A bias in a count spectrum, in general, does not trivially translate into an effect of the 
same order of magnitude in the parameter values derived from a spectral analysis. When 
fitting a source spectrum with a single power law, for instance, the combination of the 
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Fig. 75. — Bias factor for count spectra induced by not taking into account the finite energy 
resolution of the LAT. The plots refer to simulated point sources (with no background) with 
power-law spectra, for several different values of the spectral index. 


long lever arm of the high-energy data points and the inability of the spectral model to 
accommodate any curvature in the count spectrum, results in the energy dispersion having 
very little effect on the fit parameters (even though it can produce very significant residuals). 
Figure 76 (a) shows that in this setup the bias in the measured spectral index, when fitting 
down to 30 MeV, is smaller than 0.03 for any reasonable input spectrum. 
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Fig. 76. — Bias in the spectral index returned by a binned likelihood analysis with a minimum 
energy of 30 MeV (this setup is taken as a worst case scenario; increasing the minimum energy 
for the fit to 100 MeV decreases the bias by a factor of 2 or more). The plots refer to simulated 
point sources (with no background) with two different spectral shapes — a single power law 
(a) and power-law with an exponential cutoff at 1 GeV (b) — for several different values of 
the input index. Enabling the energy dispersion handling in the ScienceTools significantly 
reduces the bias in the fit parameters. 


142 


The situation can be quite different for other spectral shapes — e.g., a power-law spec- 
trum with exponential cutoff, as illustrated in Figure [76] (b). Particularly, the bias introduced 
by the energy dispersion reaches in this case ~ 0.15 for hard spectral indices. The problem 
can be studied, for a specific source model, by means of dedicated MC simulations, and 
this is strongly recommended, especially for spectral analysis below 100 MeV. Generally 
speaking, enabling the energy dispersion handling in the standard binned likelihood analysis 
significantly reduces the bias in the fit parameters. 


To quantify the magnitude of the bias introduced by ignoring energy dispersion in the 
likelihood fit we have re-analyzed both the B2 1520+31 and the PG 1553+113 ROIs with 
the energy dispersion treatment enabled using the ScienceTools (version v9r28p0). We list 


the changes with respect to the fit results with the energy dispersion ignored in Table 22 


Parameter 

B2 1520+31 
(soft) 

PG 1553+113 
(hard) 

SN 0 /N 0 

+3.9% 

+2.2% 

5T 5a 

-0.04 

+0.01 

sp 

+0.02 

- 

5F 25 /F 25 

+1.4% 

+3.3% 

$s 25 /s 25 

+2.0% 

+1.1% 


Table 22: Systematic variations arising from ignoring the effects of energy dispersion when 
performing likelihood fitting. For the spectral index (T or a) and spectral curvature (/3) we 
give the absolute variation with respect to the value obtained when ignoring the energy dis- 
persion (e.g., hr). For the flux prefactor and the integral fluxes we give the relative variation 
with respect to the value obtained when ignoring the energy dispersion (e.g., 5N 0 /N 0 ). 


MC simulations also show that the effects of systematic uncertainties in the energy 
resolution, at the level we understand the energy dispersion of the detector, are essentially 
negligible over the entire energy range when studying y-ray sources that do not have sharp 
spectral features — which accounts for the vast majority of cases of practical interest. One 
noticeable exception is the search for y-ray lines, which requires a dedicated analysis when 
evaluating limits. But even in that case the systematic uncertainty on the energy resolution 
is not a major source of concern: if the actual energy resolution was ~ 10% broader than 
that predicted by the MC the fitted signal counts would be ~ 10% lower, i.e. , not enough to 
dramatically decrease the sensitivity to a spectral line. 

The effect of systematic uncertainties on the absolute energy scale is also strongly de- 
pendent on the energy range of interest. In the ideal case of an energy-independent effective 
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area, it can be shown that, for a power-law spectrum with index h, a relative bias 6 in the ab- 
solute scale translates into a rigid shift of the spectrum itself by the amount A F/F = (T — 1)6 
(or A S 25 /S 25 — (r — 2)6 if we consider the integral energy flux). In short, the measurement 
of the spectral index is not affected — while the flux obviously is. 


As a concrete example, we estimate the effect of systematic uncertainties of the abso- 
lute energy scale on the integral fluxes i+5 and S25 of PG 1553+113 and B2 1520+31 by 
integrating the fitted source spectra and shifting the limits of integration by the uncertain- 
ties stated in § |7.3.4[ Specifically, we consider the integration ranges 96 MeV-95 GeV and 


104 MeV-102 GeV: the resulting uncertainty estimates are reported in Table 23 


Parameter 

B2 1520+31 

PG 1553+113 


(soft) 

(hard) 

SF 25 /F 2 5 

+13.4% -4.6% 

+3.4% 

-2.5% 

ss 25 /s 25 

+3.8% -1.4% 

+1.7% 

-0.6% 


Table 23: Systematic variations of integral fluxes arising from uncertainties in the absolute 
energy scale. We give the relative variation with respect to the nominal value (e.g., SF 25 /F 25 ). 


Along the same lines, measured cutoff energies for sources with curved spectra reflect the 
bias in the energy scale directly: A E c /E c = 6. These statements are useful rules of thumbs 
to estimate the order of magnitude of the effect above ~ 1 GeV, where the effective area 
is not strongly dependent on energy (Fig. 77). In fact, since we have good indications that 
we understand the energy scale of the detector at the ~ 5 % level (see § 7.3.1 and § 7.3.3) 
at these energies, the effect of systematic uncertainties in the scale is generally smaller than 
other sources of systematics in this energy range, and does not affect the measurement of 
the spectral index in the power-law case. 


Below a few hundred MeV the systematic uncertainties in the absolute energy scale 
should be carefully considered in any analysis due to the unfortunate combination of a 
steeply increasing effective area and a worsening of the energy resolution at low energies. As 
for the energy dispersion the effect is more pronounced for steep sources. For the typical 
case of a spectral index of Y « 2, a ±5% uncertainty in the energy scale is negligible down 
to ~ 100 MeV but at 50 MeV can easily translate into an additional 10-20% bias in the 
count spectrum. Again, the order of magnitude of the effect can be easily assessed on a 
case-by-case basis, depending on the spectral shape of the source under study and energy 
range of interest. 

Spectral unfolding (deconvolution) is an approach for taking into account the energy 
dispersion in the spectral fitting. Unfolding is in effect the opposite of the forward folding 


144 



Fig. 77. — Bias factor on the count spectra induced by systematic uncertainties on the 
absolute energy scale for a simulated point source (with no background) with power-law 
spectrum with T = 2. The parameter b indicates the relative bias in the absolute energy 
scale; the case b = 0 corresponds to the line at V = 2 in Figure [75} 


approach implemented in many popular astronomical data analysis packages (and used in 
the standard likelihood point source analysis implemented in the ScienceTools ), where the 
IRFs are folded with the spectral model and the model parameters are varied until the best 
match with the measured count spectrum is found — e.g., this is how XSPEC (Arnau d||199 6) 
handles the energy dispersion. 

Energy unfolding has been used for several analyses of LAT data, mainly to cross-check 
the results of the standard likelihood analysis (Abdo et al. 2010g) and especially for the lowest 


or highest energies: below r\j 100 MeV (Abdo et al. 2010b) or above 


et al. 2010b; Abdo et al. 2009e). The Bayesian approach detailed in D’Agostini (1995) is 


100 GeV ( 

Ackermann 

D’Agostini 

(1995 

) is 


that typically being used by the LAT Collaboration. 


7.5. Event Analysis-Induced Spectral Features 


Though we discuss potential event analysis-induced spectral features in the context of 
energy reconstruction, we emphasize that the effective area is also germane to this topic. 
We conclude that the event analysis and IRFs do not introduce significant artificial spectral 
features because bright y-ray sources are devoid of spectral features — particularly the Earth 
, which is the single brightest source in the LAT energy range. 


limb (§ 3.6.3 


Figure 78 shows a count spectrum of P7CLEMLV6 7 rays from the Earth limb based 
on the first two years of data. The spectrum is made with ~ 16 bins per energy decade — 
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Fig. 78. — Count spectrum of the Earth limb based on the first two years of data, including 
P7CLEAN_V6 events within three times the acceptance-averaged PSF 68% containment from 
the nominal zenith angle of the limb. Since the limb itself is typically observed at large 
incidence angle in the instrument frame, events impinging on the LAT more than 65° off- 
axis have been removed in order not to overweight the edge of the FoV relative to typical 
analyses of celestial sources. The turnover below ~ 1 GeV is mostly due to the decrease of 
the LAT acceptance, rather than to an intrinsic roll-off of the source spectrum. 


corresponding to a bin width slightly smaller than the typical LAT energy resolution — and 
its smoothness is a good qualitative indicator of the smallness of any possible instrument- 
induced spectral feature. 


In order to quantify this we compared the counts n t in the i-th energy bin with the 
value fi returned by log-parabola fit to the four nearest bins (i.e. , two bins on each side). 
We can construct two obvious metrics, namely the normalized residuals: 


n _ (m - fi) 

and the fractional residuals 

j (ni ~ fi) 
‘ f, 


(54) 


(55) 


Figure 79 shows a scatter plot of the two metrics, each point corresponding to an 
energy bin in the count spectrum in Figure [78] There are two striking features: (i) a nearly 
horizontal branch corresponding to the energy bins below ~ 1 GeV, where the statistics 
are large and the curvature of the count spectrum makes the quadratic fits less accurate 
and (ii) a nearly vertical branch, corresponding to the high-energy points, whose relatively 
larger fractional deviations (up to ~ 10%) are just due to statistical fluctuations. We note 
for completeness that there are no points farther than 5 a from the fit value above 2 GeV, 
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Fig. 79. — Single-point residuals of the count spectrum in Figure [78] with respect to a log- 
parabola fit of the four nearest points. 


which implies that we do not expect to be able to see artificial spectral features above this 
energy in any practical case. At lower energies the high counting statistics of the Earth limb 
emission allows placing an upper limit on possible narrow (i.e., less than about twice the 
energy resolution) spurious features at the level of a few percent. This method for setting 
limits on spurious line features is decreasingly sensitive to broader features, which are more 
likely to arise from correlated errors in A e g 


see 


5.6.2). 


7.6. Comparison of Pass 6 and Pass 7 


Though the underlying energy reconstruction algorithms are exactly the same, Pass 7 
differs from Pass 6 in (i) the new un-biasing stage described in § 3.2.1 and (ii) the fact that 


the LH estimator is no longer used in the event-level analysis for the final energy assignment 


see 


3.3.2). 


Figure [80] shows that the energy resolution for the P7S0URCE_V6 event class is quite 
comparable to that of P6_V3_DIFFUSE over most of the LAT energy range. The most notice- 
able difference, namely a slight worsening below 1 GeV (and especially below 100 MeV), is 
a small trade-off for the much greater Pass 7 low-energy acceptance. On the other hand 
Pass 7 has uniformly better energy resolution than Pass 6 in the energy range above 10 
GeV. 
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Fig. 80. — Comparison of the on-axis energy resolutions of the P7S0URCEJV6 

P6_V3_DIFFUSE IRFs. 


and 
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8. PERFORMANCE FOR HIGH LEVEL SCIENCE ANALYSIS 

We report here the high-level science analysis performance of the LAT for the Pass 7 
event selections. 


8.1. Point Source Sensitivity 


A detailed description of the point-source sensitivity of the LAT is given in Abdo et ah 


(2010a) and Nolan et al. (2012). In particular, following the procedure described in Abdo 


et al. (2010a), a semi-analytical estimate of the LAT sensitivity for point sources can be 


calculated for the Pass 7 event analysis. Here we calculate the LAT sensitivity to a point 
source for the P7S0URCE_V6 IRFs — as described in Abdo et al. (2010a) for P6_V3 — under the 
following assumptions: 


• power-law spectrum with T = 2; 

• diffuse Galactic emission as described in the 2-year Galactic diffuse model, publicly dis- 
tributed as gaL2yearp7v6_v0.fits ; 

• diffuse isotropic background as described in the 2-year template, publicly distributed as 
iso jplv6source.txt ; 

• no confusion with nearby sources; 

• exposure calculated for the first 3 years of nominal science operations. 

A map of the flux limit is shown in Figure [8TJ Away from the Galactic plane the 
sensitivity is rather uniform. In reality, the presence of bright point sources will affect the 
flux limit and a dedicated analysis is recommended to evaluate the limit for any particular 
circumstance. 

In Figure [82] we show the corresponding differential sensitivity curves. The curves show 
the flux limits for narrow energy ranges and illustrate the sensitivity for spectral measure- 
ments as a function of energy. We require that in each energy band the source be bright 
enough to have 5cr detection and cause at least 10 7 rays to be collected. 

Finally, we stress that along the Galactic plane the uncertainties of the Galactic diffuse 
emission can affect the detection significance of a source, and structured residuals in the 
Galactic diffuse model can be mischaracterized as point sources. Further discussion can be 
found in Fermz-LAT point source catalogs and papers about the Galactic diffuse emission 


(Abdo et al. 

2010a 

Nolan et al. 

2012 

Ackermann et al. 

2012b) 
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Fig. 81. — Flux above 100 MeV required for 5 a sensitivity for the P7S0URCE_V6 event class 
for a point source with power-law spectrum with index T = 2. The calculation assumes a 
3-year exposure. The entire sky (a) and a zoom on the Galactic center (b) are shown. 



Fig. 82. — Differential sensitivity for class P7S0URCEJV6 for a point source; calculation for 
3-year exposure, 4 bins per energy decade. Requirements are 5a sensitivity and at least 10 
counts per bin. The sensitivity is calculated at three locations in the sky: at the Galactic 
pole, at an intermediate latitude and on the Galactic plane. 


8.2. Point Source Localization 


The performance of the LAT for the localization of a point source at the detection 
threshold also can be evaluated with a semi-analytical approach. The estimated 2-year 
performance (see Nolan et al. 2012) indicates that the 95% localization radius for an isolated 
point source detected at 5a significance at high Galactic latitude ranges from ~ 0.1° for a 
hard spectral index (T = 1.5) to ~ 0.3° for a soft spectral index (T = 3.0). 

Scaling the source location determinations to more intense background levels is not 
straightforward. First of all, localization regions are usually elliptical in shape, as described 


m 


Nolan et al. (2012). Secondarily, as the astrophysical background increases, e.g., in the 


Galactic plane, for a given value of the spectral index a higher flux is needed to reach 
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the detection threshold. As a consequence the size of the localization region, which is more 
sensitive to the number of high-energy 7 rays than the point source sensitivity, will be smaller 
within the Galactic plane. This apparently counter-intuitive result is shown in Figure [83] 
where the 95% localization radius for a power-law source is shown as a function of the 
spectral index for three different locations. The improvement of the localization radius with 
time (Figure [83]) is of the order of 10% or less going from a 3-year to a 5-year exposure. 


Nolan et al. (2012) compared the positions and error regions of 7 -ray sources to the 
positions of associated multiwavelength counterparts and empirically corrected the localiza- 
tion uncertainties by multiplying the error region by a 1.1 scale factor and adding 0.005° in 
quadrature to the 95% error ellipse axes. Although the uncertainties of the Galactic diffuse 
emission are significantly larger along the Galactic plane, we did not observe any systematic 
increase in the localization offsets with respect to associated multi wavelength counterparts 
in the Galactic plane. 



Fig. 83. — 95% localization radius for class P6S0URCE_V6 for a point source at the detection 
threshold with power-law spectrum, as a function of the spectral index (3-year exposure 
for different locations in the sky). For sources in the Galactic plane detected at a given 
significance, the localization performance for those sources which are softer than the Galactic 
diffuse background (which consists of spectral components V G [2.1, 2.7]) is significantly 
degraded with respect to harder sources. 


8.3. Flux and Spectral Measurements 

The precision to which we can measure the flux and spectral index of any source depends 
primarily on the counting statistics, which in turn depends on the flux of the source. However, 
since the width of the PSF decreases with energy, high-energy 7 rays contribute more to 
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source detection, so given similar energy fluxes, we can measure harder sources more precisely. 
Finally, the flux of the source relative to nearby diffuse backgrounds and other nearby sources 
will also limit the precision of the parameter measurements in a likelihood analysis. Figure [84] 
shows the dependence of the statistical uncertainties of spectral index T and the integral 


energy flux between 100 MeV and 100 GeV (S25) for the sources in the 2FGL catalog (Nolan 


et al. 2012). In both cases the minimum statistical uncertainty decreases as roughly S- 


-0.6 
25 • 
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Fig. 84. — Statistical uncertainties on spectral index T (a) and integral energy flux above 


100 MeV S 25 (b) as a function of S 25 for all sources in the 2FGL catalog (Nolan et al. 2012). 
The dashed lines show a S^ 5 0 ' 6 dependence for comparison. 

For comparison, we combine the systematic uncertainty estimates based on the bracket- 
ing tables for the effective area (§ 5.7), the PSF (§ 6.5), and the effect of ignoring the energy 
dispersion (§ 7.4) for the two sources we studied in detail (B2 1520+31 and PG 1553+113). 
Since we have not seen any evidence of correlated biases in the different terms of the IRFs 
(see, for example § 7.1.4) we simply combine the uncertainties in quadrature, as summarize 


in Table 24). 


Quantity 

A e ft 

PSF 

Energy 

Dispersion Scale 

Total 

<5rpG1553+113 

0.05 

0.02 

0.01 

0.05 

B21520+31 

0.09 

0.07 

0.04 

0.12 

SS 25 

S 25 PG1553+113 

10% 

5% 

1% 2% 

11% 

<?-S l 2 5 

S 25 B21520+31 

8% 

6% 

2% 4% 

11% 


Table 24: Rough estimates of the magnitude of the effects of various sources of systematic 
errors on the integral energy flux and spectral index of B2 1520+31 and PG 1553+113. 


Figure 85 shows the distributions of the statistical uncertainties on T and S 25 from the 
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2FGL catalog, along with our estimates of the typical systematic uncertainties for hard and 
soft sources based on our analysis of B2 1520+31 and PG 1553+113. For the majority of 
sources in the 2FGL catalog the measurement precision is still limited by the statistical 
uncertainties, though for the brightest sources the systematic uncertainties dominate. For 
both T and S 25 the transition from statistical limitation to systematic limitation occurs for 
S -25 ~ 10 -11 erg cm -2 s -1 . 




Fig. 85. — Statistical uncertainties on V (a) and S 25 (b) for all 2FGL sources. For comparison 
the estimated systematic uncertainties for B2 1520+31 and PG 1553+113 are also shown. 


We emphasize that the results in this section were derived using bracketing IRFs de- 
signed to maximize the variation of the fit results for the particular sources under study. 
As such, they represent the systematic uncertainties of measurements on any single source. 
Since measurements of different sources share the same IRFs and associated uncertainties, 
any relative comparison between measurements is significantly more precise. 


Finally, along the Galactic plane the uncertainties of the Galactic diffuse emission can 
affect measurements of source parameters and spectral indices. As with the question of how 
these uncertainties affect the LAT point source detection sensitivity, detailed discussion is 
beyond the scope of this paper, but can be found in Abdo et al. (2010a), Nolan et al. (2012), 
and Ackermann et al. (2012b). 


8.4. Variability 

The LAT’s ability to measure the variability of any given source depends on the char- 
acteristics of both the source and the astrophysical backgrounds. 

We can use the detection significances of sources in the 2FGL catalog to derive a rough 
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estimate of the time-scales at which we can probe the variability of those sources. For 
detecting variability, the worst-case scenario is a steady source, so we ask how long would be 
required on average to detect any source at a particular threshold, assuming steady emission. 
Specifically, we can construct a metric r 2(r , the time required to achieve 2 a detection of a 
source, assuming steady emission at the average 2FGL level, which can be expressed in terms 
of the significance of the source in the 2FGL catalog (ct 2 fgl) and the amount of data used 
to construct the catalog (~ 730 days): 


r 2a ~ 4cr 2 p GL 730 days. (56) 

Figure [86] shows this estimator for all 2FGL sources, we emphasize that we are sensitive to 
more significant flux changes on timescales shorter than r 2a . 
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Fig. 86. — LAT sensitivity to variability of all 2FGL catalog sources: panel (a) shows the 
cumulative distribution of the estimator t 2ct , an estimate of the time needed to detect a 
steady source with 2 a significance; panel (b) shows r 2(J as a function of the integral energy 
flux between 100 MeV and 100 GeV (S 2 5) for both low- and high- Galactic latitude sources. 


we examined how instrumental uncertainties in the A e s and PSF 
can change over time, and could potentially induce artificial variability in the measurement 
of point source fluxes. 

More specifically, in the case of the A e ff-related variability we showed that the phase- 
gated counts excess was extremely stable, with the Fourier transform of Sn — n — h showing 
only a single small feature consistent with the 53.4-day orbital precession period. However, 
since that analysis did not treat the PSF, nor allow for the effect of changing statistics from 
differing exposures on the Fourier transform of 5n , we could not use it to predict the level 
of instrumentally induced variability we might measure in a likelihood-based analysis using 
the ScienceTools. 


In § |5.6.3 and § 6.5.3 
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Similarly, our analysis of PSF-related variability estimated how much the containment 
as defined by the ^-averaged i ? 68 and R 95 might vary over a precession period, but did not 
attempt to quantify the effect on a likelihood based analysis using the ScienceTools. In order 
to quantify these effects we have studied the Vela and Geminga (Abdo et al. 2010d, PSR 


J0633+1746) pulsars in 12-hour time bins with the ScienceTools (version 09-26-02) unbinned 
likelihood analysis and the following analysis parameters: 


1. we started with the y-ray and time interval selection criteria as for the Vela calibration 
sample (see 


3.6); for the Geminga sample we used the same 15° radius for the ROI; 


2. tightened the zenith angle requirement to 9~ < 95° and used all 7 rays with energies 
E > 70 MeV; 


3. included all 2FGL sources within 20° in our likelihood model, with same spectral parametriza- 
tions as were used in the 2FGL catalog; 

4. included models of the Galactic diffuse emission ( ring-2year.P76-vO.fits ) and the isotropic 
diffuse emission ( isotrop-2year_P76source-vl.txt ); 

5. held fixed all of the parameters of the model except for the normalization of the flux of 
the pulsar under test and isotropic diffuse emission. 


We performed this study with both the P7S0URCEJV6 and P7S0URCEJV6MC IRFs, and 
included the </> dependence of the A eS in the exposure calculations. In all cases we measured 
the pulsars to be almost consistent with having constant fluxes during the first 700 days of 
the mission. Specifically, for Geminga P(y 2 ) — 0.04, while for Vela P(y 2 ) = 6.3 x 10 -5 , but 
adding 2% error in quadrature brought the P(y 2 ) up to 0.17. Furthermore, no peak was 
visible in the Fourier spectra of the fitted fluxes between 12 hours and 700 days. On the other 
hand, the fitted normalization for the isotropic diffuse component was not consistent with 
constant, and in both cases the Fourier analysis showed a significant peak at the 53.4 day 
orbital precession period (see Fig. 87 for results of the Fourier analysis of the Vela flux 
and associated isotropic normalization). The variability in the normalization of isotropic 
background is caused by the CR background leakage not having the same 6 dependence as 
the 7 -ray effective area, so that exposure calculations used to predict the expected counts 
from the isotropic background template suffer slightly different biases during different phases 
of the orbital precession period. As we have shown here, this can be handled by leaving the 
normalization of the isotropic template free, provided that the source dominates the nearby 
isotropic background. 


By comparing results obtained on the Vela pulsar with the P7S0URCE_V6MC 
IRFs (which include the 6 dependence of the PSF) to results obtained with the 
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Fig. 87. — Fourier transforms of the normalized residuals of the Vela integral counts flux 
(a) and the isotropic normalization factor (b). Note that the figures are normalized and the 
vertical scale is expressed in units of the statistical uncertainty. 


P7S0URCE V6 IRFs (which do not), we estimate that ignoring the 9 dependence of 
the PSF causes a ~ 4% RMS variation of the flux of Vela when fitted in 12-hour 
time intervals. The effect decreases on longer time scales. Since, as stated in 
§ 6.1| , the PSF depends much more strongly on 9 than on 0, we neglect the 0 
dependence of the PSF as a potential source of instrument-induced variability. 

In summary, when the angular dependence of the IRFs is properly accounted 
for, and time variations of the CR-background leakage are absorbed into the nor- 
malization of the isotropic component in the likelihood fit, we find the level of 
instrument-induced variability to be small (< 5%) for all time scales between 12 
hours and 2 years. Since the LAT boresight follows very similar paths across the 
sky every two orbits, resulting in similar observing profiles, we believe these sys- 
tematic uncertainties are also applicable at the 3-hour time scale. On the other 
hand, for observations that do not consist of a complete orbit, and which are 
therefore more susceptible to biases in the IRFs at specific incidence angles, the 
(~ 10%) uncertainties in the effective area quoted in § 5.6.1 are more applicable. 


Finally, we note that the slow change in the light yield of the Csl logs of the 
CAL, and the corresponding change in the energy scale described in § |7.3.2 


can 


induce slow shifts in observed fluxes by shifting spectra as a function of time. 
The effect depends on the spectrum of the source under study, as well as on the 
energy band used, and differs for counts fluxes and energy fluxes, as summarized 
in Table I 
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9. SUMMARY 


Since the beginning of science operations in 2008 August, the Fermi-L AT has opened 
a new window on the 7-ray sky, providing the science community with an unrivaled set of 
7-ray data. The LAT data set covers the entire sky over the energy range from ~ 20 MeV 
to greater than 300 GeV every 3 hours. 


During these first years of the mission, the LAT Collaboration has studied the on-orbit 
performance of the LAT, the character of the flight data from the LAT, and the optimization 
of event selections and IRFs for science. This has led to significant improvements in the event 
analysis (in particular the Pass 7 version of the event analysis described in § 3.3) and in 


the science analysis, and to important clarifications of the systematic uncertainties at each 
level of the analysis. In addition, the procedures we designed and developed to perform 
these investigations can be reiterated with minimum modifications whenever a new event 
reconstruction and classification is adopted. 


We have shown that the LAT has performed extremely well: the data are of uniformly 
high quality, we have lost less than 0.5% of potential observing time (i.e, when Fermi is 
outside the SAA) to unplanned outages and instrumental issues. Furthermore, the LAT 
response has been extremely stable: we are able to use a single set of IRFs for the entire 
mission to date, with any variations in performance contributing negligibly to systematic 
uncertainties (§[2]). 

We have also confirmed that the LAT data reduction and analysis tools have performed 
very well. The combination of a configurable hardware trigger, an on-board event Liter, and 
ground based event reconstruction and analysis have allowed us to reduce the CR contam- 
ination from 5-10 kHz passing through the instrument to ~ 1 Hz in the P7S0URCE event 
sample (and ~ 0.1 Hz in the P7CLEAN event sample, §[3]) while maintaining a peak acceptance 
of over 2 m 2 sr in the P7S0URCE sample (and > 1.75 m 2 sr in the P7CLEAN sample) (§ [H]). 


We have validated the quality of our MC simulations and found only one significant 
discrepancy between the MC simulations and the flight data. Specifically, limitations in the 
pre-launch calibration algorithm of the CAL light asymmetry produced calibration constants 
that did not match the MC predictions of the spatial resolution performance of the CAL. 
However, we have quantified the effects of this discrepancy on the A eS (§ 5.4) and PSF (§ 6.2) 


and are currently assessing the improvement of the PSF for data that were reprocessed with 
improved calibration constants. 


Pass 7 data have been available for public analysis since 2011 August and provide sub- 
stantial improvements over Pass 6, primarily due to greatly increased A e g below ~ 300 MeV 


5.8) and improved modeling of the IRFs. Coupled with improved understanding of 
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the effects of energy dispersion (§ 7) this increase in effective area is opening a 
new window for analysis of LAT data for energies below 100 MeV. Table 25 lists 
the features of the components included in each of the IRF sets recommended for analysis 
with both Pass 6 and Pass 7. 



IRF Set 

A e ft 

PSF 

P 6 _V3 -TRANS I ENT 

flive, 0 

Legacy MC 

P6_V3_DIFFUSE 

flive ? 4* 

Legacy MC 

P6_V3_DATACLEAN 

flive i 4* 

Legacy MC 

P6_V11 .DIFFUSE 

flive, 4>, in-flight 

in-flight 

P7TRANSIENT_V6 

flive , 4* 

in-flight 

P7S0URCE_V6 

flivei 4* 

in-flight 

P7CLEAN_V6 

flivei 4* 

in-flight 

P7ULTRACLEAN_V6 

flive i 4* 

in-flight 

P7S0URCE_V6MC 

flive i 4* 

MC 

P7CLEANJV6MC 

flive 5 4* 

MC 


Table 25: IRF sets recommended for data analysis, and the features included in the various 
components. For A eff , f Uve and 4> refer to the live time fraction and </> dependence, respectively 
(§ 5.2.2 and § 5.2.3), and “in-flight” refers to flight-based corrections (§ 5.4). For the PSF 


Legacy MC” refers to an early version of the parametrization fit to allGamma samples 
6.1.2), “in-flight” refers to the flight-based PSF (§ 6.2), and “MC” refers to the more 


recent parametrization fit to allGamma samples (§ 6.1.1). Finally, we have used the same 


parametrization of the energy dispersion for all of the IRF sets (§ 7.1.2). 


In addition to the caveat: 


14 


and the documentation^] accompanying the data and the 
science analysis software, we provide with this paper a reference document for all currently 
known systematic issues. Table |26] provides numerical estimates of the residual uncertainties 
and refers to the sections of this paper where we detail procedures to estimate the systematics 
uncertainties for many analyses. 

Finally, in §|8j we provide details on the science performance we obtain with the Pass 7 
event analysis and IRFs. In particular, we provide estimates of the source detection sensi- 
tivity threshold (§ 8.1), the source localization performance (§ 8.2), the expected precision 


14 http://fermi. gsfc.nasa.gov/ssc/data/analysis/LAT _caveats.html 
1 J http: // fermi. gsfc . nasa.gov/ssc /data/ analysis / documentation / 
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Quantity A eff PSF Energy 

Dispersion Scale 

F25 ~ 8% (§ 

S*25 ~ 10% (§ 

r ~ 0.09 (§ 

Variability ~ 3% (§ 

Localization - 

5.7) 

5.1 

+7 

5.6) 

~ 8% (§ 
7 ) ~ 6% (§ 

) ~ 0.07 (f 

~ 3% (§ 
~ 0.005° 

6.5 

6.5 

6., 

6.5 

(§ 

~ 3% (§ 
~ 2% (§ 
5) ~ 0.04 (§ 

8.2 ) a - 

7.4 

7.4 

7/ 

+13% -5% (§ 7.4) 
+4% — 2%(§ 7.4) 

1) - 


Tabic 26: Rough estimates of the magnitude of the effects of various sources of systematic 
errors for commonly-measured 7-ray source properties. We also provide references to the 
relevant sections with more details. 


“See Nolan et al. (2012) for a discussion of the systematic uncertainties on source localization. 


and accuracy of measurements of fluxes and spectral indices 
accuracy of variability measurements (§ 0- 


|8.3), and the precision and 


As stated in § [lj the LAT team will continue to make both major improvements and 
minor refinements to many aspects of the event reconstruction, analysis and to the associated 
IRFs. We will continue to keep the 7-ray astronomy community informed of the state-of- 
the-art of our understanding of the LAT and issues relating to analyzing LAT data. 
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A. SOLAR FLARES AND BAD TIME INTERVALS 


With the beginning of the solar activity connected with Cycle 24, in early 2009, it 
became clear that Solar Flares (SFs) can cause significant X-ray pile-up in the ACD, clearly 
visible in many low-level ACD quantities (see Fig. 88). The on-board ACD veto electronics 
(see Table [TJ are fast enough that this is not causing any change in the efficiency of the 
trigger and/or that of the GAMMA filter (i.e., there is essentially no loss of 7 rays on-board). 
However the effects of the pile-up are potentially much more severe at the level of the event 
selection on the ground — where the slow signals are used — as the additional activity in the 
ACD can cause 7 rays to be misclassified as background. 



UTC 


Fig. 88. — (a) X-ray flux for the XI. 5 solar flare on 2011 March 9 measured by the GOES-15 
satellite in the energy bands 0. 5-4.0 A (3.1-24.8 keV) and 1. 0-8.0 A (1.5-12.4 keV). The 
latter is customarily used to classify solar flares based on their X-ray radiance, as indicated 
by the letters A, B, C, M, X on the right. The effect on two of the basic ACD quantities 
(see text for more details) is clearly visible in panels (b) and (c). 


The basic phenomenology is somewhat similar to the ghost-induced loss of efficiency 
introduced in § |2.1.4| and discussed in detail throughout the paper, the main difference being 


that during the most intense SFs the effect can be large enough to make the LAT essentially 
blind to 7 rays. From the point of view of the data analysis the most relevant implication 
is that there are time intervals in which the standard IRFs do not accurately represent 
the detector and therefore the results of the likelihood analysis are potentially unreliable. 
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While the LAT collaboration is considering possible modifications to the event reconstruction 
and selection aimed at mitigating the problem, these Bad Time Intervals (BTIs) are being 
systematically identified and flagged. 


Operationally a BTI is characterized by a suppression of the rate of events in the stan- 
dard 7-ray classes. The 7-ray rates intrinsically feature large orbital variations depending on 
both the geomagnetic environment (through the background leakage and the ghost effect) 
and the rocking angle of the observatory (through the change of the arc length of the Earth 
limb in the FoV). These variations can be parametrized to an accuracy of 20-30% and in 
fact are accounted for (at a similar level of accuracy) in event rates normalized to predicted 
values that are routinely accumulated in 15 s t im e bins for data monitoring purposes (see 
e.g., Figure 89). Still, the residual variations (especially during rocking maneuvers or non- 
standard pointings) make the normalized rates not directly suitable for identifying the BTIs. 
We use a two-step procedure instead, in which: 


(i) within each orbit we search for evidence of significant X-ray pile-up activity in the ACD 
in order to define a list of candidate BTIs (see Figure [88]) ; 


(ii) we search for a 
those intervals (see 


temporally coincident decrease in the transient 7-ray rate within each of 
Figure |89l). 



22:30 


23:00 


23:30 


UTC 


Fig. 89. — (a) Normalized rate of P7TRANSIENT events for the same time interval shown 
in Figure 88 and (b) integrated time loss over the candidate BTI in Figure 88 (b). Each 
data point corresponds to a 15 s time bin. The solid line is a second order polynomial fit 
to the normalized rate outside the candidate BTI and it is used to calculate the integrated 


time loss, as defined in Eqn. (Al). 


As shown in Figure 


X-ray pile-up induced by bright flares is typically visible both 
in the normalized ACD hit multiplicity and in the single tile hit occupancy — particularly in 
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tile 63, which is the largest tile on the +x side of the observatory (i.e., the side facing the 
Sun). Though the correlation with the X-ray radiance measured by GOES is far from perfect 
(with a zero-suppression threshold of ~ 100 keV the ACD is sensitive to somewhat higher 
energies than GOES), both quantities are good proxies for the increase in the X-ray flux. 
Figure |89| (a) shows the corresponding large suppression of the normalized P7TRANSIENT rate 
for the same solar flare. As shown in Figure 89 (b), by fitting this normalized rate outside 


the candidate BT1 it is also possible to define an integrated time loss Ti oss corresponding to 
a particular candidate BT1 


Tioss = > Xfi - Ti)At 


(Al) 


BTI 

(where _/) and 77 are the fitted and measured normalized rates in the i-th time bin, respec- 
tively, and At is the bin width, i.e., 15 s). Roughly speaking, if Ti oss is larger than a few 
minutes we will mark the relevant time interval as bad. 


All LAT y-ray data automatically have the DATA_QUAL field in the spacecraft pointing 
and live time history hie (spacecraft file)P^| set to 1 by the data processing system and are 
immediately exported to the FSSC where they become publicly available. The potential for 
BTIs is reviewed only after the fact. This is because such a high percentage of the LAT data 
is good and making them publicly available as quickly as possible is a priority. 



Fig. 90. — Time history of the sum of BTIs induced by SFs. The sum of time intervals 
marked as BAD as of March 2012 is ~ 600 minutes (i.e., ~ 10 hours). 

For each flare the entire information available is reviewed manually and if there is a 
significant loss of efficiency the corresponding time period is marked as BAD (specifically, 
the DATA_QUAL field is set to —1) and a new spacecraft file is generated with the appropriate 
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http: / /fermi. gsfc.nasa.gov / ssc / data/ analysis / documentation / Cicerone / Cicerone_Data/LAT_Data_Columns.html#SpacecraftFile 


data quality flag in the corresponding 30 second time bins. This new file is then exported 
to the FSSC and supersedes the original file. As of 2012 March the sum of time intervals 
marked as BAD is ~ 10 hours (Fig. 


90 ). 
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B. LIST OF ACRONYMS AND ABBREVIATIONS 

Fermi LAT Second Source Catalog 

LAT Anti-Coincidence Detector subsystem 

Active Galactic Nucleus 

Bad Time Interval 

Blind Search Pattern Recognition 

LAT imaging Calorimeter subsystem 

Carbon, Nitrogen, Oxygen [CR species] 

Charged Particle in the Field of view analysis 

Cosmic Ray 

Calorimeter-Seeded Pattern Recognition 
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LAT Calibration Unit 

Energetic Gamma-Ray Experiment Telescope 
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Fermi Science Support Center 
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Instrument Response Function 

Instrument Science Operations Center 
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Maximum LikeliHood [energy estimation algorithm] 

Monte Carlo 
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Minimum Ionizing Particle 

Parametric correction [energy estimation algorithm] 

photo-electron 
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LAT Tracker /converter subsystem 


2.1 
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C. NOTATION 


A e ft 

B ■, -B a , H ( ] x 
B, B 68 , B r 
b 
b 

(l,b) 

^68 

Ci 


Ci 

D 

<U 

E 

F 25 

Ft 

/core 

K 

k 

L 

Lext 

Lpt 

M 

No 

Ngen 

n 

n 

P 

P 

-Pall 

PfcAL 

P 

1 core 

-PcPF 

Pe 

Ptkr 

P 

R 

R 


Effective area 

Bracketing functions for A e g- 

Bracketing functions for the PSF 

Bias in the energy scale 

Distribution of CR background events mischaracterized as 7 rays 

Galactic coordinates 

68% containment radius of the PSF 

Actual containment at a given nominal containment radius 

Generic fitting constants 

Energy dispersion parametrization 

Normalized deviations 

Energy 

Integral counts flux between 100 MeV and 100 GeV 

Live time fraction 

Fraction of counts in the “core” King function for the PSF 

King function [used to model the PSF] 

A e fj bracketing smoothing constant 

Mcllwain L 

Log likelihood for fit as an extended source 

Log likelihood for fit as an point source 

Predicted counts distribution of 7 rays 

Flux prefactor 

Number of events generated [when making IRFs] 

Number of observed counts 

Expected counts 

Point-spread function parametrization 

PSF averaged over the observing profile 

Combined estimator that the event is a 7 ray 

Estimator that the event is a 7 ray, based on the CAL topology analysis 

Estimator of the quality of direction measurement 

Estimator that the event is a 7 ray, based on the CPF analysis 

Estimate of the quality of energy measurement 

Estimator that the event is a 7 ray, based on the TKR topology analysis 

Incident direction in celestial reference frame 

Rando function, used to model the energy dispersion 

Ratio of efficiency for flight data to MG simulated data 


■§0 

5XT] 

6X1 
§ [7T 
5 14~6 


3.3.3 

6X3 


••§0 

5X2] 

•■§0 

5XT] 

5X2 

6X1 

6X1 

5X1 

2X1 

6X2 

6X2| 

••§0 

5XT] 

§ [XT 
§ |XT 
5X3] 

••§0 

§ |6T 

3X8 

3X7 

3XX 

3X5 

3X2 

3X6] 

■•§0 

XX] 

5|5X 


167 


-^68 

r 

r 


' / 
s 

S25 

Sd 

S P 

s 

TS ext 

t obs 

V 

X 

X 

X 

a 

P 

P 

r 

7 

7 

e 

? Idata ) ^Imc 

(M) 

(^, 0 Z ) 

£ 

cr 

cr 

T n 


68 % containment radius of the PSF (also i? 95 ) 

Ratio of size of signal region to background region 

Ratio of 95% to 68 % containment radii of the PSF 

Normalized residuals 

Fractional residuals 

Source distribution of 7 rays 

Integral energy flux between 100 MeV and 100 GeV 

Energy dispersion scaling function 

PSF scaling function 

Light yield scaling function 

Test statistic for source extension 

Observing time 

Incident direction in the LAT frame 

Scaled angular deviation 

Scaled energy redistribution 

Rando function break point 

Spectral Index in “log-parabola” model 

Spectral curvature parameter in “log-parabola” model 

PSF energy scaling index 

Spectral index 

King function “tail” parameter 

Rando function “tail” parameter 

Relative uncertainty in A e g 

Efficiencies for flight data and MG simulated data 

Polar angle and azimuth of the incident direction in the LAT frame 

Zenith angle and Earth azimuth angle 

Unit vector along local 6,(j> directions 

Folded 0 

King function width parameter 

Rando function width parameter 

Correlation metric 


6.2.3 

5X1 

6X1 

• §|54 

• § 55 

■■§0 

5.7.1 

7X1 

••§0 

.§[53 

.m 


• ■§E 

• •§0 
6 . 1.1 

7XT 
7X2 
5X1 
XT 
6X1 
5X1 
6X1 
7X2 
.§[27 
5X1 
§ [XT 
§ I2T2 

. §^42 
5X3 
6X1 
7X2 
XX 


REFERENCES 


Abdo, A. A., et al. 2009a, Phys. Rev. D, 80, 122004 

— . 2009b, Phys. Rev. D, 80, 122004 

-. 2009c, ApJ, 696, 1084 

-. 2009d, ApJ, 707, 580 

— . 2009e, Physical Review Letters, 102, 181101 

— . 2009f, Astroparticlc Physics, 32, 193 

— . 2010a, ApJS, 188, 405 

— . 2010b, ApJ, 708, 1254 

— . 2010c, Physical Review Letters, 104, 091302 

— . 2010d, ApJ, 720, 272 

-. 2010e, ApJ, 712, 459 

— . 20101, Physical Review Letters, 104, 101101 

— . 2010g, ApJ, 716, 30 

— . 2010h, ApJ, 713, 154 

Ackermann, M., et al. 2010a, Phys. Rev. D, 82, 092004 

— . 2010b, Phys. Rev. D, 82, 092004 

— . 2012a, ApJ, 745, 144 

— . 2012b, ApJ, 750, 3 

— . 2012c, Astroparticlc Physics, 35, 346 

— . 2012d, Physical Review Letters, 108, 011103 

Ackermann, M., et al. 2012, in preparation 

Agostinelli, S., et al. 2003, Nucl. Instrum. Meth., A506, 250 

Arnaud, K. A. 1996, ASP Conf. Series, 101, 17 



- 169 - 

Atwood, W. B., et al. 2007, Astroparticlc Physics, 28, 422 

-. 2009, ApJ, 697, 1071 

Baldini, L., et al. 2007, in American Institute of Physics Conference Series, Vol. 921, The 
First GLAST Symposium, ed. S. Ritz, P. Michelson, & C. A. Meegan, 190-204 

Breiman, L. 1996, Machine Learning, 24, 123, 10.1007/BF00058655 

Breiman, L., Friedman, J., Olshen, R., & Stone, C. 1984, Classification and Regression Trees 
(New York: Chapman & Hall), 358 

D’Agostini, G. 1995, Nuclear Instruments and Methods in Physics Research Section A: 
Accelerators, Spectrometers, Detectors and Associated Equipment, 362, 487 

Efron, B., & Tibshirani, R. 1993, An Introduction to the Bootstrap (Boca Raton, FL: 
Chapman & Hall/CRC) 

Esposito, J. A., et al. 1999, ApJS, 123, 203 

Finlay, C. C., et al. 2010, Geophysical Journal International, 183, 1216 

Fruhwirth, R. 1987, Nucl. Instrum. Meth., A262, 444 

Grove, J. E., & Johnson, N. 2010, Proc. SPIE, 7732 

Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, MNRAS, 369, 655 

Kalman, R. E. 1960, Transactions of the ASME-Journal of Basic Engineering, 82, 35 

King, I. 1962, AJ, 67, 471 

Kirsch, M. G. F., et al. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 5488, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, ed. G. Hasinger & M. J. L. Turner, 103-114 

Lande, J., et al. 2012, arXiv:1207.0027 

Lee, H., et al. 2011, ApJ, 731, 126 

Lott, B., et al. 2006, Nuclear Instruments and Methods in Physics Research Section A: 
Accelerators, Spectrometers, Detectors and Associated Equipment, 560, 395 

Mattox, J. R., et al. 1996, ApJ, 461, 396 



- 170 - 

Mizuno, T., Kamae, T., Godfrey, G., Handa, T., Thompson, D. J., Lauben, D., Fukazawa, 
Y., & Ozaki, M. 2004, ApJ, 614, 1113 

Moiseev, A. A., Hartman, R. C., Orrnes, J. F., Thompson, D. J., Amato, M. J., Johnson, 
T. E., Segal, K. N., & Sheppard, D. A. 2007, Astroparticlc Physics, 27, 339 

Moskalenko, I. V., & Porter, T. A. 2007, ApJ, 670, 1467 

Nolan, P. L., et al. 2012, ApJS, 199, 31 

Pelassa, V., Preece, R., Piron, F., Omodei, N., Guiriec, S., Fermi LAT, f. t., & GBM 
collaborations. 2010, arXiv: 1002.2617 

Rando, R., & the Fermi LAT Collaboration. 2009, arXiv:0907.0626 

Read, A. M., Rosen, S. R., Saxton, R. D., & Ramirez, J. 2011, A&A, 534, A34 

Rochester, L., Usher, T., Johnson, R. P., & Atwood, B. 2010, arXiv:1001.5005 

Smart, D. F., & Shea, M. A. 2005, Advances in Space Research, 36, 2012 

Smith, D. A., et al. 2008, A&A, 492, 923 

Student. 1908, Biometrika, 6, 1 

Walt, M. 2005, Introduction to Geomagnetically Trapped Radiation, Cambridge Atmo- 
spheric and Space Science Series (Cambridge LIniversity Press) 

Weltevrede, P., et al. 2010, Publications of the Astronomical Society of Australia, 27, 64 


This preprint was prepared with the AAS P-TLX macros v5.2. 



